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COMPUTER PROGRAM FOR STIRLING ENGINE PERFORMANCE CALCULATIONS 

Roy C. Tew, Jr. 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

I. ABSTRACT 

To support the development of the Stirling engine as a possible alterna- 
tive to the automobile spark-ignition engine, the thermodynamic characteris- 
tics of the Stirling engine were analyzed and modeled on a computer. The com- 
puter model is documented. The documentation includes a user's manual, sym- 
bols list, a test case, comparison of model predictions with test results, and 
a description of the analytical equations used in the model. 


II. INTRODUCTION 

The Stirling engine is being developed as a possible alternative to the 
spark-ignition engine under the Department of Energy's Stirling Engine Highway 
Vehicle Systems Program. NASA Lewis Research Center has project management 
responsibility for the program. 

A Stirling engine performance model has been developed at Lewis to support 
both the project management activities and the Stirling engine test program at 
Lewis. An early version of the model, published in reference 1, assumed fixed 
heater and cooler tube temperatures. The model was then expanded to include 
the coolant side of the cooler and used to make predictions for comparison 
with the single cylinder GPU-3 Stirling engine test results (ref. 2). More 
recently, variable specific heats, appendix gap pumping losses and adiabatic 
connecting ducts have been included in the model, and it has been used to sim- 
ulate one of United Stirling of Sweden's P-40 (approx 40 kW) engines. This 
engine, which has four cylinders and double-acting pistons, is now being 
tested at Lewis. Some of the test results are reported in reference 3; this 
reference also compares a few of the P-40 model predictions with some of the 
test results. Additional test results and a description of the test facility 
are reported in reference 4. 

This model predicts engine performance for a given set of engine operating 
conditions (i.e., mean pressure, boundary temperatures, and engine speed). 

One of the four engine working spaces is modeled, and the resultant power is 
multiplied by four (controls models such as the one documented in reference 5 
require modeling all four working spaces). The working space model includes 
two pistons, the piston swept volumes - the expansion and compression spaces, 
three heat exchangers - heater, regenerator and cooler, and four connecting 
ducts. The pistons are positioned as functions of time according to the spec- 
ified frequency. The working space is divided into appropriately sized con- 
trol volumes for analysis. Flow resistances and heat transfer coefficients 
are calculated for each control volume at each time step over the engine cy- 
cle. Within each gas volume the continuity and energy equations are integrat- 
ed with respect to time; a simplified momentum equation (pressure drop is a 



function of a friction factor and flow rate) and an equation of state are also 
used in the calculations. 

This report documents the current version of the Lewis Stirling engine 
performance model. A user's manual, symbols list, a test case and comparison 
of model predictions with test results for the P-40 engine are included in the 
documentation. 


III. MODEL DESCRIPTION 

The United Stirling P-40 engine, for which the test case and other model 
predictions were generated, is shown schematically in figure 1. The model 
simulates the thermodynamics of one of the four engine working spaces. The 
engine parts whose dimensions define the working spaces of the engine are: 

(1) the four pistons and four cylinders, connected in a square-four 
arrangement, as shown in the lower right corner schematic of Figure 1. 

(2) the circular array of heater tubes, the heater head, which connect 
the hot ends of the cylinders (expansion spaces) and the regenerators. 

(3) the eight regenerators (two per cylinder). 

(4) the eight coolers (two per cylinder) which connect the regenerators 
with the cold ends of the cylinders (compression spaces). 

(5) the four transition regions or connecting ducts per working space 
(expansion space-heater, heater regenerator, etc.). 

The hot expansion volume over one piston (part of the blackened area in 
the lower right corner schematic of fig. 1) is connected via one quadrant of 
the circular heater tube array, two regenerators and two coolers to the cold 
compression space volume beneath an adjacent piston; this constitutes one of 
the four working spaces. The model, assuming the four working spaces 
contribute equal amounts of power, multiplies the power predicted for the one 
simulated working space by four. 

The simulated working space was divided into control volumes as shown in 
Figure 2 for the test case. The model provides for one control volume each 
for the expansion space, compression space and the four connecting ducts (or, 
optionally, the connecting duct volumes may be lumped with the adjacent 
control volumes - thus neglecting the loss due to the adiabatic nature of the 
control volumes). For the test case, 3 heater, 5 regenerator, and 3 cooler 
control volumes were used. However, the heater and cooler may be divided into 
any number of equal sized control volumes. The regenerator may be divided 
into any odd number of equal-sized control volumes (the regenerator matrix 
temperature convergence method has been checked out only for an odd number). 

In addition, two (optional) isothermal appendix gap control volumes, one 
adjacent to the expansion space and one adjacent to the compression space, are 
available to evaluate appendix gap pumping losses. For all predictions 
discussed in this report, the 17 non-isothermal plus two isothermal appendix 
gap control volumes shown in Figure 2 were used in the model. 

The basic computer model equations are applied to each of the control 
volumes. The temperatures, masses, heat transfer coefficients, flow rates, 
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etc. for each of the control volumes and interfaces (except the appendix gap 
volumes and interfaces) are represented by dimensioned variable names in the 
computer model. The numbering procedure used for control volume and interface 
variable names is defined in figure 2. The circled numbers in figure 2 
correspond to control volume variables. The numbers with solid arrows 
correspond to interface variables. Appendix gap interfaces are labeled with 
numbers 0 and 17, respectively (dashed arrows); however, appendix gap volume 
and interface variables are represented by unique nondimensioned names in the 
computer model. 

The model has recently been generalized to allow changing the number of 
heater, regenerator or cooler control volumes by resetting the appropriate 
parameters. Preliminary results of changing the number of control volumes are 
summarized in table I. It is seen that increasing the number of control vol- 
umes increases the value of the predicted power. It also increases the value 
of regenerator effectiveness and the required computer time (regenerator ef- 
fectiveness, as used in this model, is defined in table IX). 

Increasing the number of control volumes should increase the accuracy of 
the model, since variables which change continuously along the working space 
are being approximated by lumped parameters which change discontinuously from 
one control volume to another. However, increasing the number of control vol- 
umes costs additional computing time. Also the model already overpredicts 
power and efficiency for the P-40 engine with the control volume configuration 
of figure 2. Additional runs are needed to define how the number of control 
volumes affects the trade-off between computing time and accuracy. For the 
purpose of this documentation the 17 nonisothermal plus 2 isothermal control 
volumes, as shown in figure 2, were used. 

The required engine operating conditions which must be input to the model 
are - heater tube outside wall temperatures (the combustor is not modeled), 
expansion and compression space inside wall temperature’s, cooling water inlet 
temperature, cooling water flow rate, engine speed, and mean pressure. The 
cooler tube inside wall temperature is solved for by iteration but is constant 
for any one cycle. The only wall temperatures which are allowed to vary dur- 
ing a cycle are the regenerator matrix temperatures; a technique for speeding 
up the convergence of these temperatures was used to get a solution in a rea- 
sonable amount of computing time. Cylinder and regenerator housing tempera- 
tures for conduction calculations can either be inputs or can be calculated 
from heater and cooler input temperatures. 

Losses due to imperfect heat transfer and appendix gap pumping losses are 
an integral part of the cycle calculations. The appendix gap pumping calcula- 
tions assume isothermal appendix gaps as in reference 6. A cold appendix gap 
is included for the sake of generality; however, its volume is very small and 
its effect is negligible for the P-40 engine. Heat conduction and piston 
shuttle losses are calculated and are accounted for in the efficiency 
calculations. 

The pressure drop and heat transfer calculations are based on correlations 
taken from Kays and London (ref. 7). The pressure drop calculations are based 
on a simplified momentum equation which neglects gas inertia. 
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Pressure drop calculations are also decoupled from the basic thermodynamic 
calculations for the working space to neglect pressure wave dynamics. A more 
rigorous modeling of pressure drop, accounting for pressure wave dynamics, 
would require a much smaller time step for stable calculations (with the ex- 
plicit, one iteration per step, numerical integration used in this model). 

In the early version of the model reported in reference 1, one pass, con- 
sisting of about 25 engine cycles, was made through the cycle calculations. 

In the model documented here two separate passes, using 25 engine cycles each, 
are usually made through the cycle calculations. The optional second pass was 
added to improve the modeling of the effect of pressure drop on engine 
performance. 

Calculated power loss due to pressure drop is about the same whether one 
or two passes are made. However, in the second pass calculations, the effect 
of pressure drop on heat transfer to and from the engine is more accurately 
modeled; the net effect on predicted performance is to increase the basic 
power (power before pressure drop loss) and efficiency of the engine. More 
details of the method used to account for the effect of pressure drop on en- 
gine performance are discussed in appendix E. 

For hydrogen at design P-40 conditions the effect of the second pass is to 
increase predicted brake power by about 1.2 kW (3.5 percent). For helium at 
design P-40 conditions a more significant increase, about 1.9 kW (9.7 percent) 
is found. As indicated above, the significant change is in the basic power 
(before pressure drop loss) and heat transfer; the power loss due to pressure 
drop is essentially unchanged. Since the model overpredicted power with one 
pass, the above changes increase the errors in predicted power. The shape of 
the revised predicted curve (power as function of speed at constant pressure) 
does however approximate more closely the shape of the experimental curve. 

Real or ideal gas equations of state can be used for pure hydrogen or he- 
lium working gas. Only the ideal gas equation of state can be used for a mix- 
ture of hydrogen and carbon dioxide. Working gas thermal conductivity, vis- 

cosity and specific heats are functions of gas temperature. 

Current computing time for the model is about 2.5 minutes for 50 cycles 
on an IBM 370 or 3 seconds per cycle^ This is based on 500 iterations per 
engine cycle or a time step of OXIO - ^ seconds when the engine frequency is 
66.7 Hz (4000 rpm). (1000 iterations per cycle were required to give 
satisfactory accuracy when trapezoidal integration was used for the work 
integration as in the model of ref. 1; it was shown there that when the number 

of iterations per cycle was reduced from 1000 to 200, the error in the 

prediction of both power and efficiency approached 10 percent; numerical 
stability was, however, maintained. It was then found that by switching to 
the more accurate Simpson rule integration, the number of iterations required 
for good accuracy decreased from 1000 to 500.) 

The analytical model upon which the computer program is based is discussed 
in appendix A. The working gas temperature differential equation used in the 
model is derived in appendix B. The method used to numerically integrate the 
decoupled gas temperature differential equation is explained in appendix C. 

The calculation of expansion and compression space heat transfer coefficients 
is discussed in appendix D. The simplifications made in the general form of 
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the one dimensional conservation of momentum equation and the decoupled pres- 
sure drop calculations are discussed in appendix E. The symbols used in the 
FORTRAN source programs and the input and output datasets are defined in ap- 
pendix F. Predictions of the model are compared with P-40 engine data in 
appendix G. 


IV. USERS' MANUAL 
A. Overall Simulation Structure 

The overall simulation structure is shown in figure 3. The computer model 
consists of a main program, MAIN, and five subroutines - ROMBC, HEATX, XDEL, 
CNDCT, and CYCL. 

In program MAIN, a data statement specifies the number of heater, regener- 
ator, and cooler control volumes and the number of time steps per cycle. Di- 
mensions for all control volume variables are specified in MAIN; instructions 
for setting the dimensions are given in comment statements in MAIN. MAIN com- 
municates only with subroutine ROMBC. 

ROMBC reads in the basic engine parameters and uses them to calculate con- 
trol volume geometry; it also reads in engine operating conditions, option 
switches (indexes), and multiplying factors. ROMBC initializes variables and 
steps time and crank angle; at each new crank angle it recalculates the varia- 
ble volumes and calls subroutine HEATX to update the working space heat and 
mass transfer calculations. ROMBC also integrates to determine work, stores 
working space variables for plotting, and averages working space variables 
over the cycle; instantaneous values of working space variable are also writ- 
ten out during each cycle (optional). Regenerator matrix temperature correc- 
tions, to speed up convergence, and cooler tube temperature corrections are 
made in ROMBC at the end of specified cycles. Subroutine CYCL is called at 
the end of each engine cycle to make summary calculations for the cycle. When 
predictions are completed for one set of operating conditions and ROMBC does 
not succeed in reading in a new set of input data, execution returns to MAIN 
for program termination. 

Subroutine HEATX updates pressure, heat transfer, gas temperatures, regen- 
erator matrix temperatures, gas flow rates and sums heat transfers over the 
cycle for use in energy balance and efficiency calculations. HEATX also calls 
subroutine XDEL to make a new calculation of engine pressure drop loss and 
subroutine CNDCT to calculate heat conduction losses. 

Subroutine XDEL calculates pressure drop for tube and wire screen fric- 
tion, tube 45, 90, and 180 degree turns, flow path contractions and expan- 
sions; subroutine calling arguments specify the type of pressure drop to be 
calculated and the flow geometry. 

Subroutine CNDCT calculates conduction losses through the cylinder 
housing, piston and the regenerator housing and also shuttle losses. 

Once per cycle, subroutine CYCL calculates the net heat into the engine by 
adding conduction, shuttle losses, etc., to the heat transferred into the 
working space over the previous cycle. Mechanical friction loss is calculat- 
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ed. Then net heat out is calculated by adding conduction, shuttle, appendix 
gap pumping and mechanical friction losses to the heat transferred out of the 
working space over the previous cycle. CYCL writes out summary results at the 
end of each engine cycle (optional). After the last engine cycle, CYCL cal- 
culates, auxiliary losses and brake power and efficiency; it then outputs an 
overall summary of operating conditions and performance results. 

The input data is read into ROMBC. The output data is written from either 
CYCL or ROMBC. The form of the input and output data will be discussed in the 
following two sections. 


B. Program Setup 

Array dimensions for all control volume variables are specified in the 
main program, MAIN. Several indexes which affect the choice of these dimen- 
sions are set in a data statement in MAIN; NH, NR, and NC specify the number 
of control volumes alotted to the heater, regenerator and cooler, respective- 
ly. The index, ISCD, is set equal to 1 to use separate control volumes for 
the following connecting ducts: expansion space-heater, heater-regenerator, 
regenerator-cooler, and cooler-compression space. If ISCD = 0 then the con- 
necting duct volumes are lumped with adjacent control volumes. The index, 
NITPC, specifies the number of time steps per engine cycle (normally = 500). 

The engine geometry is defined by reading the engine parameters into sub- 
routine ROMBC. For the test case, the P-40 engine parameters shown in 
table II were read into ROMBC via NAMELIST/ENGINE/. For convenience, the 
engine parameters of table II are defined in table III. To set the model up 
for another engine would require changing these engine parameters, the 
variable volume equations in subroutine ROMBC, and the mechanical and 
auxiliary loss equations in subroutine CYCL; the function definition, WINT, 
used in the work integration in ROMBC, would also need changing to be 
consistent with new variable volume equations. Also, the calls to the 
pressure drop subroutine XDEL from subroutine HEATX should be checked to see 
if the types of pressure drop calculations specified are appropriate for the 
new engine. 

The model option switches and multiplying factors and the engine operating 
conditions are also defined by reading the appropriate parameters into ROMBC. 
For the test case these parameters, as shown in table IV, were read in via 
NAMELIST /STRLNG/ and NAMELIST/INDATA/. The parameters of table IV are defined 
in tables V and VI. 

Multiple runs for a given engine can be made by adding sets of input data 
(table IV data), sequentially. After a run is complete the program tries to 
read a new set of input data; if another set of data is not found, the run 
terminates. 

The model options and multiplying factors shown in table V are discussed 
below: 

The parameter REALGS is set equal to 1 to use a real gas equation of state 
or equal to 0 to use an ideal gas equation. FACT1 and FACT2 are empirical 
factors used in the procedure for speeding up convergence of regenerator ma- 
trix temperatures. The current values, 0.4 and 10, respectively, have yielded 
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satisfactory results for all simulations attempted. The index, NOCYC, speci- 
fies the number of engine cycles to be calculated per pass; this is usually 
set at 25. However, there have been cases when a particular combination of 
operating conditions and engine parameters required as many as 40 cycles to 
get satisfactory convergence. (Convergence indicators are the percent errors 
in the engine and regenerator energy balances. These will be discussed later 
under Output - Test Case.) NSTRT specifies the cycle number at which the 
regenerator matrix and cooler temperature convergence procedures are turned 
on; this is usually set equal to one. The index, NOEND, specifies the cycle 
number at which the regenerator matrix and cooler temperature convergence 
procedures are turned off; this index is set at five less than NOCYC. For the 
last five cycles, the matrix energy equation alone determines regenerator 
matrix temperatures. This constitutes a check to see if the convergence 
procedure arrived at a temperature profile consistent with the basic matrix 
energy equation. (If it did not, then the percent error in the regenerator 
energy balance will increase during the last five cycles.) The index, 

MWGAS, = 2 to use hydrogen working gas or = 4 to use helium. RHCFAC, HHCFAC, 
and CHCFAC are multiplying factors for regenerator, heater and cooler heat 
transfer coefficients, respectively, for use in sensitivity studies. Set in- 
dex IPCV = 0 to make a second pass through the calculations or = 1 to elimi- 
nate the second pass. The first pass calculations include a correction of 
engine power and an approximate correction of heat into and out of the engine 
for the effect of pressure drop; the second pass provides a more accurate cal- 
culation of the effect of pressure drop on heat transfer and power. FMULT and 
FMULTR are overall pressure drop and regenerator pressure drop multiplying 
factors, respectively. Set index IMIX = 1 to use a mixture of hydrogen and 
carbon dioxide as the working gas; the mixture is defined by the volume frac- 
tion of hydrogen, VH2, which is set next after IMIX, as shown in table V. If 
IMIX = 1, then REALGS should be set equal to 0 and MWGAS should be set equal 
to 2. IMIX = 0 for pure hydrogen or helium. Set index IPUMP = 1 to include 
the piston-cylinder "appendix" gap pumping loss; IPUMP = 0 omits the pumping 
loss calculation. Set index ICOND = 1 to calculate cylinder and regenerator 
housing temperatures for conduction calculations from input hot end tempera- 
tures (TM(1) and TM(4)) and the coolant inlet temperature (TH20IN). If 
ICOND = 0, then the specified input values of cylinder and regenerator housing 
temperatures are used. The remaining indexes in table V are discussed in the 
next section. Output Options. 

The pressure drop subroutine, XDEL, is set up to calculate pressure drop 
due to tube friction, wire mesh friction, expansions, contractions, and 45, 

90, and 180 degree turns. The desired option is specified by setting the sub- 
routine calling argument, K, to the appropriate value of KTYPE as defined in 
comment statements in subroutine XDEL (KTYPE, in comment statements = K). 

C. Output Options 

The last five entries in table V define the output options. Three differ- 
ent sets of output data are available as shown in tables VII to IX. Table VII 
has been used primarily for debugging purposes and its shortened form, 
table VIII, has been used very little since the output of table IX was added. 

If the index, IOUT (in table V), is set equal to 1, then the value of the 
index, JIP, determines whether the data of table VII or VIII is output. If 
JIP = 0, then summary data for each cycle, as shown in table VII is output. 
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If JIP = 1, the shortened form, table VIII, is output (Symbols and units used 
in these two tables are defined in the symbols list in appendix F). If 
IOUT = 0, then neither of these two tables is produced. 

At the beginning of table VII, the input parameters and some of the cal- 
culated control volume parameters are output in NAMELIST write format. The 
remaining portion of the table contains summary data for each cycle. 

For example, during the first cycle, variables are written out at 
TIME = 0.0, the beginning of the cycle, and TIME = 0.0150 secs., the end of 
the cycle. The number of time steps between these variable printouts is 
specified by the index, IPRINT (table V); thus with 500 time steps per cycle, 

I PRINT = 500 yields the two lines of printout shown for all except the last 
5 cycles (per pass). For the last 5 cycles the number of time steps between 
variable printouts is changed to I PRINT /25 (= 20 for the test case); varia- 
tions over the cycle of gas temperatures, pressures, Reynolds Numbers and gas 
flow rates are shown in these expanded printouts for the last 5 cycles; the 
last cycle in table VII shows flowrates at intervals of 20 time steps (or 
25 printouts per cycle). 

Most of the quantities in table VII that are determined by summing or av- 
eraging over each time step in a given cycle (such as QIN, heat in per cylin- 
der per cycle, and QOUT, heat out per cylinder per cycle) are also printed out 
with definitions in table IX. There are some additional working space vari- 
able average and maximum values available in the output of table VII (and not 
in IX) that may be of interest. On the bottom half of the last page of 
table VII are average gas and metal temperatures (TGACYC, TMCYC), average and 
maximum heat transfer coefficients (HACYC, HMX), and, average and maximum heat 
fluxes (Q0AAVG, QOAMX); these values are shown for expansion and compression 
spaces, connecting ducts, and the hot and cold end control volumes for each 
heat exchanger. Expansion space values are at the extreme left and compres- 
sion space values are at the extreme right. Near the bottom of the page, brake 
power per cylinder, AUXPWR, and brake efficiency, AUXEFF, are found; these can 
be checked against the corresponding values in table IX for consistency. 

The index, ITMPS (table V), can be set equal to 1 to print variable tem- 
peratures at each time step for debugging purposes. The index, MAPLOT, = 1 to 
store cycle variables for plotting or = 0 to skip the storage procedure; the 
particular variables to be stored are specified in the FORTRAN coding of sub- 
routine ROMBC. 

Table IX, the Final Summary Printout, is written out after the final cycle 
for each set of input data; no provision is made for switching off this out- 
put. Table IX includes a summary of the engine operating conditions and 
predicted performance. The predicted performance includes an engine energy 
balance. For example the brake power plus total heat rate from the engine 
should equal the heat rate to the engine for a perfect energy balance. The 
values printed out in table IX show that the energy balance error for the test 
case, on a 4 cylinder basis, is: 
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BRAKE TOTAL HEAT RATE \ / 

POWER + FROM ENGINE , \ „ 1nn f 37.239 kw + 102.100 kW n 
TOTAL HEAT RATE TO ENGINE ~ ~ l 138.931 kW 

x 100 = 0.294 percent 

Brake and indicated engine efficiencies are also shown. Table IX also shows 
heat flows into and out of the engine separated into various parts. Regenera- 
tor heat flows, per cent error in regenerator energy balance and two measures 
of regenerator effectiveness are shown; the two measures of regenerator effec- 
tiveness are defined in the table (The two most important criteria in deter- 
ming whether the simulation has converged satisfactorily to a solution are the 
percent error in the engine energy balance and the percent error in the 
regenerator energy balance). Overall pressure drop is shown separated into 
parts by component. Maximum and minimum pressures and pressure ratios are 
shown for expansion and compression spaces. 

D. Program Execution 

Table X shows which Read/Write FORTRAN statement unit numbers, in sub- 
routines ROMBC and CYCL, were linked with the test case input and output 
data. For example, item 1 in table X indicates that a read statement in sub- 
routine ROMBC used unit number 4 to read in the input data of table II. 

The P-40 simulation test case was initiated by using system commands to: 

1. Link the READ unit #'s with the appropriate sets of input data as shown 
in table X. 

2. Link the WRITE unit fs with the desired locations for the tables of 
output data. 

3. Execute the main program, MAIN. 


V. OUTPUT— TEST CASE 

A sample run was made using the input data shown in tables II and IV. The 
engine parameters specified in table II, the variable volume equations in sub- 
routine ROMBC, and the mechanical and auxiliary loss equations in subroutine 
CYCL set the model up to simulate the United Stirling P-40 engine. The oper- 
ating conditions specified in table IV correspond to a NASA Lewis P-40 experi- 
mental run which approximated the design operating conditions of the engine 
(15 MPa (2175 lbf/in 2 ), 66.78 Hz (4000 rpm) ) . The effective outside heater 
tube temperature was estimated to be 930 K (1672° R); the cooling water inlet 
temperature was 323 K (580° R), and the cooling water flow rate was 
0.860 liter/sec (13.6 gal/min). 

The previously discussed tables of output data (VII to IX) were generated 
using the above test case operating conditions. Table IX shows that brake 
power and efficiency predicted for the P-40 at the specified operating 
conditions are 37.2 kW and 0.268, respectively. This efficiency does not ac- 
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count for external heat system (combustor) losses. Predicted heat flows, 
pressure drop losses and pressure ratios are also shown in table IX. 

The index, MAPLOT, was set equal to 1 (table IV) for the test case to 
store the cycle variables for plotting. A separate plotting program (not doc- 
umented in this report) using the IBM 370 graphics package was used to read in 
the stored data and make plots. The results are shown in figures 4 to 9. 
Expansion and compression space pressures are shown as a function of crank 
angle in figure 4. Expansion and compression space volumes are shown in 
figure 5. Expansion and compression space gas temperatures are shown in 
figures 6(a) and (b), respectively. They are also shown plotted to the same 
scale in figure 6(c). Gas flow rates out of the expansion space and into the 
compression space are shown in figure 7; flow is assumed positive from the 
expansion space toward the compression space (flow rates are also calculated 
at each interface between control volumes). Overall engine pressure drop 
(expansion space pressure minus compression space pressure) is shown in 
figure 8. P-V diagrams for expansion and compression space are shown 
superimposed in figure 9. 


VI. CONCLUDING REMARKS 

Testing of the United Stirling P-40 engine at NASA-Lewis has provided en- 
gine performance data for comparison with the predictions of this model. A 
comparison of a small portion of this data with predictions of this model is 
discussed in appendix G. 

There are several possibly significant losses which have not been included 
in the model and so have not been evaluated for the P-40 engine. Papers by 
Kangpil Lee (refs. 8 and 9) suggest that cyclic heat transfer or hysteresis 
loss due to heat exchange between working gas and cylinder walls may be sig- 
nificant for engines as small as the GPU-3. Also, no attempt has been made to 
evaluate losses due to non-uniform flow distribution or regenerator matrix 
"by-pass flow" (that is, flow near the regenerator housing that does not ex- 
change heat effectively with the matrix and thus degrades the effectiveness of 
the regenerator) 

The model's qualitative ability to predict variations in the state of the 
working gas over the cycle and to predict performance trends has been useful 
in helping to understand operation of the engine, to plan the experimental 
program, and to study sensitivity to various engine and working gas parame- 
ters. The model was recently used to generate a performance map for a 37 kW 
(50 hp) scaled down version of a United Stirling 67 kW (90 hp) engine design. 
The performance map was then input to a vehicle driving cycle code which was 
used to predict the fuel economy of a 37 kW Stirling engine powered vehicle. 

Predictions made with helium working gas tend to result in regenerator 
effectivenesses that appear too high unless an adjustment is made in the slope 
of the temperature variation across the regenerator control volumes. A typi- 
cal adjustment made for helium gas is to multiply the quantity, DTGASL, in 
subroutine HEATX by 0.95 (the model is set up to make this adjustment automat T 
ically when helium working gas is specified). 
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With relatively minor modification the model should be able to simulate 
Stirling cycle refrigerators or heat pumps. For the refrigerator application 
it would be necessary to define the working gas properties over a lower range 
of operating temperatures. 
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VII. APPENDIXES 


APPENDIX A: Analytical Model 

A rigorous simulation of the gas dynamics in the working space of a 
Stirling engine would require solution of a set of partial differential equa- 
tions. The simulation problem can be simplified by assuming the flow is es- 
sentially one-dimensional and dividing the working space into control volumes; 
a set of ordinary differential equations is then solved for each of the con- 
trol volumes; this is the approach is used here. 

Each of the control volumes shown in figure 2 is a special case of the 
generalized control volume shown in figure 10. The generalized control volume 
includes mass flow across two surfaces, heat transfer across one surface, and 
work interchange between gas and piston. Only the expansion and compression 
space control volumes are variable. 

The basic equations used to model the thermodynamics of the gas in each 
control volume are conservation of energy, mass, momentum and an equation of 
state. These equations are used to determine temperature, pressure and mass 
distributions within the working space at a particular time. 

The energy, mass and state equations, as written for the generalized con- 
trol volume shown in figure 10, are as follows, where three formulations of 
the equations of state are shown: 

Conservation of energy (for negligible change in kinetic energy across the 
control volume): 


< KC v T > 

= hA(T w - T) + 

( C W.T. - C WTJ 
P i 1 1 P 0 0 o' 

_ P dV 
v dt 

Rate of 

Rate of 

Rate of enthalpy 

Rate of 

change of 

heat transfer 

flow across 

work done 

i nternal 

across 

boundary of 

by gas in 

energy of 

boundary of 

control 

control 

control 

control 

volume 

volume 

volume 

volume 




( 1 ) 


Conservation of mass: 


4 " w. -w 

dt i o 


( 2 ) 


Equation of state: 


PV = MRT for ideal gas 

PV = MR[T + 0.02358 P] for hydrogen-real gas (3) 
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PV = MR[T + 0.01613 P] for helium-real gas 
where 

A heat-transfer area of control volume 

C p ,C v heat capacities at constant pressure and volume 
h heat-transfer coefficient 

M mass of gas in volume 

P pressure 

R gas constant 

T bulk or average temperature of gas in volume 

T-j,T 0 temperatures of gas flowing across surfaces i and o, 

respectively (in fig. 10) 

T w temperature of metal wall adjacent to heat-transfer area, A 

t time 

V volume 

W-j,W 0 flow rate across surfaces i and o, respectively 

(The real-gas equations of state were developed from data in ref. 10) 

Several assumptions are inherent in the use of these equations: 

(1) Flow is one dimensional. 

(2) Heat conduction through the gas and the regenerator matrix along the 
flow axis is neglected. The thermal conductivity of the regenerator matrix is 
assumed to be infinite in calculating the overall gas-to-matrix heat-transfer 
coefficient. 

(3) Kinetic energy can be neglected in the energy equation. 

(4) The time derivative term in the momentum equation is neglected (see 
appendix E). 

In appendix B it is shown that equations (1) and (2) and the ideal-gas 
equation of state can be used to derive the following differential equation: 


MC 


dT 
p dt 


= hA(T w - T) 


(W„ - W,)C„ r + (C M.T, - C„ III 
' 0 1 ' P ' p.j 1 1 P 0 0 o' 


+ V 


dP 

dt 


dC 

~ MT dt^ (4) 


The same result is obtained if either of the real-gas equations of state are 
used in the derivation. This equation says that the bulk or average gas tem- 
perature of a control volume is a function of the following four processes: 

(1) Heat transfer across the boundary from the wall 

(2) Gas flow across the boundary 

(3) Changes in pressure level 
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(4) Changes in working gas specific heat 

Equation (4) can be solved for the temperature derivative to get: 


dT _ hA , T 
it = w 


- T) 


("« 




W.T. 

Pi 1 1 


- C 

"RUT 


W T \ 

>o 0 °/ + V 
MC_ 


dP 

it 


T ^ 
^dt 


( 5 ) 


The approach used in numerically integrating equation (5) (suggested by 
Jefferies, ref. 1) was to decouple the four processes that contribute to the 
temperature change and solve for the temperature change due to each process 
separately. This approach allows a trade-off between computing time and 
accuracy of solution with little concern for numerical instabilities. The 
approach is suggested by representation of equation (5) in the following form: 


dT = dT + dT + dT + £T 

iT total itf due to dt due to dt due to dt due to 

heat transfer mixing pressure change in 

change specific 

heat 


where 


dT T dC p 

dt due to change TT" dt 

in specific heat p 


(5a) 


dT = V dP 

dt due to change MC dt 
in pressure p 


dT 

dt due to 
mixing 


(V w i) 

M 


T + 



1,T, - C W T n \ 

11 P 0 0 °) 



(5b) 


(5c) 


d T _ hA /t 

cTt due to " MC ' w 

heat transfer p 


T) 


(5d) 


In appendix C it is shown that equations (5a), (b), and (d) can be in- 
tegrated in closed form and that equation (5c) can be numerically integrated. 
When the results of appendix C are modified slightly to show just how they are 
used in the model, the resulting expressions are: 
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(5a 1 ) 
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(5b 1 ) 


(5c 1 ) 


(5d 1 ) 


where the superscripts t and t+At denote values of the variables at times t 
and t+At. The subscript S denotes the value of the temperature after it has 
been updated for the effect of change in the specific heat. The subscript SP 
denotes the value of the temperature after it has been updated for the effect 
of change in specific heat and pressure. The subscript SPM denotes the value 
of the temperature after it has been updated for the effects of change in 
specific heat, pressure and mixing. No subscript (as on the left side of 
eq. ( 5d 1 ) ) denotes the value of the temperature after it has been updated for 
all four effects-change in specific heat, pressure, mixing, and heat transfer 
to or from the metal. 

Discussion of Equations in Order of Calculation Procedure 

The equations considered so far have been derived and discussed with ref- 
erence to the generalized control volume of figure 10. In the computer model 
these equations are applied to each of the control volumes shown in figure 2. 
Thus temperatures, masses, heat-transfer coefficients, flow rates, etc., are 
all subscripted with an index for the non-isothermal control volumes. The 
index varies from 1 to NCV for variables that are averages for the control 
volumes and from 1 to NCV-1 for variables defined at the interfaces between 
control volumes (as shown in fig. 2). The equations discussed in this section 
include these indexes. The discussion of the equations follows the steps shown 
in the outline of calculation procedure in figure 11. 

Update Time and Crank Angle, ROMBC (Step 1, Fig. 11) : 

Time is an independent variable input to the computer model. The assump- 
tion of constant frequency means that both crank angle and time are updated by 
fixed steps at the beginning of each iteration. There are NITPC (= 500, 
usually) fixed time and crank angle steps per engine cycle. 
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Compression and Expansion Space Volumes, ROMBC (Step 2, Fig. 11) : 

Equations for expansion and compression space volumes as a function of 
crankshaft angle are as follows: 


V 


e 




cos a) + L 


( L - 


V 


c 





^ 1- - (r si " «) 2 )] + V e, clearance (6) 

- L ^- - Al !• - (f si " (« + f)) J 

+ V c clearance, 


where V e , clearance includes the hot appendix gap volume and V c , clearance 
includes the cold appendix gap volume - only if the appendix gap pumping 
losses are not calculated (that is, only if the index, IPUMP = 0) 

where 

V e expansion-space volume 

V c compression-space volume 

Ap piston cross-sectional area 

Ap r piston minus piston rod cross-sectional area 

r crank radius 

L piston rod length 

o crank angle 


Thermal Conductivity and Viscosity Equations (Step 3) : 

The equations used in subroutine ROMBC to calculate gas thermal conductiv- 
ity and viscosity (for hydrogen, helium and carbon dioxide) assume both 
quantities vary linearly with temperature. The equations are derived from 
data in reference 11. The mixture equations used to determine the con- 
ductivity and viscosity for a mixture of hydrogen and carbon dioxide are based 
on information in references 12 and 13. 


Pressure (Step 4) : 

The pressure P is calculated by 


p t + At 


» R 


NCV 

£ 

i=i 




ri 


NCv 

^ V t + Atpt+A t 
1=1 


if IPUMP = o 
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M* T, + 
hgp ngp 


if I PUMP = 1 


(MjTj) 


= R 


+ M t t 
cgp'cgp 


V. Fi +At + 
hgp 1 


/..t+Atpt+At, .. pt+At 
lv I r I ' cgp nCV 


I 


where 

P pressure at center of regenerator 

R gas constant 

M control volume mass 

T control volume average temperature 

V control volume 

F ratio of pressure in control volume to pressure at center of 

regenerator 

I index denoting which of control volumes 1 through NCV is under 

consideration 

hgp subscript denoting isothermal hot gap control volume 
cgp subscript denoting isothermal cold gap control volume 

Equations (8) are obtained by first writing the ideal-gas equation for the 
Ith control volume 


Pi Vj = (F I P)V I = MjRTj 


where 

Pj pressure in Ith control volume 

then summing over each of the control volumes under consideration and solving 
for P, the pressure at the center of the regenerator. 

Equations (8) indicate that if appendix gap pumping losses are included in 
the model (by setting IPUMP = 1), then the summation is over each of the 1 
through NCV plus the two isothermal appendix gap control volumes. If appendix 
gap pumping losses are not included (IPUMP = 0) then the appendix gap volumes 
are lumped with expansion and compression space clearance volumes and the sum- 
mation is only over the 1 through NCV control volumes. The appendix gap 
pumping loss model used is based on the one in reference 6. 

If the real-gas equation of state for hydrogen is used and the same pro- 
cedure is followed, the result is: 



if I PUMP = 0 


NCV 


,t+At 


= R 


1=1 




NCV 

E 

1=1 


p t+At v t+At _ 0.02358 R 


Wv 


e m 


1=1 


(9) 


NCV 


M 


Thgp + £ (Mfr‘) ♦ M‘ gp T cgp 


hgp hgp 

HUJ 


= R 

r /rt+Atw . / C t + At,,t+At\ , ct + At„ \ 

[( F l %P E ( F I V I ) +F NCV V cgpj 


- 0.02358 R (F^Xgp + E (^‘M) + C^cgp)] 1 ' 


if I PUMP = 1 


Equations (8) and (9) are both included in the model. Also an equivalent 
real-gas equation for helium is included. An index, REALGS (table IV), speci- 
fies whether a real or ideal equation is to be used. 


The variable, F, represents an array of pressure ratios with different 
values for each control volume at each time step over an engine cycle. F is 
used above to evaluate the effect of the decoupled pressure drop calculations 
on pressure level at the center of the regenerator. F is defined as follows: 


During the first pass through the calculations (NOCYC cycles), each ele- 
ment of F = 1. At the end of the first pass, pressure drop information for 
each control volume at each time step is stored in F; each element represents 
the ratio of pressure at a particular control volume and increment of time to 
pressure at the center of the regenerator. Thus the array appears as: 


F l,l F l,2 F 1,NRC-1 1 F 1,NRC+1 — F 1,NCV 

F 2,l F 2,2 F 2, NCR-1 1 F 2,NRC+1 •“* F 2,NCV 

NITPC 

1 ROWS 

F-j k = . (time 

increases) 


f NITPC,1 f NITPC,2 F NITPC,NRC-1 1 f NITPC,NRC+1--* f NITPC,NCV 

NCV COLUMNS 1 
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where 

NCV = 17 

NRC = 9 

time steps per engine cycle, NITPC = 500 for the sample runs. 

This array of pressure drop information is used in the second pass (NOCYC 
additional cycles). 

The second pass, using the array F, was incorporated to get a more accu- 
rate evaluation of the effect of the decoupled pressure drop on the thermody- 
namic calculations. The method is discussed in more detail in appendix E. 

Update Gas Specific Heats (Step 5) : 

The equations used to calculate gas specific heats in subroutine HEATX 
assume a quadratic variation with temperature (except for helium whose specif- 
ic heat is essentially constant over the temperature range of interest). The 
equations are derived from data in reference 11. 

Update Temperatures for Effect of Change in Specific Heat (Step 6) : 

Introducing the control volume index, I, into equation (5a‘), the form of 
the equation used to correct gas temperatures for the effect of changes in 
specific heat is: 


.t+At T t P I 
I,S “I pt+At 

Pi 


Update Temperatures for Effect of Change in Pressure (Step 7): 


( 10 ) 


t+At 

I,sp 



pt+Atpt+At 
ptpt 



(ID 


This equation is commonly used to relate temperature and pressure for an adia- 
batic fixed-mass process. 

Mass Distribution (Step 8) : 

On the first pass through the calculations the mass distribution is deter- 
mined by assuming that the mass redistributes itself in accordance with the 
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new volumes and temperatures in such a way that pressure is uniform throughout 
the working space. The pressure, P, throughout the working space is derived 
from the perfect-gas law as follows: 

The perfect-gas law for the Ith control volume can be written 


Summing over the NCV control volumes (and the two isothermal appendix gaps 
control volumes if IPUMP = 1) 


NCV 

NCV 

f i v i 
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R - NCV Z ,, 
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Now substituting for P/R into the perfect-gas equation for the Ith control 
volume gives 


M. 4. F t V t 
„ total I I 

M i = m t — 


z 

1=1 


f i v i 


IPUMP = 0 


The form of this equation used in the model to calculate the working gas mass 
in each of the NCV control volumes is: 
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p t + At..t + At 

■"I V I 



(where Tj + ^ represents T updated for changes in specific heat and pressure 


but not for mixing and heat transfer). 

The preceding equation calculates the new mass distribution for the case 
of a perfect gas. The following equation, which can be derived in the same 
manner, is used to approximate the real properties of hydrogen for the case of 
no appendix gap pumping loss: 


M 


t+At 

I 


pt+At,.t + At 

r I V I 



(13) 


A similar equation that approximates the real properties of helium is included 
in the model. 

When appendix gap pumping losses are included (IPUMP = 1) then the summa- 
tions used in deriving the equivalents of equations (12) and (13) are over the 
1 through NCV plus the two isothermal appendix gap control volumes. 


Flow Rates (Step 


Once the new mass distribution is known, the new flow rates are calculated 
from the old and new mass distributions according to 


o o 


if IPUMP = 1 


and 


(14) 


21 



Wj = - 1 - - ^ ■ I + W I _ 1 I = l.NCV if IPUMP = 1 

or I = 1,10-1 if IPUMP = 0 

where 

M 0 = mass in hot appendix gap if IPUMP = 1 (not used if IPUMP = 0) 

W 0 = FHGP, flow from hot gap to expansion space, if IPUMP = 1 (not used if 
IPUMP = 0 

W NCV = FCGP, flow to cold qap from compression space, if IPUMP = 1 (not used 
if IPUMP = 0) 

Wj = is the flow rate at the Ith interface between control volumes. 

Update Temperatures in Each Control Volume for Effect of Gas Flow Between 
Control Volumes (Step lOj ; 

The following equation (modification of eq. (5c')) was used to update 
temperature for the mixing effect following gas flow between control volumes: 


T t+At 

'i 
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^1 T t+At + 

t+At I,SP 
I 



1-1 


t+At t+At r t+At, ,t+At 

' 1-1 9 I -1 " C Pl ^1 


t+At^ 


M t+At r t+At 
"i L p 


(15) 

I = 1,NCV (IPUMP = 1) 


(where the temperature Tj sp has already been updated for specific heat 
change and pressure change) . 

The temperature of the fluid flowing across the interface has been given a 
new variable name, e, to better distinguish it from the average control vol- 
ume temperature, T, and to keep the subscripts as simple as possible. The 
procedure used to update the temperature, e, for each interface is now 
defined. 

The temperature of the fluid flowing across a control volume boundary is 
just the bulk temperature of the control volume from which the fluid came - 
for flow from the expansion space, heater, cooler, or compression space con- 
trol volumes or hot and cold appendix gap volumes. This is a reasonable as- 
sumption for these volumes since the actual temperature gradient across each 
is expected to be relatively small. In a five-control-volume regenerator, 
however, the temperature gradient is not small. One option would be to in- 
crease the number of control volumes in the regenerator. However, to save 
computing time, an alternative approach was used. It was assumed that a tem- 
perature gradient existed across each volume in the regenerator. The magni- 
tude of the gradient was assumed to be equal to the corresponding regenerator 
metal gradient. 
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A schematic of a regenerator control volume is shown in figure 12(a). 

Flow across both interfaces is, for now, assumed to be in the direction shown 
(which is defined to be the positive flow direction). The cross-hatched area 
represents the portion of the fluid that will flow across interface I during 
the time step, At. The assumed temperature profile of the control volume is 
characterized in figure 12(b). The vertical dashed line in figure 12(b) 
defines the temperature at the left boundary of the fluid that will flow 
across interface I during At. If Tj is defined as the average temperature 
of control volume I and aTj equals one-half the change in temperature across 
the control volume, then Tj - aTj is the temperature of the fluid at 
interface I and 


W. At 

T i _ aT i + 2 At i 


is the temperature of fluid at the vertical dashed line, (figure 2 shows the 
numbering methods used for the control volumes and the interfaces between con- 
trol volumes). 

The temperature of the fluid that flows across an interface during At is 
assumed to be equal to the average temperature of that fluid before it crosses 
the interface. The average temperature of the fluid in the cross-hatched area 
of figure 12(a) is then 


1 

2 




Wj At 



+ (Tj - ATj ) 



W T At 

aT I + AT I 


Therefore, for the flow directions shown in figure 12(a), the updated tempera- 
tures of the fluid that crosses the interfaces during At are 
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If the flow direction is reversed at both interfaces, then 


(16) 
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( 17 ) 
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Heat Transfer Coefficients (Step 11): 


The heat transfer coefficient calculations for heater and cooler are 
derived from figure 7-1, page 123 of reference 7; heat transfer coefficient 
calculations for the regenerator are derived from figure 7-9, page 130 of the 
same reference. The assumption used in calculating heat transfer coefficients 
for the expansion and compression spaces are discussed in appendix D. 


Update Temperature in Each Gas Control Volume for Effect of Heat Transfer 
Between Gas and Metal (and Determine Heat Transfer Between Gas and Metal) 
(Step 12jl 

This temperature update for control volumes 1 through NCV is accomplished 
by using the following equation (a modification of equation ( 5d ' ) ) : 


t+At _ T t+At / T t T t+At \ , 

I “ I,SPMV W, I " I,SPM/ L 1 
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h t+At A 
h I A I 

M t+At r t+At 

M T 

i p, 


I At 


1=1, NCV 


where Tu j is the wall temperature of the Ith control volume. Note that, no 
matter now large the heat-transfer coefficient, the gas temperature cannot 
change more than the aT between the wall and the gas. Thus this calculation 
cannot cause the solution to become unstable, but it can lead to significant 
inaccuracies if the time increment. At is made too large. 

The heat transferred between gas and metal is then calculated from: 


n t+At / X t+At T t+At \ M t+At r t+At 

Q i - -( T I - t i,spm) h i c p r 


1=1, NCV 


(18) 


so that heat transfer from gas to metal is defined to be positive. 

The appendix gap control volumes are assumed to be isothermal. Three 
steps are used to calculate the heat transfer between the cylinder wall and 
the appendix gap working gas that would be required to maintain constant gap 
temperature. 
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For the hot appendix gap: 


1. The change that would occur in appendix gap temperature due to pressure 
change if there were no heat exchange with the cylinder wall is: 


aT, 


hgp T hgp 


pt + At F t+At 

ptpt 


hgp 



F]_ is used since there is assumed to be no pressure drop between the hot gap 
and the expansion space. 


2. The net change that would occur in appendix gap temperature due to 
pressure change and mixing, if there were no heat exchange, is calculated by 
adding an additional "mixing" term to the above expression when there is flow 
from the expansion space to the appendix gap. That is: 
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where 
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hgp is a subscript denoting quantities within the hot gap control volume 

hgpj is a subscript denoting variable values at the flow interface between 

the expansion space and the appendix gap. 


3. The rate of heat exchange with the cylinder wall required to maintain 
an isothermal appendix gap is then calculated to be: 


hgp 


aT M t+At r 

AT hgp M hgp C p 
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A similar set of calculations is made for the cold appendix gap. 
Regenerator Metal Temperature (Step 13) : 

The equation used to update the metal temperatures in the regenerator con- 
trol volumes is 
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I=NR1,NRL 


(19) 


MjC 




where Qj is the rate of heat transfer between gas and metal. This is inte- 
grated numerically by setting 
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( 20 ) 


where 


Mj mass of metal in Ith volume 

C thermal capacitance of metal 

At time increment 

For most regenerators the thermal capacitance of the metal is so much 
larger than the thermal capacitance of the adjacent gas volume that an exces- 
sive number of engine cycles (from the point of view of computing time) are 
required for the metal temperatures to reach steady state. Therefore, it is 
necessary to apply a correction to the metal temperatures after each cycle to 
speed up convergence. The method used is discussed in a later section. 

Pressure-drop Calculations (Step 14) : 

Since the pressure-drop calculations have been decoupled from the heat and 
mass transfer calculations, pressure drop need not be re-calculated over every 
cycle. Pressure drops are re-calculated only every third cycle. Thus the 
indicated work calculation is corrected using the most recently calculated 
loss. 

A general form of the conservation of momentum equation for one- 


dimensional 

flow is: 
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where 


p density 

v velocity of flow 

f friction factor 

Df, hydraulic diameter 
P pressure 

t time 

x distance 

In appendix E it is shown that by combining the continuity and momentum 
equations and then neglecting the time derivative term in the resulting 
equation, the following equation results: 


vdv + 



dx + — = 0 
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( 22 ) 


This equation can be integrated over a length L for the special cases of 
adiabatic or isothermal flow processes (the two extremes). When the re- 
sulting adiabatic and isothermal expressions were applied to the P-40 heat 
exchangers (by setting index, K, appropriately in the call to subroutine 
XDEL from subroutine HEATX for heat exchanger pressure drop calculations), 
the contribution of the vdv term was negligible for the two extremes. By 
neglecting the vdv term, the expression for pressure drop is reduced to 


1 vl dx + dP _ 0 

2 D h QX p U 

or applying the differential equation (23) over a finite length L 
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where aP is the pressure drop over length L (using the adiabatic or isother- 
mal forms of the pressure drop equation with the vdv term retained requires 
an iterative solution procedure which increases computing time by about 
20 percent). 

A modification of this equation can also be used to account for the 
effect of expansions and contractions in flow area. The form of the modi- 
fied equation is: 



( 25 ) 
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It is applied at each area change in the flowpath between the expansion and 
compression spaces. At a particular point where an area change occurs, K is 
a function of the two areas and the direction of flow (since an expansion for 
one flow direction is a contraction when the flow reverses). The term, K ,is 
calculated in accordance with the procedure given in references 14 and 15. 
Values of K are also specified to account for pressure drop due to tube 
bends. 

The types of pressure drop calculations that can be made in subroutine 
XDEL are specified by the calling argument, K, and are defined in comment 
statements in the subroutine (KTYPE, in comment statements = K). 

For the heater and cooler control volumes the friction factor, f, is de- 
termined from equations based on figure 7-1, page 123 of reference 7. The 
friction factor for the regenerator is derived from figure 6.3-1, page 6-35 of 
reference 16. 

With the pressure level, P, known (assumed to be the pressure at the cen- 
ter of the regenerator) and the aP's across each of the control volumes calcu- 
lated, the pressures needed in the work calculations, P e and P c , can be 
calculated as follows: 


NRC-1 



NRC+1 


Near the end of the first pass the pressure drop information for each con- 
trol volume over a complete cycle is incorporated into the array of pressure 
ratios (discussed under step 4) for use in the second pass. 

Heat Conduction From Hot End to Cold End of Engine and Shuttle Loss 
(Step 15) : ~ 

Three separate paths were considered in the calculation of heat conduction 
losses from the hot end to the cold end of the engine: 

(1) Through each of the regenerators 

(2) Through the cylinder wall 

(3) Through the wall of the piston from the hot space to the cold space 
The effect of temperature on metal conductivity was accounted for. 
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The piston picks up heat from the cylinder at the hot end of its stroke 
and loses heat to the cylinder at the cold end of its stroke. This shuttle 
loss is calculated by using the following equation from reference 17: 


n _ KttDS 2 aT 

^shuttle SCI 


(26) 


where 

K thermal conductivity of gas 

D piston diameter 

S stroke 

aT temperature difference across displacer length 

C clearance between displacer and cylinder 

L displacer length 

The conduction and shuttle losses are calculated once per cycle. The cal- 
culations could be made just once per run except that the conduction through 
the piston is assumed to depend on the average gas temperature (The conduction 
through the piston is sufficiently small for the P-40 that a once per run cal- 
culation would yield very little error.) 

Sum Up Heat Transfers Between Gas and Metal for Each Component (Step 16) : 

The basic heat into the working space per cycle is the sum of the net heat 
transfer from metal to gas in the heater and expansion-space control volumes 
over the cycle. The basic heat out of the working space per cycle is the sum 
of the net heat transfer from the gas to the metal in cooler and compression- 
space control volumes per cycle. Since it is assumed that there are no losses 
from the regenerator matrix, the net heat transferred between gas and metal in 
the regenerator over a cycle should be zero. This net heat transfer in the 
regenerator over the cycle is the most convenient criterion for judging when 
convergence of regenerator metal temperatures has been achieved. 

The net heat into the engine is the basic heat (as defined above) plus 
conduction and shuttle losses. The net heat out of the engine is the basic 
heat out plus conduction, shuttle, appendix gap losses, mechanical losses and 
auxiliary power losses. Conduction, shuttle, appendix gap and mechanical 
losses (and any heat transfer out via the compression space) are assumed to 
pass into the cooling water but not through the cooler tubes (there are 
cooling water flow passages in contact with the cylinder). It is arbitrarily 
assumed that the auxiliary power requirement does not increase the heat load 
on the cooling water but is dissipated via convection and radiation to the 
surroundings. 
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to: 


Work Calculations (Step 17) : 

The indicated work, neglecting pressure drop loss, is calculated according 




W =f n P(d V + DV ) 
e c 


(27) 


The indicated work, accounting for pressure drop loss, is calculated according 
to: 

W =t/)(P e dV e + P C dV c ) (28) 

From the volume equations for the P-40 engine, (6) and (7) it is found that 


dV = A r sin a 
e p 


1 + 


j- COS a 


^1 sin 


da 


dV c - - A pr r si " (“ + f) 


T COS 


1 + 


(“ + f) 


^ 1 - (f sin(« 


da 


or, defining 


F(P,A,d) = PAr sin m 1 + 


then 


W = CD (P dV + P dV ) 
' e e c c' 


£ cos t 


^ - (f 5,n< ) 


= ^ F (Pe> A p'«) tF ( P c-V- a *?) 


da 


f (a) da 


30 



The above integration over a cycle was accomplished numerically using 
Simpson's rule integration, that is: 



(a)da 


= [f(« 0 ) + 4f (a x ) + f(a 2 )J 

A number of additional work calculations were made to separate the work 
loss due to pressure drop for each of the components and for the end effects. 

The chart in table XI shows how the various pressure and work para- 
meters were made equivalent to arrays to allow reducing the number of pro- 
gramming steps required for the calculations; this chart is included only as 
an aid in following the FORTRAN programming steps of subroutine ROMBC. 

Is Cycle Complete? (Step 18) : 

The number of iterations made during the current engine cycle is 
checked to see if the cycle is complete. If the cycle is not complete, then 
the model loops back to step 1 and another iteration is begun. 

Convergence Method for Regenerator Metal Temperatures (Step 19) : 

The correction to regenerator matrix temperatures between cycles to 
speed up convergence (suggested by Jefferies, ref. 1) is made as follows: 


aT t = 


TIME=NAt 

E 

TIME=0 


< T w,I " V 


1 - 


h I A I 

M I C P, 


At 


M, 


TIME=NAt 

E 

TIME=0 


I=NR1 ,NRL 


1 - 


Vi 

M I C p 


I At 


(29) 


where 

N number of iterations per cycle 

aTj weighted average difference between wall and gas temperature over 
cycle for Ith regenerator control volume 

Tj Ith gas temperature (instantaneous average over control volume) 

Tyj^j Ith wall temperature 


Then let 
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I-NR1.NRL 


AT j = aTj^qlq x FACT^ + aT j x FACT2 

where FACT1 = 0.4 and FACT2 = 10.0 are the factors that were found to work 
best when the method was originally developed. (These factors were not re- 
optimized for the P-40 engine). The final step in the correction is: 


t w,NR1,NEW = t w,NR 1,0LD ~ ( RC 1,1 x aT NR1 + RCi j2 x aT nr1+1 + .... + RC 1>NR x aT nrl ) 
t w,NR1+1,NEW = t w,NR1+1,0LD - ( r C2,1 x aT NR i + RC2,2 X aT NR i+i + + RC2.NR x aT nrl ; 


t w,NRL,NEW = T w,NRL,0LD ~ ( rc NR,1 x aT NR1 + rc NR,2 x aT NR1+1 + + rc NR,NR x aT nrl) 

where the coefficients RCj^ are calculated as follows: 


For k > i, RC.j ^ = 


(NR + 1 - k)NR 
NR 

E k 

k=l 


if i = 1 


= i x RC 


l,k 


if i ^ 1 


for k < i, RC^ = RC^* 

For NR=5 (that is, 5 regenerator control volumes) the coefficients generated 
by the above equations are: 


RC, 


i,k - 


5 4 2 1 

1 

3 3 3 3 

4 8 4- 2 

2 

3 3 3 3 

1 2 3 2 1 

2 4 8 4 

2 

3 3 3 3 

12 4 5 

1 

3 3 3 3 
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Caclulated Indicated Power and Efficiency (Step 20) : 

Indicated efficiency is defined to be the indicated work divided by the 
net heat into the engine (per cycle). 

Calculate Mechanical (Plus Leakage) Loss (Step 21) : 

The mechanical loss calculations for the engine are based on informa- 
tion obtained from United Stirling. 

The mechanical loss per engine (4 cylinders) is assumed to be: 


M.L. 


12 - 8 tt 


(P + 5) 
20 


where 

M.L. mechanical loss/cylinder in KW 

N engine speed 

Nd design engine speed 

P mean pressure in MPa 

This "mechanical loss" is also assumed to include loss due to leakage. A plot 
generated with this equation is shown in figure 13. 

Calculate Auxiliary Losses and Brake Power and Efficiency (Step 22) : 

A plot of the auxiliary power requirement is shown in figure 14. The only 
auxiliary power requirement assumed to change significantly with the mean 
pressure level is that of the combustion blower. The auxiliary power require- 
ment for mean pressures between 15 and 4 MPa is obtained by interpolating be- 
tween the two curves. The lower curve is assumed to define the minimum auxil- 
iary power requirement. 

Brake power is defined to be indicated power minus mechanical friction and 
auxiliary losses. Brake efficiency is defined to be the brake power divided 
by the net heat rate into the engine. The net heat rate into the engine is 
defined to be the net heat transfer from metal to gas in the heater and expan- 
sion space plus conduction and shuttle losses. 

Convergence Method for Cooler Tube Temperatures (Step 23) : 

The cooler tube temperature is a function of cooling water inlet tempera- 
ture and flow rate, and the rate of heat out through the cooler. Since the 
rate of heat out through the cooler is a function of cooler tube temperature, 
an iterative procedure is required to solve for cooler tube temperature. 

The cooler tube temperature is updated every third cycle during the same 
period that the regenerator matrix temperature convergence procedure is opera- 
tive. The procedure used is outlined as follows: 
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T H 2 0 " T H 2 0,IN + 

g H 2 0,0UT 

2C p,h 2 o^h 2 o 

where 


t h 2 o 

average coolant temperature 

t h 2 0’ in 

inlet coolant temperature 

q h 2 o, out 

rate of heat out through coolant 

c p,h 2 o 

coolant specific heat 

w h 2 o 

coolant flow rate 


Calculate water side heat transfer coefficient using the following two steps 
to incorporate a fouling factor: 


1. h 2 = 


12k H 2 0 

• 35D o N Re* 55 N pr - 333 


2. h = 


1 . 


-J— + 1.8 

h l 


where 

h heat transfer coefficient, incorporating a fouling factor 

k 

H^O thermal conductivity of coolant 

Nr 6 coolant Reynolds number 

Np r coolant Prandtl number 

D 0 cooler tube outside diameter 
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The fouling factor 1.8 has units sec-ft 2- 8 R/Btu (0.881 cm^- # K/w) 
Calculate water side and cooler tube thermal resistances: 


R H 2 0 “ hA 


1. 

HT N T 


12 LOG 


_0D 

3 ID 


TUBE 


2irL N T k 
e T ss 


where 


AHT 

n t 

JID= 

^SS 


J 0D 


water side heat transfer area per tube 
no. of cooler tubes 

effective heat transfer length of cooler tube 
cooler tube inside and outside diameters, respectively 
cooler tube thermal conductivity 


Then the cooler tube temperature is updated as follows: 


t new = t h 2 o + 1,287 E ~ 3 q h 2 o,out “( r h 2 o + R TUBE^ 


T NEW " 0,5T NEW + 0,5T 0LD 


where 


t NEW 

new tube temperature 

Told 

old tube temperature 

Qh^.out 

rate of heat out through cooler 

(i) 

engine frequency 


The compression space wall temperature, TMfjQy, is then set as follows: 


™ncv - t h 2 o,in + ^ho t h 2 o,in^ 

Have the Specified Number of Cycles (NOCYC) Been Completed? (Step 24) : 

A check is made to see if the specified number of cycles (NOCYC) has been 
completed. If not the model loops back to step 1. If, yes, then the proce- 
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dure continues to the next step, 25, provided the two pass option, IPCV = 0, 
was specified (If a one pass option was chosen, IPCV = 1, then the procedure 
jumps to step 27, skipping steps 25 and 26). 

Is This the Second Pass Through NOCYC Cycles? (Step 25) : 

A check is made to see if the the second pass was just completed. If not, 
the second pass is begun (step 26). If, yes, then the procedure continues to 
the final step, 27. 

Second Pass Calculations (Step 26) : 

Time is reset to zero. The pressure drop information from the first pass 
is used in making working space thermodynamic calculations (instead of using a 
uniform pressure throughout the working space for these calculations) when the 
procedure loops back to step 1 and begins the second pass iterations. 

Final Step (Step 27): 

When the second pass is completed, if IPCV was set equal to 0, or the 
first pass is completed, if IPCV was set equal to 1, then the summary of pre- 
dictions shown in table IX is written out. The model then attempts to read in 
a new set of operating conditions; if succesful, the entire calculation proce- 
dure of figure 11 is repeated; when no new operating conditions are found, the 
simulation is terminated. 
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APPENDIX B 

DERIVATION OF GAS TEMPERATURE DIFFERENTIAL EQUATIONS 


The basic gas volume equations used in the derivation are: 

it (MC T) = h A ( T - T) + (C W.T. - C W T ) - P 
dt ' v ' ' w ' ' i l p o o o' dt 


PV = MRT 

Expanding the first term of equation (1) gives: 


: vS +c v T f + HT itr - hfl < T „ - T > + < c p"i T i- 


V° T »> 


P dV 

v dt 


Differentiating equation (3) gives: 


mr £ + rt h" p ^ + v £ 


Letting R = C p - C v in the first and second terms of equation (B2) and 
solving for 


C T — 

V dt 


yields 


r T dM M , r r vdT , r T dM D dV dP 
V dt - M(C p - C v>dt + V dt " P ^t ■ V dt 


Substituting the right side of equation (B3) for the second term of equation 
(Bl) yields: 


X& +m(c p'X>& +c p t $ - p X- 


v f + MT w- 


= hA(T - T) + (C n W.T. - C n WTJ - P 
W ' p. 1 1 p_ o o' 
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dC 


MC 


£ + C p T $ - v % + MT *r - “(Tw - T > + < C p,“i T i' - C P "o T o> 


(B4) 


Using equation (2) to substitute for dM/dt in equation (B4) and also 
substituting 


dC dC 

dF = It- in {B4) 


then solving for MC p gives 


ht hp dC 

MC ttf = hA(T - T) + (W n - W. )C T + (C W.T. - C W T ) + V - MT -r^- (4) 
p dt w ' o i p i i p Q o o dt dt 

which is the equation used in the model to solve for gas temperature. 
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APPENDIX C 


INTEGRATION OF DECOUPLED TEMPERATURE EQUATIONS 
Temperature time derivative due to change in specific heat: 


-nr due to 
change in 
specific heat 


, dC 

_I E 

C p dt 



Integrating - 


LN T 


t+At 

= - LN C 
+ P 


t+At 

t 


rt+At 


LN 



C 

T t+At T t p 

~ pt+At 

P 


(5a) 


Temperature time Derivative Due to Pressure Change ; 
The equation 


dT _ V dP 

dt due to - MC dt 
change in P 
pressure 


(5b) 


can be integrated in closed form (if it is assumed that Cq is constant 
over the time increment) by solving the equation of state for V/M and sub- 
stituting in equation (5b). 


V RT 
PV = MRT » J = f- 

Substituting 
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dT _ RT dP 

It " PIT dt 

P 


dT _ R_ dP 
T = C p P 

Assuming Cp is constant over At and integrating - 



For the second pass it is assumed that P is the pressure in the cen- 
ter of the regenerator and a pressure ratio factor (ratio of pressure at the 
control volume of interest to pressure in the center of the regenerator) is 
introduced to better account for the influence of pressure drop on the heat 
transfer calculations. Introducing the array of pressure ratios, F, into 
the above equation the result is: 


yt + At 



R 

r t+At 

L P 


(Each element of this array, F, is set equal to 1 during the first pass.) 


Temperature Time Derivative Due to Mixing : 


dT 

dt 


(ki - W.) C T + (C W.T. - C W T ) 
\ o i/ P V p.j i l p Q o oj 

due to = MUZ 

mixing p 


(5C) 
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Using numerical integration let 


T t+4t . T* + 41 At 

at 


r t+At 


/..t+At ..t+At\pt+At T t . /pt+At..t+At T t+At 

. ("o - W i ) C p T * ( C p, W i T i 

= H ^ - 


Substituting 


pt+At..t+At T t+At\ 

C P 0 "o T o ) 


M t+At c t+At 

P 


At 


( U t + At . H t+At) _ M* - M t +4t 


At 


T 


t+At 


M 


t+At 


r 1 - + 



Pi 1 


+At pt+Atyt+A^t+At 

, - C p 0 "o T o 


) 


M t+At^t+At 


At 


Temperature Time Derivative Due to Heat Transfer : 


dT 

dt 


hA / t \ dT hA 
due to ~ MC u w " u * T - T = MC 
heat transfer P w p 


dt 


(5d) 


Assume T^ is constant over the time step for the purpose of inte- 
grating the left side with respect to time. This is a reasonable assumption 
since T w changes much more slowly than T due to the relatively large 
heat capacity of the metal. It was also assumed that h and M were con- 
stant over the time step to allow integration of the right side of the 
equation. 



T) 


t+At 


hA 

W~ 

P 


t+At 

t 


LN 


(T - T) 
v w ' 


t+At 


< T w - T >‘ 


hA 

MET 


At 


hA 

WT 


At 


<T„ - T) 


t+At 


< T w - T > e 
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jt+At _ j 


W 


( T w- T > e 


hA 

t 


At 


T t+ At = + (T w - T l )\l 


- hA 


At 


- e 


This equation says that, as the time step is made larger, the gas temperature 
approaches the wall temperature asymptotically. Thus using large time steps 
cannot cause instabilities because of excessive change in gas temperature. 
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APPENDIX D 


EXPANSION AND COMPRESSION SPACE HEAT TRANSFER COEFFICIENTS 

EXPANSION SPACE 


This analysis assumes perfect insulation between the combustion gas and 
the expansion space wall. Heat transfer between the expansion space wall 
and the working space gas is a combination of radiation and convection. For 
radiation: 



and 



(Dl) 


where 

a Boltzmann constant 

F emissivity times view factor 

T w wall temperature 

T gas temperature 

Q rate of heat flow 

A heat transfer area 

h ra d radiation heat-transfer coefficient 

The overall F is assumed to be 0.7 

The convection heat-transfer coefficient is: 


conv 


0.023(Re) 0,8 (Pr) 0 * 4 ~ 

n 


Re > 10 000 


or 


(D2a) 


conv 


0.023(Re)°’ 8 (Pr)°' 4 jj- 


,0.07 


1 + 


2100 < Re < 10 000 


where L is the maximum distance from the cylinder head to the displacer, 
and 


h CO nv = 1-86(GRAETZ) 0,333 GRAETZ > 10; Re < 2100 (D2b) 

h 

or 
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h conv " 5,0 Btu/hr “ ft2 ~ R ^ 28 - 4 m 2 * “ -j GRAETZ < 10; Re < 2100 (D2c) 

where Graetz number = Re X Pr X DH/L. The value in equation (D2c) is an 
assumed cutoff point (close to the natural convection coefficient). For the 
combined heat-transfer coefficient the values obtained from equations (Dl) 
and (D2) are added 


Compression Space 

Since the radiation effect is small in the compression space, only con- 
vection heat transfer is considered. Equation (D2) is used for the calcula- 
tion. It is assumed that the wall temperature is known. Without detailed 
analysis or test data to identify this wall temperature, it seems reasonable 
to assume that it is about equal to the average compression space gas tem- 
perature over the cycle. The net result is that very little heat transfer 
takes place in the compression space and the compression-space process is 
essentially adiabatic. 
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APPENDIX E 

MOMENTUM EQUATION AND DECOUPLED PRESSURE DROP CALCULATIONS 

MOMENTUM EQUATION 


A general form of the conservation of momentum equation for one- 
dimensional flow is: 






= 0 


Rate of 
accumulation 
of momentum 
per unit 
volume 


Rate of 
momentum 
gain by 
convection 
per unit 
volume 


Rate of 

momentum 

gain by 

viscous 

transport 

per unit 

volume 


Rate of 
momentum 
gain by 
pressure 
force per 
unit volume 


(El) 


Expanding the first and second terms of equation (El) yields: 


p ix + y jp + y ? .(.pyJ> + pV ix + JL 

p at at ax p ax 2D. 

h 


f ..2 + aP 


pv 


ax 


= 0 


(E2) 


By the continuity equation: 




Therefore the second the third terms of equation (E2) can be eliminated to 
yield: 


p 



„ 3V + f 
pV Tx ID^ 


pv 



(E3) 


The first term in equation (E3) is neglected in the model. Neglecting this 
term and multiplying the resulting equation by ax/p yields: 


vav 



2 


v ax 



(E4) 


Note that at zero flow the second and third terms of equation (E3) are zero, 
so that it reduces to 


45 



3 V 

7C 


If. = 0 

ax 


in which case the time derivative term is responsible for any pressure 
drop. The significance of this time derivative term could be investigated 
by the use of a model which uses the complete momentum equation such as that 
of Uriel i (ref. 18) or that of Schock (ref. 19). 


Decoupled Pressure Drop Calculations: 

The pressure drop calculations are decoupled from the basic thermody- 
namic calculations for the working space; this decoupling of pressure drop 
allows use of a larger time step (and less computing time) than would other- 
wise be possible with the explicit, one iteration per time step, numerical 
integration used in the model. 

The effect of pressure drop on the thermodynamic calculations is ac- 
counted for as follows: 

(1) Engine work, heat in and heat out per cycle are calculated assuming 
no pressure drop. Pressure variation with time over the cycle is the same 
at all control volumes in the working space. Gas flow rates are, therefore, 
dependent only on the variable volumes and the fluctuations in the working 
space gas temperatures. 

(2) Using the "no pressure drop" gas flow rates calculated in step (1) 
above, pressure drops are calculated across each control volume; the 
pressure variation with time at the center of the regenerator is the same as 
in step (1). 

(3) The pressure drops calculated in step (2) are summed up from the 
expansion space to the center of the regenerator and from the center of the 
regenerator to the compression space at each time step. Pressure variations 
for expansion and compression space are corrected for pressure drop. 

(4) The Ap corrected pressure variations in the expansion and 
compression spaces are used to recalculate expansion space work, com- 
pression space work and total engine work per cycle. The difference between 
the works calculated in (1), assuming no Ap, and the Ap corrected works 
yield: 


(a) total work loss per cycle due to Ap 


(b) work loss from expansion space to the center of the regenerator 
(hot side of working space) due to Ap. 

(c) work loss from center of regenerator to compression (cold side 
of working space) due to Ap. 
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(5) The heat into the engine is now corrected for Ap by assuming: 


Heat 

into 

engine with 
aP 


Heat 

into 

engine without 
aP 


Work loss 
due to aP 
on hot side of 
working space 


The heat out of the engine is corrected by assuming: 


Heat 
out of 
engine with 
AP 


Heat 
out of 
engine 


Work loss 
+ due to aP 

without on cold side of 
aP working space 


All calculations up to this point are completed during the first 
NOCYC = 25 cycle pass through the model calculations. The Ap corrections 
discussed above were the only ones made in the model of reference 1. An 
improved correction for the effect of ap has been incorporated into the 
model by adding the following steps to the above: 

(6) During the last cycle of the first pass, store the pressure varia- 
tions over the cycle, corrected for Ap, for each control volume. A conven- 
ient way to do this is to store the ratio of pressure at the control volume 
to pressure at the center of the regenerator, for each control volume at 
each time step over the cycle. Thus an array of pressure ratios is created 
which documents the effect of Ap on the pressure variations at each control 
volume over the cycle. 

(7) A second pass (of NOCYC = 25 more cycles) is made through the model 
calculations. This time, instead of assuming a uniform pressure variation 
with time throughout the working space in making the thermodynamic calcula- 
tions, the array of pressure ratios is used to infer pressure variation, 
corrected for Ap, at each control volume. A calculation of engine work with 
no Ap (that is doing the expansion and compression space work integrations 
using the pressure at the center of the regenerator) is still made to allow 
a calculation of work loss due to Ap. Now, however, the correction to the 
heat transfer in and out over the cycle (as in (5)) is no longer necessary; 
this is because the effect of pressure drop (via the pressure variations at 
each control volume over the cycle) is now an integral part of the heat and 
mass transfer calculations at each time step. 

The second pass through the calculations does not significantly change 
the calculated work loss due to Ap. However, the heat into the heater (and 
the expansion space work) increases more than the heat out of the cooler 
(and the compression space work). Thus there is a net increase in the basic 
work and power (that is, work and power before Ap loss) calculated for the 
engine. For hydrogen at design P-40 conditions (15 MPa, 4000 RPM) the ef- 
fect of the second pass was to increase brake power by 1.2 kW (3.5%). For 
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helium at design P-40 conditions the increase is about 1.9 kW 
(9.7 percent). Since the model usually overpredicts power, the additional 
correction increases the error in predicted power at the design point; it 
did, however, cause the shape of the predicted curve to approximate more 
closely that of the experimental curve. 

In developing the second pass correction, it was found that recalcula- 
ting the pressure ratio array, F, at the end of the second and then making a 
third pass had negligible effect on the predicted performance. No attempt 
was made to optimize the correction procedure to get minimum computing 
time. For example, it may be possible to use fewer cycles during the first 
pass and get the same accuracy. 
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APPENDIX F: SYMBOLS USED IN FORTRAN SOURCE PROGRAMS, INPUT 
DATA SETS, AND OUTPUT DATA SETS 


A Ratio of inlet and outlet areas for flow coefficient calculation 

AC Compression space work per increment of crank angle, ft-lbf/rad 

( J/rad) 

ACAN Heat conduction area for external insulation container, in 
(CM 2 ) 

ACDUC Compression space work, with only cooler pressure drop, per 
increment of crank angle, ft-lbf/rad (J/rad) 

ACDUE Compression space work, with only end-effects pressure drop, per 
increment of crank angle, ft-lbf/rad (J/rad) 

ACDUR Compression space work, with only regenerator pressure drop, per 
increment of crank angle, ft-lbf/rad (J/rad) 

ACONDD Heat conduction area through piston, in 2 (cm 2 ) 

AC02 Coefficient in quadratic equation for specific heat of carbon 
dioxide, Btu/lbm-°R (J/kg-K) 

ACP Compression space work with no pressure drop per increment of crank 
angle, ft-lbf/rad (J/rad) 

ACS Array of control volume cross-sectional flow areas, 

2 2 
in /regenerator flow path (cm /regenerator flow path) 

ACSCOM Effective compression space cross-sectional flow area, used for end 

O O 

effects pressure drop calculation, in (cm) 

ACSEXP Effective expansion space cross-sectional flow area, used for end 

p p 

effects pressure drop calculation, in (cm ) 

ACSO Alternate storage array for array ACS 

AE Expansion space work per increment of crank angle, ft-lbf/rad 

( J/rad) 

AEALT Expansion space work per increment of crank angle if all pressure 
drop is calculated relative to the pressure in the compression 
space, ft-lbf/rad (J/rad) 

AEDUE Expansion space work, with only end effects pressure drop 

considered, per increment of crank angle, ft-lbf/rad (J/rad) 

AEDUH Expansion space work, with only heater pressure drop considered, 
per increment of crank angle, ft-lbf/rad (J/rad) 

AEDUR Expansion space work , with only regenerator pressure drop 
considered, per increment of crank angle, ft-lbf/rad (J/rad) 
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AEFH20 Effective cooling water flow area through coolers per cylinder, 
in (cm ) 2 2 

AEP Expansion space work, assuming no pressure drop, per increment of 
crank angle, ft-lbf/rad (J/rad) 

AHE Temperature independent term in equation for specific heat at 
constant volume for helium, in-lbf /Ibrn-'R (J/kg-K) 

O 

AHT Array of control volume heat transfer areas, in /regenerator flow 
path (cm 2 /regenerator flow path) 

AHTCTW Heat transfer area of one cooler tube on the water side, ft u 
(cm 2 ) 

AHTO Alternate storage array for array AHT 

AH2 Temperature independent term in equation for specific heat at 
constant volume for hydrogen, in-lbf/lbm-°R ( J /kg— K ) 

AIN Inlet flow area, ft 2 (cm 2 ) 

ALPHA Crank angle, rad 

AMIN Ratio of inlet and outlet areas 

ANGLE Dummy variable used in work integral function definition, rad 

AOUT Outlet flow area, ft 2 (cm 2 ) 

2 2 

AP Piston cross-sectional area, in^ (cm ) 

APCMAX Crank angle at which maximum compression space pressure occurs, deg 

APCMIN Crank angle at which minimum compression space pressure occurs, deg 

APEMAX Crank angle at which maximum expansion space pressure occurs, deg 

APEMIN Crank angle at which minimum expansion space pressure occurs, deg 

APMAR Piston cross-sectional area minus piston rod cross-sectional area, 
in 2 (cm 2 ) 

2 2 

AR Piston rod cross-sectional area, in (cm ) 

2 

AREA Dummy variable used in work integral function definition, in 

AREAIN inlet flow area, in 2 (cm 2 ) 

AREAOT Outlet flow area, in 2 (cm 2 ) 

AS Array of variables equivalent to works per increment of crank angle in 

COMMON /ASET/, ft-lbf/rad (J/rad) 

ATPC Total work per increment of crank angle when pressure drop is 

calculated relative to compression space pressure, ft-lbf/rad (J/rad) 
AUXEFF Engine efficiency, including auxiliary losses 
AUXFP1 Auxiliary power requirement per cylinder, ft-lbf /cycle (J/cycle) 

AUXFP4 Auxiliary power requirement for engine, ft-lbf/cycle (J/cycle) 


50 



AUXHP4 Auxiliary power requirement for engine, hp (kW) 

AUXKW4 Auxiliary power requirement for engine, kW 

AUXLOS Auxiliary power requirement per cylinder, hp (kW) 

AUXPWR Engine brake power (with auxiliary power requirement subtracted), hp 
(kW) 

p p 

AVGPC Time averaged compression space pressure, Ibf/in (N/cm ) 

p p 

AVGPE Time averaged expansion space pressure, Ibf/in (N/crrr) 

AVGPMP Time averaged pressure at center of regenerator, MPa 

p p 

AVGWSP Time averaged pressure at center of regenerator, lbf/in (N/cm ) 
AVPCMP Average compression space pressure, MPa 
AVPEMP Average expansion space pressure, MPa 

A1 Temperature independent term in equation for specific heat at 

constant volume, in-lbf /lbm-°R (J/kg-K) 

p p 

B Coefficient, real gas equation of state, Ibf/in (N/cm ) 

BASICP Indicated power plus pressure drop loss, per cylinder, hp (kW) 

BC02 Coefficient in quadratic equation for specific heat of carbon 
dioxide, BTU/lbm-°R 2 (J/kg-K 2 ) 

BETA Crank angle +PI/2, rad 

BHE Sensitivity of specific heat at constant volume to temperature for 
helium, in-lbf /lbm-*R 2 (J/kg-K 2 ) 

BH2 Sensitivity of specific heat at constant volume to temperature for 
hydrogen, in-lbf /lbm-°R 2 (J/kg-K 2 ) 

BPFP1 Engine brake power per cylinder, ft-lbf /cycle (J/cycle) 

BPFP4 Engine brake power, ft-lbf/cycle (J/cycle) 

BPHP1 Engine brake power per cylinder (=AUXPWR), hp (kW) 

BPHP4 Engine brake power, hp (kW) 

BPKW1 Engine brake power per cylinder, kW 

BPKW4 Engine brake power, kW 

BRKEFF Engine brake efficiency (not including effect of auxiliaries) 

BRKP Engine brake power per cylinder (not accounting for auxiliaries 
losses), hp (kW) 

B1 Sensitivity of specific heat at constant volume to temperature, 

in-lbf /lbm-°R (J/kg-K) 

CANIR Insulation container inside radius, in (cm) 

CANOR Insulation container outside radius, in (cm) 

CCMPDV Cooler-compression space connecting duct volume, in 2 (cm 2 ) 
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CC02 Coefficient in quadratic equation for specific heat of carbon 
dioxide, Btu/lbm-°R 3 (J/kg-K 3 ) 

CDEDV Cooler dead volume per cylinder, in 3 (cm 3 ) 

CFACTR Function of average working space pressure used in calculating 
auxiliary loss, dimensionless 

CHCF Cooler heat transfer coefficient multiplying factor 

CHCFAC Cooler heat transfer coefficient multiplying factor 

CHE Coefficient in equation for specific heat of helium, Btu/lbm- R 

(J/kg-K 3 ) 

o *3 

CH2 Coefficient in equation for specific heat of hydrogen, Btu/lbm- R 
(J/kg-K 3 ) 

CH2P Monatomic thermal conductivity of hydrogen, Btu/in-sec °R (W/cm-K) 
CH2PP Internal thermal conductivity of hydrogen, Btu/in-sec °R (W/cm-K) 
CLRLOD Cooler tube length to diameter ratio 

CMIXP Monatomic thermal conductivity of mixture of hydrogen and carbon 
dioxide, Btu/in-sec-*R (W/cm-K) 

CMIXPP Internal thermal conductivity of mixture of hydrogen and carbon 
dioxide, Btu/in-sec- *R (W/cm-K) 

CMPSCL Compression space clearance volume, in 3 (cm 3 ) 

CNDH20 Thermal conductivity of water, Btu/ft-sec-°R (W/cm-K) 

CNDSS Cooler tube (stainless steel) thermal conductivity, Btu/ft-sec-°R 
(W/cm-K) 

COEF Pressure drop coefficient, dimensionless 

COEFX Pressure drop coefficient, dimensionless 

COND Array of control volume gas thermal conductivities, Btu/m-sec-°R 
(W/cm-K) 

CONDT Heater tube thermal conductivity, Btu/in-sec- °R (W/cm-K) 

CONDTB Conduction length, top to bottom, of external insulation container 
(if used), in (cm) 

C02P Monatomic thermal conductivity of carbon dioxide, Btu/in-sec-°R 
(W/cm-K) 

C02PP Internal thermal conductivity of carbon dioxide, Btu/i n-sec-°R 
(W/cm-K) 

CP Array of control volume interface specific heats at constant 

pressure, Btu/lbm-°R (J/kg-K) 


52 



CPA Array of control volume specific heats at constant pressure, 

Btu/lbrn-^R (J/kg-K) 

CPAO Alternate storage array for array CPA 

CPCGP Specific heat in cold appendix gap, Btu/lbm-°R (J/kg-K) 

CPCGPI Specific heat of gas crossing interface between cold appendix gap and 
compression space, Btu/lbm-°R (J/kg-K) 

CPCYC Time increments (iterations) per engine cycle 
CPHGP Specific heat in hot appendix gap, Btu/lbm-°R (J/kg-K) 

CPHGPI Specific heat of gas crossing interface between hot appendix gap and 

expansion space, Btu/lbm- # R (J/kg-K) 

CPH20 Cooling water specific heat, Btu/lbm-'R (J/kg-K) 

CPM Regenerator matrix specific heat, Btu/lbm-°R (J/kg-K) 

CR1 Initialization constant 

CR2 Initialization constant 

CR3 Initialization constant 

CTBID Cooler tube inside diameter, in (cm) 

CTBL Cooler tube length, in (cm) 

CTBOD Cooler tube outside diameter, in (cm) 

CTBPCN Cooler tubes per cylinder 

CV Array of control volume specific heats at constant volume, Btu/lbm-°R 

(J/kg-K) 

CVF Function for calculating specific heat at constant volume, Btu/lbm-°R 

(J/kg-K) 

CYLDMB Cylinder distance between middle and bottom wall temperatures, in (cm) 

CYLDTM Cylinder distance between top and middle wall temperatures, in (cm) 

CYLIR Cylinder housing inside radius, in (cm) 

CYLORB Cylinder outside radius at bottom temperature, in (cm) 

CYLORM Cylinder outside radius at middle temperature, in (cm) 

CYLORT Cylinder outside radius at top temperature, in (cm) 

DALOSS Design auxiliary loss — four cylinders, hp (kW) 

DE Hydraulic diameter, ft (cm) 

p p 

DELP Array of pressure drops across control volumes, Ibf/in (N/cm ) 

p p 

DELPCL Pressure drop across cooler, Ibf/in (N/cm ) 

DELPHI Pressure drop across heater, Ibf/in^ (N/cm^) 

p p 

DELPRG Pressure drop across regenerator, Ibf/in (N/cm ) 

DELTIM Engine cycle period, sec 
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DELTM Change in regenerator matrix temperature from one control volume to 
the next, °R (K) 

DFLOSS Design mechanical friction loss, hp (kW) 

DFREQ Design engine frequency, Hz 

DGAPDV Piston-cylinder gap dead volume, in (cm ) 

DH Array of control volume hydraulic diameters, in (cm) 

DHO Alternate storage array for array DH, in (cm) 

DHX Hydraulic diameter, in (cm) 

DISPD Piston diameter, in (cm) 

DISPRD Piston rod diameter, in (cm) 

O O 

DNSTY Array of control volume gas densities, lbm/in (kg/cm ) 

O O 

DP Pressure drop, 1 bf /in (N/cm ) 

DPCLR Cooler pressure drop, 1 bf /i (N/cm 2 ) 

p p 

DPECLD End effects pressure drop, cold side of engine, 1 bf /i n (N/cm ) 

p p 

DPEHOT End effects pressure drop, hot side of engine, lbf/in (N/cm ) 

p p 

DPFRIC Total pressure drop excluding end effects, lbf/in (N/cm ) 

DPHTR Heater pressure drop, lbf/in 2 (N/cm 2 ) 

p p 

DPRCLD Regenerator pressure drop, cold side, lbf/in (N/cm ) 

p p 

DPRHOT Regenerator pressure drop, hot side, lbf/in (N/cm ) 

DPSI Crank angle increment, rad 

DPSUM Total pressure drop, lbf/in 2 (N/cm 2 ) 

DPX Array of pressure drops, lbf/in 2 (N/cm 2 ) 

DRPM Design engine speed, rpm (Hz) 

DSPGAP Gap width between piston and displacer, in (cm) 

DSPHGT Piston height (used in piston-cylinder gap dead volume calculation), 
in (cm) 

DSPWTH Piston wall thickness, in (cm) 

DTCGP Change in cold appendix temperature that would occur due to change in 
pressure ( if appendix gap were not isothermal), °R (K) 

DTCYL Cylinder housing temperature difference between thermocouple 
locations, when calculated by code, °R (K) 

DTGA Change in control volume gas temperature due to heat transfer between 
gas and metal, °R (K) 

DTGASL One-half of the assumed change in gas temperature across the 
regenerator control volume, °R (K) 
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DTHGP Change in hot appendix gap gas temperature that would occur due to 
change in pressure level, °R (K) 

DTIME Time increment, sec 

DTM Array of regenerator matrix temperature corrections, *R (K) 

DTREG Regenerator housing temperature difference between thermocouple 

locations, when calculated by code, °R (K) 

DWCMP Increment in compression space work for one crank angle increment, 
ft-lbf (J) 

DWEXP Increment in expansion space work for one crank angle increment, 
ft-lbf (J) 

E Crank eccentricity (was used in rhombic drive simulation) 

ECTBL Effective cooler tube length for heat transfer, in (cm) 

EFFTOT Engine indicated efficiency 

EHTBL Effective heater tube length for heat transfer, in (cm) 

EID Engine identification (alphanumeric) 

ENFCTR Enthalpy flow from cooler to regenerator per cycle, Btu (J) 

ENFHTR Enthalpy flow from heater to regenerator per cycle, Btu (J) 

ENFRTC Enthalpy flow from regenerator to cooler per cycle, Btu (J) 

ENFRTH Enthalpy flow from regenerator to heater per cycle, Btu (J) 

O O 

EN3PM Rate at which working space volume is swept out, in J /min (cm /min) 
ETYPE Engine type 

q q 

EXPHDV Expansion space-heater connecting duct volume, in (cm ) 

q q 

EXPSCL Expansion space clearance volume, in cm ) 

F Array of flow rates at control volume interfaces, lbm/sec (kg/sec) 

FACT1 Coefficient used in regenerator matrix temperature convergence method 

FACT2 Coefficient used in regenerator matrix temperature convergence method 

FAVG Array of average flow rates for each control volume, lbm/sec (kg/sec) 

FAVG2 Average control volume gas flow rate, lbm/sec (kg/sec) 

FCGP Flow rate between compression space and cold appendix gap, lbm/sec 
(kg/sec) 

FCND1,FCND11,FCND12 Functions of mass fractions and thermal conductivity of 
pure gases used in calculating thermal conductivity of mixture of 
gases 

FCND2,FCND21,FCND22 Functions of mass fractions of pure gases used in 
calculating thermal conductivity of mixture of gases 
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FCNPP1 ,FCNPP2 Functions of mass fractions and thermal conductivity of pure 
gases used in calculating thermal conductivity of mixture of gases 
FCTR Dimensionless function of heat transfer between gas and metal 
FDEN Parameter used in calculating matrix of coefficients for regenerator 
matrix temperature convergence method 
FHGP Flow rate between expansion space and hot appendix gap, Ibm/sec 
(kg/sec) 

FICLR Gas flow rate at hot end of cooler per cylinder, Ibm/sec (kg/sec) 

FIHTR Gas flow rate at hot end of heater per cylinder, Ibm/sec (kg/sec) 

FIJ,FIJ1,FIJ2 Functions of mass fractions and viscosities of pure gases for 
calculating viscosity of a mixture of gases 
F IK Array of — control volume pressure/pressure at center of 
regenerator — (values for each time increment over cycle) 

FIKS Alternate storage array for array F IK 

FIPCV On-off switch used to modify calculation of heat out and heat into 
engine 

FIREG Gas flow rate at hot end of regenerator per cylinder, Ibm/sec (kg/sec) 
FLFP1 Friction loss per cylinder, ft-lbf/cycle (J/cycle) 

FLFP4 Friction loss for engine, ft-lbf/cycle (J/cycle) 

FLHP4 Friction loss for engine, hp (kW) 

FLKW4 Friction loss for engine, kW 

FL0H20 Cooling water flow flow rate per cylinder, Ibm/sec (kg/sec) 

p 

FLOPUA Cooling water flow rate per unit cross-sectional area, Ibm/sec-in 

p 

(kg/sec-cm ) 

FLOW Absolute value of gas flow rate per regenerator flow path, Ibm/sec 
(kg/sec) 

FLOWIN Gas flow rate per regenerator flow path, Ibm/sec (kg/sec) 

FMA21 Function of inlet and outlet Mach numbers, dimensionless 
FMULT Multiplier for overall pressure drop 
FMULTR Multiplier for regenerator pressure drop 

FNUM Parameter used in calculating matrix of coefficients for regenerator 
matrix temperature convergence method 

FOA Estimate of effectiveness of radiation heat transfer from metal to 
gas in expansion space (emissivity * view factor) 

FOCLR Gas flow rate at cold end of cooler per cylinder, Ibm/sec (kg/sec) 

FOEXP Gas flow rate at exit of expansion space, Ibm/sec (kg/sec) 
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FOHTR Gas flow rate at cold end of heater per cylinder, Ibm/sec (kg/sec) 

FOREG Gas flow rate at cold end of regenerator per cylinder, Ibm/sec 

(kg/sec) 

FREQ Engine speed, Hz 

FRIN Gas flow rate at hot end of regenerator per cylinder, Ibm/sec (kg/sec) 
FRLOSS Friction loss per cylinder, hp (kW) 

FO Variable used in definition of function for Simpson rule integration 

(value of integrand at two time increments before current time) 

FI Variable used in definition of function for Simpson rule integration 

(value of integrand at one time increment before current time) 

F2 Variable used in definition of function for Simpson rule integration 

(value of integrand at current time) 

GAMMA Ratio of gas specific heats (CP/CV) 

GAMMA1 =GAMMA for adiabatic flow, =1.0 for isothermal flow 
GPMH20 Cooling water flow rate per cylinder, gal/min (gm/sec) 

GRAETZ Dimensionless number for calculating convection heat transfer in 
expansion space 

GRAV Constant, 32.2 lbm-ft/lbf-sec 2 (1.0 kg-M/N-sec 2 ) 

? © 

H Array of gas to metal heat transfer coefficients, Btu/sec-in - R in 

subroutine HEATX (W/cm -K) units converted to Btu/sec-ft - R in 
subroutine ROMBC 

HA Array of — heat transfer coefficients * heat transfer area — between 

gas and wall), Btu/sec-°R (W/K) 

HACYC Array of average heat transfer coefficients over the engine cycle, 
Btu/sec-ft 2 -°R (W/cm 2 -K) 

HAWC Dimensionless ratio — (heat transfer between gas and wall per deg of 
temperature difference) /(control volume heat capacity) 

HCONV Convection heat transfer coefficient in expansion space, 

Btu/sec-in 2 -°R (W/cm 2 -°K) 

HDEDV Heater dead volume per cylinder, in^ (cm^) 

HFACT Dimensionless factor used in calculating heat transfer coefficients 

HHCF Dimensionless factor used in calculating heater heat transfer 

coefficients 

HHCFAC Heater heat transfer coefficient multiplying factor 
HLOD Array of heat transfer coefficient function values for different tube 
length/diameter ratios, dimensionless 
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HMX Array of maximum values of heat transfer coefficients over the cycle, 
Btu/sec-ft 2 -°R (W/cm 2 -K) 

HRAD Effective heat transfer coefficient for radiation heat transfer in 

o 0 p 

expansion space, Btu/sec-in - R (W/cm-K) 

O O 

HRDV Heater-regenerator connecting duct volume, in (cm°) 

HTABL Table of values of heat transfer correlation, dimensionless 
HTBID Heater tube inside diameter, in (cm) 

HTBL Heater tube length, in (cm) 

HTBOD Heater tube outside diameter, in (cm) 

HTBPCN Number of heater tubes per cylinder 
HTRLOD Heater tube length over diameter ratio 

HWATR1 Clean tube, cooler tube to water heat transfer coefficient, 

Btu/sec-ft2-°R (W/cm 2 -K) 

HWATR2 Fouled tube, cooler tube to water heat transfer coefficient, 
Btu/sec-ft 2 -°R (W/cm 2 -K) 

I Index 

ICOND Index; =1 calculate cylinder and regenerator housing temperatures 
from TM( 1) , TM(4), and TH20IN, =0 use input values for housing 
temperatures 
IDEX Index 

IDRUN Alphanumeric identifier for input operating conditions 
IJK Index 

IK Index 

IMIX Index, =1 to calculate performance for mixture of hydrogen and carbon 
dioxide 

=0 to calculate performance for pure hydrogen or pure helium 
I0UT Index used as on-off switch for portion of output (that which goes 
into Tables VII and VIII) 1-on, 0-off 
IP Index used to control printout 

IPCV Index: =0 makes first pass through calculations using uniform 

pressure in calculating flow rates. Then, make second pass through 
calculations using pressure array, FIK, (created in first pass) in 
calculating flow rates. = 1 means eliminate second pass. 

IPLOT Counter used in storing data for plotting 
IPRINT Index used to control printout 

IPRNTO Index used to control printout 
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IPRNT2 Index used to control printout 

I PUMP Index: =1 means pumping loss due to piston cylinder gap is included 

=0 means pumping loss not included 

IP1 Index 

IRE Index 

IREV Index 

ISCD Index; =1 for separate connecting duct volumes 

=0 to lump connecting duct volumes with adjacent control 

volumes 

ISIMP Index 

I START Index 

ITER Counts total numer of iterations (time increments) since beginning of 
run 

ITMPS Index: =1 to print temperature arrays at each time increment (for 
check out) 

=0 don't print temperature arrays at each time increment 
(normal ly=0) 

ITPCYC Number of iterations (time increments) per cycle 

ITR Counts iterations (time increments) since beginning of cycle 

ITRM1 ITR-1 

I VAR Number of iterations in 5 sec 
J Index 

JCYCLE Index, counts number of cycles 
JI Index 

JIP Index: >0 for short form printout in stored dataset 

=0 long form printout 

JIP1 Index 

JJ Index 

JM Index 

JN Index 

K Index 

KK Index 

KI Index 

KIDEX Index used in updating cooler tube temperature 
KJK Index 

KTRIG Index used in updating cooler tube temperature 
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KTYPE Index, specifies type of pressure drop calculation to be made 
KWRITE Index 
L Index 

M Index 

MAPLOT Index: =1 means store data for plotting 

=0 means don't store data for plotting 
MWGAS Index: =2 for hydrogen working gas 

=4 use helium working gas 

N Index 

NA Index 

NC Number of cooler control volumes 

NCL Index number of last (nearest the compression space) cooler control 
volume 

NCLP1 NCL+1 

NCOND Index used to prevent the conduction subroutine from being called 
more than once per cycle 

NCS Index number of compression space control volume 

NCV Total number of control volumes 

NCVM1 NCV-1 

NCVP1 NCV+1 

NCVP2 NCV+2 

NCVP3 NCV+3 

NCVP4 NCV+4 

NCI Index number of first (nearest regenerator) cooler control volume 
NC1M1 NCI-1 

NC1P1 NC1+1 

NES Index number of expansion space control volume 

NH Number of heater control volumes 

NHC Index number of center heater control volume if there are an odd 
number of heater control volumes (=NH1 + (NH-l)/2) 

NHL Index number of last heater control volume 

NHLP1 NHL+1 

NH1 Index number of first (nearest expansion space) heater control volume 
NH1M1 NH1-1 

NH1P1 NH1+1 

NITPC Number of time increments per cycle 
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NOCYC Number of engine cycles 

NOEND Number of cycle at which regenerator matrix temperature convergence 
procedure is to be turned off 

NPASS Index automatically set by program on basis of input value of index, 
IPCV. 

If IPCV=0, then NPASS is set =2 to get two passes through 
calculations. 

If IPCV=1, then NPASS is set =1 to get one pass only through 
calculations. 

NPLOTS Number of variables to be plotted 

NR Number of regenerator control volumes 

NRC Index number of center regenerator control volume 

NRCM1 NRC-1 

NRCP1 NRC+1 

NRL Index number of last regenerator control volume 
NRLM2 NRL-2 

NRLP1 NRL+1 

NR1 Index number of first (nearest heater) regenerator control volume 
NR1M1 NR 1-1 

NR1M2 NR1-2 

NR1M3 NR1-3 

NSTRT Number of cycle at which regenerator matrix temperature convergence 
procedure is turned on 

OLDTIM Time at end of cycle, to use in calculating period, sec 
OMEGA Engine frequency, Hz 

p 

P Pressure at center of regenerator, lbf/in (MPa) 

p 

PC Pressure in compression space, lbf/in (MPa) 

PCDUC Pressure in compression space when only cooler pressure drop is 

p 

accounted for, lbf/in (MPa) 

PCDUE Pressure in compression space when only end effects pressure drop is 

p 

accounted for, lbf/in (MPa) 

PCDUR Pressure in compression space when only regenerator pressure drop is 

p 

accounted for, lbf/in (MPa) 

p 

PCMAX Maximum compression space pressure, lbf/in (MPa) 

p 

PCMIN Minimum compression space pressure, lbf/in (MPa) 

PCMNP Minimum compression space pressure, Mpa 
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2 

PCMP Alternate storage location for compression space pressure, lbf/in 
(MPa) 

PCMXP Maximum compression space pressure, Mpa 

PCSUM Summation of compression space pressures used in calculating time 

p 

averaged compression space pressure, lbf/in (MPa) 

p 

PCV Array containing control volume pressures, lbf/in (MPa) 

PD Desired mean pressure level, lbf/in 2 (MPa) 

PDFP4 Engine pressure drop loss, ft-lbf/cycle (J/cycle) 

PDHP4 Engine pressure drop loss, hp (kW) 

PDKW4 Engine pressure drop loss, kW 

PE Expansion space pressure, lbf/in 2 (MPa) 

PEALT Expansion space pressure when pressure drop is calculated relative to 
pressure in compression space (instead of pressure at center of 

O 

regenerator), lbf/in (MPa) 

PEDUE Expansion space pressure when only end effects pressure drop is 

p 

considered, lbf/in c (MPa) 

PEDUH Expansion space pressure when only heater pressure drop is 
considered, lbf/in 2 (MPa) 

PEDUR Expansion space pressure when only regenerator pressure drop is 
considered, lbf/in 2 (MPa) 

p 

PEMAX Maximum expansion space pressure, lbf/in (MPa) 

p 

PEMIN Minimum expansion space pressure, lbf/in (MPa) 

PEMNMP Minimum expansion space pressure, MPa 

PEMXMP Maximum expansion space pressure, MPa 

PERREB Percent error, engine energy balance 

PESUM Summation of expansion space pressures used in calculating time 

p 

averaged expansion space pressure, lbf/in (MPa) 

p 

PEXP Alternate storage location for expansion space pressure, lbf/in^ 

(MPa) 

PHASE Angle by which the compression volume lags the expansion volume, deg 
PI Constant=3. 14159265 

PIN Pressure at control volume inlet, lbf/in 2 (MPa) 

PI02 PI/2 

PI04 PI/4 

PLOT Array in which variables to be plotted are stored 

PMEAN Alternate storage location for mean pressure, lbf/in 2 (MPa) 
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POLD Value of reference pressure at time increment previous to current 
time, lbf/in^ (MPa) 

POUT Pressure out of control volume, lbf/in (MPa) 

PR Prandtl number, dimensionless in subroutine HEATX — or — 

pressure ratio (POUT/PIN) in subroutine XDEL 
PRATAV Pressure ratio— (PEMAX+PCMAX) /(PEMIN+PCMIN) 

PRATC Pressure ratio — PCMAX/PCMIN 

PRATE Pressure ratio — PEMAX/PEMIN 

PREGER Percent error in regenerator energy balance 

PRH20 Prandtl number for cooling water flow 

PRIN Pressure at hot end of regenerator, lbf/in (MPa) 

PROSTY Regenerator matrix porosity 

p 

PROUT Pressure at cold end of regenerator, lbf/in (MPa) 

PRSUM Summation of pressures at center of regenerator used to calculate 

p 

time averaged working space pressures, lbf/in (MPa) 

PS Array of variables equivalent to pressures in COMMON /PSET / 

PSI Crank angle, radians in subroutine ROMBC — or — conversion constant, 
1/144=0.006945 ft^/in^ in subroutine XDEL 
PSIDEG Crank angle, deg 

PSIOLD Value of crank angle at time increment before current time, rad 

PWRFP1 Engine indicated power per cylinder, ft-lbf /cycle (J/cycle) 

PWRFP4 Engine indicated power, ft-lbf /cycle (J/cycle) 

PWRHP Engine indicated power per cylinder, hp (kW) 

PWRHP4 Engine indicated power, hp (kW) 

PWRKW1 Engine indicated power per cylinder, kW 

PWRKW4 Engine indicated power, kW 

PXIN Pressure, lbf/in^ (MPa) 

Q Array of control volume heat transfers from gas to metal, Btu/sec (W) 

QADD Control volume heat transfer between gas and metal for one time 

increment, ft-lbf (J) 

QAPGAP Appendix gap loss per cylinder, ft-lbf/cycle (J/cycle) 

QBTUPS Rate of heat out through cooling water per cylinder, Btu/sec (W) 

QCAN Rate of heat conduction through external insulation container, 

Btu/sec (not used in P40-model) 

QCGP Rate of heat transfer between wall and cold appendix gap, Btu/sec (W) 

QCGPS Appendix gap loss, cold end of piston, ft-lbf/cycle (J/cycle) 
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QCLEXF Heat out through cooling water per cylinder, ft-lbf /cycle (d/cycle) 

QCLOUT Heat out through cooling water, excluding heat generated by 

mechanical friction, ft-lbf/cycle (J/cycle) 

QCLRSV Alternate storage location for net heat out through cooler per cycle, 
ft-lbf/cycle (J/cycle) 

QCNDCL Heat conducted through cyclinder housing, ft-lbf/cycle (J/cycle) 
QCNDCN Heat conducted through insulation container, ft-lbf/cycle (J/cycle) 
QCNDD Heat conducted through piston walls, ft-lbf/cycle (J/cycle) 

QCNDRI Heat conducted into hot end of regenerator housing, ft-lbf/cycle 
(J/cycle) 

QCNDRO Heat conducted out of cold end of regenerator housing, ft-lbf/cycle 
(J/cycle) 

QCNDTI Heat into engine via conduction (includes shuttle) per cylinder, 
ft-lbf/cycle (J/cycle) 

QCNDTO Heat out of engine via conduction (includes shuttle) per cylinder, 
ft-lbf/cycle (J/cycle) 

QCOLN Cooler heat transfer from metal to gas for one time increment, ft-lbf 

(J) 

QCOLP Cooler heat transfer from gas to metal for one time increment, ft-lbf 

(J) 

QCOM Compression space heat transfer for one time increment, ft-lbf (J) 
QCOMP Net heat transferred from gas to metal in compression space, 
ft-lbf/cycle (J/cycle) 

QCOMPN Heat transferred from metal to gas in compression space, ft-lbf/cycle 
(J/cycle) 

QCOMPP Heat transferred from gas to metal in compression space, ft-lbf/cycle 
(J/cycle) 

QCONDD Rate of heat conduction through piston wall, Btu/sec (W) 

QCOOL Net cooler heat transfer for one time increment, ft-lbf (J) 

QCOOLN Heat transferred from gas to metal in cooler, ft-lbf/cycle (J/cycle) 

QCOOLP Heat transferred from metal to gas in cooler, ft-lbf/cycle (J/cycle) 

QCOOLR Net heat transferred from gas to metal in cooler, ft-lbf/cycle 
(J/cycle) 

QCREG Rate of heat transfer through regenerator housing, Btu/sec (W) 

QCRIN Heat conduction rate into hot end of regenerator housing, btu/sec (W) 
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QCROUT Heat conduction rate out of cold end of regenerator housing, Btu/sec 

(W) 

QCYL Rate of heat transfer through cylinder housing, Btu/sec (W) 

QCYL1 Heat conduction rate from hot end to middle of cylinder housing, 

Btu/sec (W) 

QCYL2 Heat conduction rate from middle to cold end of regenerator housing, 
Btu/sec (W) 

QEIN Net heat rate to engine per cylinder, ft-lbf/cycle (J/cycle) 

QEOUT Net heat from engine per cylinder (includes heat out via cooling 

water plus auxiliary losses), ft-lbf/cycle (J/cycle) 

QEX Expansion space heat transfer for one time increment, ft-lbf (0) 

QEXP Net heat transferred from gas to metal in expansion space, 

ft-lbf/cycle (J/cycle) 

QEXPN Heat transferred from metal to gas in expansion space, ft-lbf/cycle 
(J/cycle) 

QEXPP Heat transferred from gas to metal in expansion space, ft-lbf/cycle 
J/cycle) 

QHEAT Heater heat transfer for one time increment, ft-lbf (J) 

QHEATN Heat transferred from metal to gas in heater, ft-lbf/cycle (J/cycle) 

QHEATP Heat transferred from gas to metal in heater, ft-lbf/cycle (J/cycle) 

QHEATR Net heat transferred from gas to metal in heater, ft-lbf/cycle 
(J/cycle) 

QHETN Heater heat transfer from metal to gas for one time increment, ft-lbf 

(J) 

QHETP Heater heat transfer from gas to metal for one time increment, ft-lbf 

(J) 

QHGP Rate of heat transfer between cylinder wall and hot appendix gas, 
Btu/sec (W) 

QHGPS Appendix gap loss, hot end of piston, ft-lbf/cycle (J/cycle) 

QIN Heat into engine per cylinder (accounts for heating effect of 

pressure drop loss in hot end of engine), ft-lbf/cycle (J/cycle) 

QINB Heat into engine via heater and expansion space per cylinder, 
ft-lbf/cycle (J/cycle) 

QINFP4 Heat into engine (accounts for heating effect of pressure drop loss 
in hot end of engine), ft-lbf/cycle (J/cycle) 
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QINHP4 Heat into engine (accounts for heating effect of pressure drop loss 
in hot end of engine), hp (kW) 

QINKW1 Heat into engine per cylinder (accounts for heating effect of 
pressure drop loss in hot end of engine), kW 
QINKW4 Heat into engine (accounts for heating effect of pressure drop loss 
in hot end of engine), kW 

QOA Array of control volume heat transfer rates per unit area, 
Btu/sec-in 2 or Btu/sec-ft 2 (W/cm 2 ) 

QOAMX Array of control volume maximum heat transfer rates per unit area, 
Btu/sec-ft 2 (W/cm 2 ) 

QOAAVG Array of control volume average heat transfer rates per unit area, 
Btu/sec-ft 2 

QOAMX Array of control volume maximum heat transfer rates per unit area, 
Btu/sec-ft 2 (W/cm 2 ) 

Q0TFP4 Heat out through cooling water for engine, ft-lbf /cycle (J/cycle) 

Q0THP1 Heat out through cooling water per cylinder, hp (kW) 

Q0THP4 Heat out through cooling water for engine, hp (kW) 

Q0TKW1 Heat out through cooling water per cylinder, kW 

QINKW4 Heat out through cooling water for engine, kW 

QOUT Net heat flow to coolant (larger than heat flow to cooler by 

mechanical losses), ft-lbf/cycle (J/cycle) 

QOUTB Heat out through cooler and compression space per cylinder, 
ft-lbf/cycle (J/cycle) 

p 

QRAD Rate of radiation heat transfer in expansion space, Btu/ft -hr 
(W/cm 2 ) 

QREGE Regenerator heat transfer for one time increment, ft-lbf (J) 

QREGEN Net heat flow from gas to metal in regenerator, ft-lbf/cycle 
(J/cycle) (should be close to zero for convergent solution) 

QREGN Heat flow from metal to gas in regenerator, ft-lbf/cycle (J/cycle) 

QREGP Heat flow from gas to metal in regenerator, ft-lbf/cycle (J/cycle) 

QRGN Regenerator heat transfer from metal to gas for one time increment, 

ft-lbf (J) 

QRGP Regenerator heat transfer from gas to metal for one time increment, 
ft-lbf (J) 

QSHTL Piston shuttle loss, ft-lbf/cycle (J/cycle) 

QSHTTL Rate of heat loss via piston shuttle, Btu/sec (W) 
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R Gas constant, in-lbf /lbm-*R (J/kg-K) 

RC Array of dimensionless coefficients used in regenerator matrix 

temperature convergence method 

3 

RCDV Volume of regenerator-cooler connecting ducts per cylinder, in 
(cm 3 ) 

RCRANK Crank radius, in (cm) 

3 3 

RDEDV Regenerator dead volume per cylinder, in (cm ) 

RE Reynolds number 

REALGS =1 for real gas equation of state, =0 for ideal gas equation of state 

RECD1 Reynolds number in expansion space-heater connecting duct (based on 

average of inlet and outlet flow rates) 

RECD2 Reynolds number in heater-regenerator connecting duct 

RECD3 Reynolds number in regenerator-cooler connecting duct 

RECD4 Reynolds number in cooler-compression space connecting duct 
RECMP Reynolds number at entrance of compression space 

REEXP Reynolds number at exit of expansion space 

REFF1 Measure of regenerator effectiveness (ENFRTH/ENFHTR) 

REFF2 Measure of regenerator effectiveness (ENFCTR/ENFRTC) 

REGDMB Regenerator housing distance between middle and bottom temperature 

measurement locations, in (cm) 

REGDTM Regenerator housing distance between top and middle temperature 
measurement locations, in (cm) 

REGID Regenerator inside diameter (matrix diameter), in (cm) 

REGIR Regenerator housing inside radius, in (cm) 

REGL Regenerator matrix length, in (cm) 

REGORB Regenerator housing outside radius, bottom, in (cm) 

REGORM Regenerator housing outside radius, middle, in (cm) 

REGORT Regenerator housing outside radius, top, in (cm) 

REGPCN Number of regenerators per cylinder 
REH20 Reynolds number for cooling water flow rate 

REIC Reynolds number in cooler control volume nearest the regenerator 

REIH Reynolds number in heater control volume nearest the expansion space 

REIR Reynolds number in regenerator control volume nearest the heater 
REOC Reynolds number in cooler control volume nearest the compression space 

REOH Reynolds number in heater control volume nearest the regenerator 

REOR Reynolds number in regenerator control volume nearest the cooler 
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RETABL Array of Reynolds number values corresponding to values of the heat 
transfer correlation in array HTABL 
REYNO Array of working gas control volume Reynolds numbers 

p o 

R6AREA Regenerator cross-sectional flow area with no matrix, in (crrr) 

RHCFAC Regenerator heat transfer coefficient multiplier 
RHCPCV Regenerator matrix heat capacity per control volume (per cylinder) 
Btu/°R (J/K) 

RH0H20 Density of water, lbm/ft^ (gm/cm^) 

RH20 Resistance to heat flow from cooler tube to cooling water, sec-°R/Btu 
(K/W) 

o o 

RMDEN Regenerator matrix metal density, Ibm/in (gm/cm ) 

O Q 

RO Working gas density, lbm/ft (gm/cm ) 

RODL Connecting rod length, in (cm) 

ROX Working gas density, lbm/in (gm/cm ) 

RP Gas constant, btu/lbm- 0 R (J/kg-K) 

RPM Engine speed, rpm (Hz) 

RPMF Speed for auxiliary requirement calculation (assumes design speed of 
4000 rpm), rpm (Hz) 

RPMFN Speed for auxiliary loss requirements correction, thousands of rpm 
RTUBE Cooler tube wall resistance, sec-°R/Btu (*K/W) 

RWIRED Regenerator matrix wire diameter, in (cm) 

SAVET Save time required for first nocyc cycles before resetting time=0 

SET Array of variables equivalent to work variables in COMMON /RESET/ 

SRULE Function used in performing Simpson rule integration 
STORE Array used to store cycle quantities to allow calculation of averages 
over five cylcles, various dimensions 
STROKE Piston stroke, in (cm) 

SUM Summation over gas control volumes of (pressure * volume)/ 

(gas constant * temperature), dimensionless 
SUMDEN Array used in calculating time weighted average of difference between 
regenerator matrix and gas temperatures, lbm (kg) 

SUMNUM Array used in calculating time weighted average of difference between 
regenerator matrix and gas temperatures, lbm-°R (kg-K) 

SUMV Summation of (control volumes * array of pressure ratio factors), 
in^ (cm^) 
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SUMWF Summation of (control volume gas inventories * pressure ratio 
factors), Ibm (kg) 

SUMWT Summation of (control volume gas inventories * control volume gas 
temperatures), Ibm-deg R (Kg-°K) 

T1 Expansion space gas temperature, °R (K) 

T2 Compression space gas temperature, °R (K) 

TCAN Array of external insulation container temperatures for conduction 
calculations, °R (K) 

TCAVG Average cooler tube temperature, # R in subroutine 

CYCL — or — temperature used in heat conduction calculation, °R in 
subroutine CNDCT (K) 

TCAVG1 Temperature used in heat conduction calculation, 8 R (K) 

TCAVG2 Temperature used in heat conduction calculation, *R (K) 

TCAVGC Average cooler tube temperature, °C 

TCAVGF Average cooler tube temperature, *F (*C) 

TCAVGK Average cooler tube temperature, "K 

TCGP Temperature used in appendix gap loss calculation at cold end of 
piston, °R (K) 

TCGPI Temperature of gas crossing interface between cold appendix gap and 
compression space, *R (K) 

TCLRM Inside wall cooler tube temperature, °R (K) 

TCLRMO Alternate storage location for TCLRM, °R (K) 

TCYL Array of cylinder housing temperatures used in conduction 
calculations °R (K) 

TG Array of control volume interface gas temperatures, °R (K) 

TGA Array of control volume gas temperatures, *R (K) 

TGACYC Array of time averaged control volume gas temperatures, °R (K) 

TGAO Alternate storage array for array TGA, °R (K) 

TGAVG Average temperature of piston walls, °R (K) 

TGCMPA Time averaged compression space gas temperature, *R (K) 

TGCYC Array of time averaged control volume interface gas temperatures, °R 

(K) 

TGEXPA Time averaged expansion space gas temperature, °R (K) 

THAVG Average heater tube temperature, °R (K) 

THAVGC Average heater tube temperature, °C 

THAVGF Average heater tube temperature, *F (°C) 
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THAVGK Average heater tube temperature, °K 

THCNDG Average thermal conductivity of gas in gap between piston and 
cylinder wall, Btu/in-sec-°R (W/cm-K) 

THCND1 Thermal conductivity, Btu/i n-sec-°R (W/cm-K) 

THCND2 Thermal conductivity, Btu/in-sec-°R (W/cm-K) 

THCOND Thermal conductivity, Btu/i n-sec-°R (W/cm-K) 

THGP Temperature used in appendix gap loss calculation, hot end of piston, 
°R (K) 

THGPI Temperature of gas crossing interface between hot appendix gap and 
expansion space, °R (K) 

TH20AV Average cooling water temperature, °R (K) 

TH20IN Cooling water inlet temperature, °R (K) 

TH20NC Cooling water inlet temperature, °C 

TH20NF Cooling water inlet temperature, “F (C) 

TH20NK Cooling water inlet temperature, K 

TIME Time since beginning of first engine cycle, sec 

TM Array of control volume wall temperatures, °R (K) 

TMA Array of control volume wall temperatures, °R (K) 

TMAO Alternate array for array TMA, °R (K) 

TMCYC Array of time averaged wall temperatures, °R (K) (current model lets 
only regenerator wall temperatures vary over the cycle) 

TMEXP Expansion space wall temperature, °R (K) 

TMHBR Back row heater tube outside wall temperature, °R (K) 

TMHFR Front row heater tube outside wall temperature, °R (K) 

TMIX Array of control volume gas temperatures, after mixing and before 

heat transfer, °R (K) 

TQ0FP4 Total heat out of engine (heat out through cooling water plus 
auxiliary loss), ft-lbf/cycle (J/cycle) 

TQ0HP4 Total heat out of engine (heat out through cooling water plus 
auxiliary loss), hp (kW) 

TQ0KW4 Total heat out of engine (heat out through cooling water plus 
auxiliary loss), kW 
TR Temperature ratio (out/in) 

TRAVG1 Average temperature, top half of regenerator housing, °R (K) 

TRAVG2 Average temperature, bottom half of regenerator housing, °R (K) 

TRIN Gas temperature at hot end of regenerator, °R (K) 
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TROUT Gas temperature at cold end of regenerator, °R (K) 

TRO Regenerator housing temperature, hot end, used in conduction 
calculation, *R (K) 

TR1 Regenerator housing temperature, middle, used in conduction 
calculation, °R (K) 

TR2 Regenerator housing temperature, cold end, used in conduction 
calculation, °R (K) 

TWOPI Constant, 2 * PI 

UTOTAL Internal energy content of working space gas control volumes, ft-lbf 

(J) 

O 

V Array of working space gas control volumes, in'-ycyl inder 

O 

(cm /cylinder) 

VAR Array of variables equivalent to the variable in COMMON /CYC/ 

VCAPGP Appendix gap volume at cold end of piston, in^ (cm^) 

VCLC Net compression space clearance volume (includes cold appendix 

o o 

volume), in (cnr) 

VCLE Net expansion space clearance volume (includes hot appendix gap 
volume), in^ (cm^) 

VC02 Volume fraction of carbon dioxide 

VEL Gas flow velocity, ft/sec (cm/sec) 

VELHD Velocity head, lbf/ft^ (N/cm^) 

VHAPGP Hot appendix gap volume, in^ 

VH2 Volume fraction of hydrogen 

VIS Array of control volume gas viscosities, Ibm/in-sec 

VISC Viscosity, lbm/ft-sec (kg/cm-sec) 

VISC02 Viscosity of carbon dioxide, Ibm/in-sec (kg/cm-sec) 

VISH2 Viscosity of hydrogen, Ibm/in-sec (kg/cn>-sec) 

VISH20 Average cooling water viscosity, lbm/ft-sec (kg/cm-sec) 

VISX Viscosity, Ibm/in-sec (kg/cm-sec) 

VO Alternative storage array for working space gas control volumes, 

in^ (cm^) 

VOVRT Summation over the array of gas control volumes 

of — volume/(temperature*pressure ratio (i.e. FIK)), in 4 -lbf /lbm-’R 
(cm 4 -N/kg-K) 

VR Velocity ratio (out/in) 

O O 

VTOTL Total working space volume, in (cm ) 
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3 

VTOTLO Alternate storage location for total working space volume, in 
(cm 3 ) 

W Total working space gas inventory per cylinder, Ibm (kg) 

WALT Total work per cycle if pressure drops are calculated relative to 

compression space pressure, ft-lbf /cycle (J/cycle) 

WALTO Alternate storage location for WALT, ft-lbf/cycle 
WCDUCO, WC0UC1, WCDUC2 Compression space work per time increment if only 
cooler pressure drop is considered, three consecutive time 
increments, ft-lbf (J) 

WCDUEO, WCDUE1, WCDUE2 Compression space work per time increment if only end 
effects pressure drop is considered, three consecutive time 
increments, ft-lbf (J) 

WCDURO, WCDURl, WCDUR2 Compression space works per time increment if only 

pressure drop considered is that in the cold half of the regenerator, 
three consecutive time increments, ft-lbf (J) 

WCF Correction factor for gas inventory to get desired mean pressure, 
dimensionless 

WCMPN Negative compression space work per cycle, ft-lbf/cycle (J/cycle) 

WCMPP Positive compression space work per cycle, ft-lbf/cycle (J/cycle) 

WCPCO, WCPC1, WCPC2 Compression space work per time increment, three 
consecutive time increments, ft-lbf (J) 

WCPO, WCP1, WCP2 Compression space work per time increment assuming no 
pressure drop, three consecutive time increments, ft-lbf (J) 

WDTGA (change in control volume gas temperature * control volume gas 
inventory), *R-lbm (K-kg) 

WFCTR (dimensionless function of control volume heat transfer * control 
volume gas inventory), lbm (kg) 

WEALTO, WEALT1, WEALT2 Expansion space work per time increment if pressure 
drop is calculated relative to compression space pressure, three 
consecutive time increments, ft-lbf (J) 

WEDUEO, WEDUE1, WEDUE2 Expansion space work per time increment if only end 
effects pressure drop is considered, three consecutive time 
increments, ft-lbf (J) 

WEDUHO, WEDUH1, WEDUH2 Expansion space work per time increment if only heater 
pressure drop is considered, three consecutive time increments, 
ft-lbf (J) 
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WEDURO, WEDUR1, WEDUR2 Expansion space work per time increment is only 
regenerator pressure drop is considered, three consecutive time 
increments, ft-lbf (J) 

WEPEO, WEPE1, WEPE2 Expansion space work per time increment, three 
consecutive time increments, ft-lbf (J) 

WEPO, WEP1, WEP2 Expansion space work per time increment assuming no pressure 
drop, three consecutive time increments, ft-lbf (J) 

WEXPN Negative expansion space work per cycle, ft-lbf/cycle (J/cycle) 

WEXPP Positive expansion space work per cycle, ft-lbf/cycle (J/cycle) 

WG Array of control volume gas inventories, Ibm (kg) 

WGCGP Cold appendix gap inventory, Ibm (kg) 

WGCGPO Cold appendix gap inventory at time increment previous to current 

value, Ibm (kg) 

WGHGP Hot appendix gap inventory, Ibm (kg) 

WGHGPO Hot appendix gap inventory at time increment previous to current 
value, Ibm (kg) 

WGOLD Array of control volume gas inventories, at one time increment before 
current value, Ibm (kg) 

WINT Work integral function 

WLALT Work loss at expansion space (due to pressure drop) when pressure 

drop is calculated relative to reference pressure in compression 
space, ft-lbf cycle (J/cycle) 

WLALTO Alternate storage location for WLALT 

WLCMC Work loss at compression space due to cooler pressure drop, 
ft-lbf/cycle (J/cycle) 

WLCME Work loss at compression space due to end effects pressure drop, 
ft-lbf/cycle (J/cycle) 

WLCMR Work loss at compression space due to regenerator pressure drop, 
ft-lbf/cycle (J/cycle) 

WLEALT Work loss at expansion space due to pressure drop when the pressure 
drop is calculated relative to the compression space pressure, 
ft-lbf/cycle (J/cycle) 

WLEXE Work loss at expansion space due to end effects pressure drop, 
ft-lbf/cycle (J/cycle) 

WLEXH Work loss at expansion due heater pressure drop, ft-lbf/cycle 
(J/cycle) 
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WLEXR Work loss at expansion due to regenerator pressure drop, ft-lbf /cycle 
( J/cycle) 

WRKBAS Indicated work + pressure drop work loss, per cylinder, ft-lbf /cycle 
(J/cycle) 

WRKCMP Compression space work per cycle, ft-lbf /cycle (J/cycle) 

WRKEXP Expansion space work per cycle, ft-lbf/cycle (J/cycle) 

WRKLC Work loss at compression space due to cooler pressure drop, per 

cylinder, ft-lbf/cycle (J/cycle) 

WRKLCM Work loss at compression space due to pressure drop, ft-lbf/cycle 
(J/cycle) 

WRKLCO Alternate storage location for WRKLC, ft-lbf/cycle, (J/cycle) 

WRKLE Total work loss due to end effects pressure drop, ft-lbf/cycle 
(J/cycle) 

WRKLEO Alternate storage location for WRKLE, ft-lbf/cycle (J/cycle) 

WRKLEX Work loss at expansion space due to pressure drop, ft-lbf/cycle 
(J/cycle) 

WRKLH Work loss at expansion space per cylinder due to heater pressure 
drop, ft-lbf/cycle (J/cycle) 

WRKLHO Alternate storage location for WRKLH, ft-lbf/cycle (J/cycle) 

WRKLR Total work loss due to regenerator pressure drop, ft-lbf/cycle 
(J/cycle) 

WRKLRO Alternate storage location for WRKLR, ft-lbf/cycle (J/cycle) 

WRKLT Total work loss due to pressure drop, ft-lbf/cycle (J/cycle) 

WRKLTO Alternate storage location for WRKLT, ft-lbf/cycle (J/cycle) 

WRKTOT Indicated work per cycle, ft-lbf/cycle (J/cycle) 

WTPCO, WTPC1, WTPC2 Total work per time increment when pressure drop is 
calculated relative to compression space pressure, ft-lbf/cycle 
(J/cycle) 

WO Array of variables equivalent to works in COMMON /TIMEO/ 

W1 Array of variables equivalent to works in COMMON /TIME 1 / 

W2 Array of variables equivalent to works in COMMON /TIME2/ 

X Array (two dim.) of piston positions, in (cm) 

XC02 Mass fraction of carbon dioxide 

XCPA Specific heat at constant pressure, Btu/lbm-°R (J/kg-K) 

XCV Specific heat at constant volume, Btu/lbm- # R (J/kg-k) 

XGAM Ratio of specific heats (CP/CV) 
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XH2 Mass fraction of hydrogen 

XL Array of control volume flow lengths, in (cm) 

XLG Length, ft (cm) 

XLGTH Length, in (cm) 

XLO Alternate storage array for XL, in (cm) 

XMAO Estimate of outlet Mach number 

XMA1 Inlet Mach number 

XMA2 Outlet Mach number 

ZERO Constant =0.0 

ZMC02 Molecular wt. of carbon dioxide 

ZMH2 Molecular wt. of hydrogen 

ZMMIX Molecular weight of mixture of hydrogen and carbon dioxide 
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APPENDIX G: COMPARISON OF PREDICTIONS WITH TEST DATA 


Predicted P-40 engine brake power and efficiences are compared with the 
results of engine tests made at Lewis Research Center in figure 15. The 
tests were made with auxiliaries powered by the engine. The efficiencies 
shown are overall efficiencies. The efficiency predicted by the computer 
program does not account for the combustor efficiency. Thus it was neces- 
sary to use an assumed combustor efficiency to adjust the predictions of the 
computer program. The combustor efficiencies calculated from the Lewis P-40 
test data were all about 90 percent for the test points shown. When the 
predicted efficiencies were multiplied by 0.90, the upper predicted 
efficiency curve was obtained . However, information obtained from United 
Stirling suggests the P-40 combustor efficiency may be closer to 80 percent 
for the range of operation shown. When the predicted efficiencies were 
multiplied by 0.80, the lower predicted efficiency curve was obtained. 

The regenerator effectiveness (average of REFF1 and REFF2 - defined in 
the symbols list) was about 0.996 for the predictions of figure 15. When 
the computer program was modified to yield a regenerator effectiveness of 
about 0.990 (by multiplying DTGASL by 0.96, in subroutine HEATX) the pre- 
dictions were as shown in figure 16. 
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TABLE I. - EFFECT OF CHANGING THE NUMBER OF 
CONTROL VOLUMES IN THE HEAT EXCHANGERS 


CONTROL VOLUME 
CHANGE 


CHANGE IN TOTAL 
NO. OF CONTROL VOLUMES 

EFFECT ON 
ENGINE POWER 

REGENERATOR, 5 » 

7 

17 » 19 

« + 0.5 kW 

HEATER, 3 > 5 


17 » 19 

= + 0.25 kW 

COOLER, 3-5 


17 » 19 

= + 0.25 kW 

REGENERATOR, 5 > 
HEATER, 3 » 5 
COOLER, 3 > 5 

7 

17 » 23 

= + 1.0 kW 
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TABLE II. - INPUT DATA— ENGINE PARAMETERS 


8ENGINE EID=8HP40 ,ETYPE=1., DISPD=2.165, DISPRD= . 4724 , DSPGAP=. 01575, 

DSPHGT=3.53, RODL=3.937, RCRANK= . 7874 , E=0.0, PHASE=90., 

HTBOD= . 1772 , HTBID=.1181, HTBPCN=18., HTBL=11.05, EHTBL=9.913, 

REGPCN=2 . , REGID=2 . 244 , REGL=1.535, RWIRED=1 . 969E-3, 

PROSTY=0 . 58 » RMDEN=0 . 282, CPM=0.11, CTB0D = 5 . 906 E-2 , 

CTBID=3. 937E-2,CTBPCN=384. , CTBL=3.150, ECTBL=2.646, 

CNDSS=2.778E-3,CPH20=1.0, RH0H20=62.4, AEFH20=5.00, EXPSCL= . 1452, 

EXPHDV= . 2691, HRDV=1.017, RCDV=.3112, CCNPDV=1 . 843 , CMPSCL=. 1379, 

CYLORT=l . 36 , CYL0RM=1.31, CYL0RB=1.25, CYLDTM=0 . 787 , CYLDMB=0 .787 , 

DSPWTH=0 . 075 , REG0RT=1.42, REGORM=l . 37 , REGORB=1.30, REGDTM=0 . 594 , 

REGDMB=0 . 594 , STROKE=l . 575 , CANOR=0.0, CANIR=0.0, CONDTB=0.0, 

DFREQ=66 . 67 , DFLOSS = 17 . 165 , DALOSS = 10.06 

SEND 
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SYMBOL 

EID 

ETYPE 

DISPD 

DISPRD 

DSPGAP 

DSPHGT 

RODL 

RCRANK 

E 

PHASE 

HTBOD 

HTBID 

HTBPCN 

HTBL 

EHTBL 

REGPCN 

REGID 

REGL 

RWIRED 

PROSTY 

RMDEN 

CPM 

CTBOD 

CTBID 

CTBPCN 

CTBL 

ECTBL 

CNDSS 


TABLE III: SYMBOL DEFINITIONS FOR 

ENGINE PARAMETER INPUT DATA 
(NAMELIST /ENGINE/) 

DEFINITION 

Alphanumeric engine identifier 
Numeric engine identifier 
Piston diameter, in (cm) 

Piston rod diameter, in (cm) 

Piston-cylinder gap, in (cm) 

Displacer height, in (cm) 

Connecting rod length, in (cm) 

Crank radius, in (cm) 

Eccentricity (not used) 

Angle by which compression volume lags expansion 
volume, deg 

Heater tube outside diameter, in (cm) 

Heater tube inside diameter, in (cm) 

Number of heater tubes per cylinder 
Heater tube length, in (cm) 

Length of heater tube effective in heat transfer, in (cm) 
Number of regenerators per cylinder 
Regenerator inside diameter, in (cm) 

Regenerator matrix length, in (cm) 

Regenerator matrix wire diameter, in (cm) 

Regenerator matrix porosity 

Regenerator matrix metal density, — ) 

in J \cnT / 


Regenerator matrix specific heat, 

lb m 

Cooler tube outside diameter, in (cm) 

Cooler tube inside diameter, in (cm) 

Number of cooler tubes per cylinder 
Cooler tube length, in (cm) 

Length of cooler tube effective in heat transfer, in (cm) 
Cooler tube thermal conductivity. 
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CPH20 


RH0H20 


AEFH20 

EXPSCL 

EXPHDV 

HRDV 

RCDV 

CCMPDV 

CMPSCL 

CYLORT 

CYLORM 

CYLORB 

CYLDTM 

CYLDMB 

DSPWTH 

REGORT 

REGORM 

REGORB 

REGDTM 

REGDMB 

STROKE 

CANOR 


Btu f - watt j 

ft - sec - °R 'em - k / 

Cooling water specific heat, 

Btu l joules \ 

1 bm - °R \ gm - k / 

Density of water, — 3 /^ 3 ^ 

ft y cm J 


9 

Effective cooling water flow area per cylinder, in 
(cm 2 ) 

Q O 

Expansion space clearance volume, in (cm ) 

Expansion space - heater connecting duct volume, irr 
(cm 3 ) 

Heater-regenerator connecting duct volume, in J (cm 0 ) 
Regenerator-cooler connecting duct volume, in 3 (cm 3 ) 
Cooler-compression space connecting duct volume, in 
(cm 3 ) 

O O 

Compression space clearance volume, in (cm ) 

Cylinder housing outside radius, top, in (cm) 

Cylinder housing outside radius, middle, in (cm) 

Cylinder housing outside radius, bottom, in (cm) 

Cylinder housing conduction length, top to middle, in (cm) 
Cylinder housing conduction length, middle to bottom, in 
(cm) 

Piston wall thickness, in (cm) 

Regenerator housing outside radius, top, in (cm) 

Regenerator housing outside radius, middle, in (cm) 
Regenerator housing outside radius, bottom, in (cm) 
Regenerator housing conduction length, top to middle, in 
(cm) 

Regenerator housing conduction length, middle to bottom, in 
(cm) 

Piston stroke, in (cm) 

Insulation container outside radius (not used) 
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CANIR 
CONDTB 
DFREQ 
DFLOSS 
DA LOSS 


Insulation container inside radius (not used) 
Insulation container conduction length (not used) 
Design engine frequency, Hz (rpm) 

Design mechanical friction loss, hp (kW) 

Design auxiliary power requirement, hp (kW) 
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TABLE IV. - MODEL OPTION SWITCHES AND MULTIPLYING FACTORS 
AND ENGINE OPERATING CONDITIONS 


8STRLNG REALGS=1 . , FACT1=0.4,FACT2=10.0, 

NOCYC=25,NSTRT=1,NOEND=20,MWGAS=2, 

RUCFAC=1.0,HHCFAC=1.0,CHCFAC=1.0,IPCV=0,FMULT=1.0,FMULTR=1.0, 

INIX=0, VH2=0 . 99, IPUMP=1 , ICOND=0 , 

IOUT=1,JIP=0,IPRINT=500,ITMPS=0,MAPLOT=1 SEND 

8IHDATA IDRUN = 12HREAD S123BR , P = 2 17 1 . , 0MEGA = 66 . 78 , TMEXP=16 A3 . , 

TMHFR=1672. ,TMHBR=1672. ,TCYL=1643. ,1AA1. ,1239. ,TCAN=1191. ,999. ,TR0=1A72. , 
TR1=121A. ,TR2=956. , GPMH20=13 . 63 , TH20IN=580 . 1 SEND 
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TABLE VII. - PARTIAL LISTING OF LONG FORM OF OUTPUT 
(PRODUCED WHEN IOUT = 1, JIP = 0) 

* RUN IDENTIFICATION XX 
READ ttl23BR , 


CALCULATED CONTROL VOLUME AND ENGINE PARAMETFJ G , AND INPUT DATA *x 
X FT ARAM 

1 ■'‘RDV- 2 1733355358357 
r DrDV= 7 0 A 20762338519 
CEEDV- 1 6725263637283 
DGAPDV = 0 3S090003926925 

DI’O = 2X2 1650 , 3X0.11310 , 2.2660, 5X0 .27190952380952D-02, 0.22660D01, 3X0 . 3937 0D-01 
2-0 169260D01, 13X0.0 

ACSO= 3 6813379069113, 1.8606689536556, 3x0 . 98539G65512928D-01 , 0 . 39569007266381D01 
5-0 22933626216501D01, 0 . 3956 90 07266381D0 1 , 3X0 . 27373602598861D0, 0 . 17530335352219D01 
0 35060b70706639D01, 13X0.0 

AHTO = 2-0 0, 3X11.033826006996, 0.0, 5X1035.9A69775373, 0.0, 3X20.9A5228AA0128 
15-0.0 

XLO= 0 0,0 73098626215699D-01, 3X0 . 36833333333333D01 , 0 . 12857665588833D0 

5-0 3 0 7 0 D 0 , 0.39363591850931D-01, 3X0.1050D01, 0 . 52566 022356121D0 , 16X0.0 

VO- 0.0, 0.26910, 3X0.72627852861190, 1.0170, 5X1.6086192667706, 0.31120 

3^0.69036165657609, 1 3630, 16X0.0 

V C L E - 0 52610003926925 

VCLC = 0 17599000392693 

CLRLOD- 80 010160020320 

1 1 T R L 0 D - 93.566775613887 

PCAREA= 3.9569007266381 

tV- 3 6813379069113 

AR = 0.17527083666763 

APMAR- 3 5060670706639 

VH APGP = 0.38090003926925 

'<CAPGP= 0 . 33090003926925D-01 

THGP = 1373 6666666667 

TCGP = 666.60896661755 

qiigps= o.o 

qcgps= o.o 

iruru>= l 

icond= o 

X END 


xx INPUT DATA — ENGINE PARAMETERS XX 
2 ENGINE 

EID= -0 6 7377992506120 D28 

E7YPE= 1.0 

DI5rD= 2.1650 

DIGrRD= 0.67260 

D5PGAP- 0 1 57 50D-0 1 

DGP1IGT = 3.530 

RODL = 3.9370 



rCRANK= 0 . 78740 
E = 0 0 
r»IASE= 90.0 
illPOD = 0 17720 

ii r a i d= o.n8io 

IITBPCM = 18.0 
HT3L = 11.050 
EI!T3L = 9.9130 
REGPCN= 2.0 
REGID= 2.2940 
F’.GL= 1 5550 
RWIRED= 0 19690D-02 
rR05TY= 0.580 
RMDEN= 0 2320 

crn= o no 

C f BOD- 0.59060D-01 
CTBID= 0.39370D-01 
CT 3PCN = 384.0 
CTBL= 3.150 
ECTBL = 2.6460 
CND 35= 0 27780D-02 
CPH20= 1.0 
RII0H?0 = 62.40 
AEHI20 = 5.0 
E. V PSCL = 0.14520 
EXrtlDV= 0.26910 
l!RDV = 1.0170 
P.CDV = 0.31120 
CCr;PDV= 1 8430 
CmPSCL = 0.13790 
CYLORT = 1.360 
CYL ORM= 1.310 
CYLORB= 1.250 
CYL DT!1= 0.7870 
CYLDI-!B= 0 7870 
psruni: 0 . 750D-01 
REGORT = 1 420 
RECORM= 1.370 
REC0RB= 1.30 
RCGDTM= 0.5940 
RCGBMB= 0.5940 
E T R0KE = 1.5750 
CAHOR= 0.0 
CANIR= 0.0 
C0NDTB= 0.0 
DPREQ= 66.670 
DTLOSS= 17.1650 
DALOES= 10.060 
l END 


XV INPUT DATA— OPTION SWITCHES AND MULTIPLYING FACTORS XX 


3STRLNG 



k 


REALGS= 1.0 
FACT 1 - 0.40 
TACT2= 10 0 
liCCYC= 25 
NSTRT= 1 
N0CND= 20 
IHJGAS= 2 
RIICF AC= 1.0 
HHCrAC= 1.0 
CHCrAC= 1.0 
IPCV = 0 
FMUL T = 1.0 
FMULTR= 1.0 
K1IX = 0 
VH2= 0 990 
IPUMP= 1 
ICOIiD= 0 
IOUT= 1 
J IP- 0 
IPPINT = 500 
ITMPS= 0 
MAF LOT = 1 
3FPD 


'X INPUT DATA— ENGINE OPERATING CONDITIONS xx 
SINDATA 

IDR'JN= -641351223, 1081864690, -205334165 
r= 2635.5934534073 
ui OriEGA= 66.780 
TMEXP= 1643 0 
TI1UFR= 1672.0 
TMHBR= 1672.0 

TCYL- 1643.0, 1441.0, 1239.0 

TCAN= 1191 0, 999.0 

T RO - 1472.0 

TP.1= 1214.0 

TR2= 956.0 

GPMII20 - 13.630 

TII20IN= 580.10 

J END 


xx CYCLE BY CYCLE SUMMARY DATA XX 


x CYCLE NO. 1 x 

TIME ANGLE X(l) X(2) 
0 0000 0.0 0.00 0.71 

0 0150 360.0 0.00 0.71 

TC-- 1542. 1580. 

TGA= 1544. 1543. 1581 

TM= 1643. 1672. 1672 


TRIN TROUT T1 

0.00 0.00 1709. 

1491.76 718.73 1544. 

1645. 1492. 

1648. 1643. 1413. 

1672. 1672. 1413. 


T2 PE PC PRIN 

801. 2635.59 2635.59 0.00 

737. 2601.96 2601.96 2601.96 

719. 704. 

788. 718. 703. 716 

789. 646. 646. 646 


PROUT FRIN FROUT 
0.00 0.00000 0.00000 
2601.96-0.05553-0.23852 
734. 737. 

734. 737. 

646. 696. 



t 


0.380 REFF1= 0.971 REFF2= 


QIN = 497 260 QOUT = 386.541 WRKEXP= 407.993 WRKCMP =-21 9 . 0 57 WRKTOT= 188.937 EFFTOT= 

QRCGEU= 54.846 PKRHP= 22.940 FP.EQ= 66.780 WRKLT= 0.000 UTOTAL= 0.000 
EN3Pn= 6310.710 AVGWSP= 2305.225 „ „„„ 

7EXPN= -6.329 QEXPP= 0.021 QHEATN=-490 . 952 QHE.‘TP = 0.000 QCOOLN= -0.131 QCOOLP= 333.692 QCOMPN= -0.570 QCOMPP- 

QREGN=-1805. 030 QREGP = 1859 .876 QCNDRI= 0.000 QCNERC- 0.000 

QCNDCL = 0 . 000 QCNDD= 0.000 QCNDCN= 0.000 TGEaPm = 14 98 . 87 3 T GCMP A= 736.856 QSHTL= 0.000 

QINB = 497 260 WRKBAS = 138.937 WALT = 188.937 QOUTa= 334.777 WRKLEX= 0.000 WRKLCM= -0.000 WLEXR= 0.000 


0.963 


' >J L C f* 1 R = 


-0.000 WL EXE= 0.000 WLCME= -0.000 l.'RKLK- 


0.000 WRKLH : 


0.000 WRKLC= 


0.000 WRKL E= 0.000 WLALT= 


0.000 


» CYCL 
TIME 
0 015 
0 . 029 
TG = 
TCA = 

TM = 

QIN = 4 
CREGEN 
EN3PM= 
QEXPN= 
QREGH= 
QCNDCL 
QINB = 
WLCMR= 

< CYCL 
line 
0 029 
0 04 4 
, TG= 

T GA = 
T!1= 

QIN = 4 
QREGEN 
EN3TM= 
Q EXPN= 
QRESN= 
QCHDCL 
QINB = 
ULt MR = 


NO. 2 X 
ANGLE 
360.0 
360.0 
1569. 
1571. 
1643. 


PE PC PRIN PROUT FRIN FROUT 
2601.96 2601.96 2601.96 2601.96-0.05553-0.23852 
2583.98 2583.98 2583.98 2583.98-0.05611-0.23860 
718. 701. 718. 714. 

717. 700. 705. 718. 714 

646. 646. 646. 646. 696 

0.447 REFF1= 


X(l) X( 2 ) TRIN TROUT T1 T2 

0.00 0.71 1491.76 718.73 1544. 737. 

0.0 0 0.71 1489.0 7 717.78 15' r l 714. 

1603 1658. 1489. 

1570. 1603 1663. 1656. 1411. 794 

1672. 1672. 1672. 1672. 1411. 795 

21.420 QOUT = 332.610 WRKEXP= 402.730 WRKCMP =-21 4 . 44 5 WRKTOT= 188.285 EFFTOT = 

= 2.964 PURHP = 22.861 FREQ= 66.780 WRKLT = 0.000 UTOT AL= 0.000 

-5!259 1 QEXPP = r ' P Q 2 003 ’ QHEATN=-416 . 164 0HE/.TP= 0.000 QCOOLN= -0.665 QCOOLP= 283.912 QCOMPN= -1.268 QCOMPP : 
-1824.564 QREGP=1827 .528 QCNDRI= 0.000 CCNDRO= 0.000 

0.000 QCNDD= 0 . 00 0 QCNDCN= 0.000 TGEXPA = 1515 . 042 TGCMPA= 696.121 QSHTL = 

421 420 WRKBAS = 188.285 WALT = 188.285 Q0UT5= 232.702 WRKL EX= 0.000 WRKLCM = 


0.991 REFF2= 0.978 


0.000 

-0.000 WLEXR= 


0.000 


-0.000 WL EXE= 0.000 WLCME= -0 . 000 WP.KLR= 0.000 WRKLH= 0.000 WRKLC= 0.000 WRKLE- 0.000 WLALT= 0.000 


NO. 3 x 

ANGLE X(l) X(2) 
360.0 0.00 0.71 

360.0 0.00 0.71 

1578. 1609. 

1579. 1578. 1610. 

1643. 1672. 167; 


PE PC PRIN PROUT FRIN FROUT 
2583.98 2583.98 2583.98 2583.98-0.05611-0.23860 


2573.47 2573.47 2573.47 
717. 700. 

716. 699. 700 

646. 646. 646 


TRIN TROUT T1 T2 

1489.07 717.78 1571. 714. 

1486.70 717.41 1579. 703. 

1662. 1487. 

1668. 1660. 1409. 797 

1672. 1672. 1409. 798. .... 

06 087 QOU T = 308.442 WRKEXP= 398. 8S2 WRKCMP = -2 1 1 . 725 WRKTOT= 187.157 EFFTOT= 0.461 REFF1 - 0.996 REFF2= 0.981 

= -16.195 PWRHP= 22.724 FREQ= 66.780 WRKLT= 0.000 UTOTAL= 0.000 

-^SO^QEXPP^ 215 0^074 QH E ATN = -4 02 . 1 0 1 Q H E A T P = 0.838 QCOOLN= -1.336 QCOOLP= 259.344 QCOMPN= -1.697 QCOMPP = 

-1825 480 QREGP = 18 0 9 . 285 QCNDRI= 0.000 QCNDRO= 0.000 

= 0.000 QCNDD= 0.000 QCNDCN= 0.000 TGEXP A = 1 526 . 426 TGCMPA= 680.010 QSHTL= 0.000 

406.087 WRKBAS = 187.157 WALT = 187.157 QOUTB= 256.701 WRKLEX= 0.000 URKLCM= -0.000 WLEXR= 


2573.47-0.05629-0.23836 
710. 703. 

710. 703. 

646. 696. 

0.461 REFF1= 0 


-0 000 


WL EXE= 0 . 000 WLCME= -0.000 WP.KLP.= 0.000 WRKLH= 0.000 WRKLC= 0.000 WRKLE= 


0.000 

0.000 WLALT= 


0.000 


‘ CYCLE NO 4 X 

TIME ANGLE XC 1 J X(2) TRIN TROUT T1 T2 

0 0449 360.0 0.00 0.71 1486.70 717.41 1579. 703. 


PE PC PRIN PROUT FRIN FROUT 
2573.47 2573.47 2573.47 2573.47-0.05629-0.23836 


x/ METAL TEMPERATURES FOR CONDUCTION CALCULATIONS XX 
2 IHTEMP 

T CT L = 1643 0, 1441.0, 

" T C AN= 1191.0, 999.0 
TRO= 1472 0 
TR1= 1214.0 
7 R2= 956.0 
SEND 


.786 


.722 


.389 


1239.0 



0.0599 360.0 0.00 0.71 1484.87 716.85 1581. 698. 2562.99 2575.32 2563.88 2569.98-0.05633-0.23828 

TG= 1580. 1611. 1663. 1485. 717. 699. 706. 698. 

TO A" 1581. 1580. 1612. 1669. 1661. 1403. 797. 716. 698. 697. 706. 698. 

Ill- 1643. 1672. 1672. 1672. 1672. 140S. 799. 646. 646. 646. 646. 696. 

9IN = 390.941 QOUT = 308.637 WRKEXP= 376.737 WRKCMP=-215 . 681 WRKTOT= 161.056 EFFTOT= 0.412 REFF1= 0.997 REFF2= 0.984 
CRFGEH= -17.064 PWRHP= 19.555 FREQ= 66.780 URKLT= 25.379 UTOTAL= 0.000 
CN3FM- 6310.710 A VGIJSP = 2269.566 

CEXPtt= -4.793 OEXPP= 0.105 QHEATN=-398 . 994 QHEATP= 1.619 QCOOLN= -1.934 QC0OLP= 245.702 QCOMPN= -1.926 QCOMPP= 0.251 
QPEGN =-1819. 351 QP.EGP = 1802 . 287 QCNDRI= 5.178 QCHDRO= 5.178 

CCNDCL = 1.3S3 QCNDD= 0.353 QCNDCN= 0.000 TGEXPA=1529 . 822 TGCMPA= 672.399 QSHTL= 1.939 

OINB = 402.063 URKBAS = 186.435 WALT = 183.283 QOUTB= 242.092 WRKLEX= 19.981 WRKLCM= 5.398 WLEXR= 3.767 

ULCUR= 2.663 WLEXE= 9.274 WLCME= 1.312 URKLR= 6.430 WRKLH= 6.940 WRKLC= 1.423 WRKLE= 10.587 WLALT= 22.227 

-< CYCLE NO 5 x 

TINE ANGLE X(l) X(2) TRIN TROUT T1 T2 PE PC PRIN PROUT FRIN FROUT 

0 0599 360.0 0.00 0.71 1484.87 716.85 1581. 698. 2562.99 2575.32 2563.88 2569.98-0.05633-0.23828 

0.0749 360.0 0.00 0.71 1494.37 721.81 1585. 697. 2568.83 2568.83 2568.83 2568.83-0.05616-0.23831 

TG= 1584. 1615. 1666. 1494. 722. 705. 707. 697. 

TGA = 1585. 1584. 1615. 1671. 1664. 1416. 794. 721. 704. 700. 707. 697. 

TI1 = 1643. 1672. 1672. 1672. 1672. 1416. 796. 646. 656. 656. 646. 594. 

QIN = 398.079 QOUT= 256.862 WRKEXP= 395.727 WRKCMP=-209 . 382 WRKTOT = 160.966 EFFTOT = 0.404 REFF1= 1.003 REFF2= 1.000 

QREGEN= -14.004 PWRHP= 19.544 FREQ= 66.780 WRKLT= 25.379 UTOTAL= 0.000 
CN3PM= 6310.710 AVC-WSP = 2269.482 

QLXP N= -4.668 QEXPP= 0.124 QHEATN=-386 . 7 91 QI1EATP= 2.115 QCOOLN= -5.736 QCOOLP= 197.080 QCOMPN= 0.000 QCOMPP= 3.533 
QRCGN=-1S29.721 QREGP =18 1 5 . 7 1 7 QCNDRI= 5.178 QCNDRO= 5.178 

CCNCCL= 1.388 QCNDD= 0.353 QCNDCN= 0.000 TGEXPA=1533 . 308 TCCMPA= 669.501 QSHTL= 1.939 

OINB = 389.221 URKBAS = 186.345 WALT = 186.345 Q0UTB= 194.877 WRKL EX= 0.000 WRKLCM= 0.008 WLEXR= 0.000 

WLCnR= 0.003 WLEXE= 0.000 WLCME= 0.003 WRKLR= 6.430 WRKLH= 6.940 WRKLC= 1.423 WRKLE= 10.587 WLALT= 0.000 


AVERAGE VALUES OVER LAST 5 CYCLES *x 

OIN= 422.757 QOUT= 318.618 WRKEXP= 396.414 WRKCMP=-214 . 058 WRKTOT= 177.280 EFFTOT= 0.421 REFF1= 0.991 REFF2= 0.981 
QREGEN= 2.109 PWRHP= 21.525 FREQ= 66.780 WRKLT= 10.152 UTOT AL= 0.000 
EN3PM= 6310.710 AVGW5P= 2281.762 

CEXrN= -5.189 QEXPP= 0.065 QHEATN = -419. 000 QHEATP = 0.914 QCOOLN= -1.960 QCOOLP= 263.946 QCOMPN= -1.092 QCOMPP= 1.336 

QREGN=-1820 .829 QREGP=1822 . 939 QCNDRI = 2.071 QCNDRO= 2.071 

OCNDCL = 0.555 QCNDD= 0.141 QCNDCN= 0.000 TGEXPA = 1520 . 694 TGCMPA= 690.978 QSHTL= 0.776 

Q I N 3 = 423.210 URKBAS = 187.432 WALT = 186.801 Q0UTB= 262.230 WRKLEX= 3. 996 WRKLCM= 1.081 WLEXR= 0.753 

WLCMR= 0.533 WLEXE= 1.855 WLCME= 0.263 WRKLR= 2.572 WRKLH= 2.776 URKLC= 0.569 WRKLE= 4.235 WLALT= 4.445 


xx AVG. TEMPS.— LAST 5 CYCLES, AVG. H. T. COEFS. & HEAT FLUXES— LAST CYCLE XX 


TGCYC= 1526. 1524. 1524. 1485 713. 702. 683. 673. 


TGACYC= 

1533. 

1523. 

1542. 

1542. 

1506 . 

1408. 

790. 

708. 

691. 

681. 

678. 

670 . 

ti;cyc= 

1643. 

1672. 

1672. 

1672. 

1672. 

1407 . 

790 . 

646 . 

656. 

656 . 

646 . 

594. 

H A C Y C = 

0.057 

0.000 

0.835 

0.836 

0.000 

1 .965 

1.924 

0.000 

0.640 

0.689 

0.000 

0.077 

HMX = 

0 401 

0.000 

1.467 

1.488 

0.000 

5.034 

2.848 

0.000 

1.054 

1.124 

0.000 

0.248 

QOAAVG= 

6.00 

0.00 

72.41 

74.46 

0.00 

4.55 

4.46 

0.00 

21.83 

19.25 

0.00 

6.02 

OOAMX= 

24.19 

0.00 

140.12 

127.07 

0.00 

10.33 

8.71 

0.00 

46.45 

59.67 

0.00 

29.20 



< CYCLE 
TIME 
0 0 7 4 9 
0 0398 
TG- 
TCA = 

IM = 


NO. 6 * 
ANGLE 
360 0 
360.0 
1585. 

1585 
1643. 


XC1) X( 2 ) TRIN TROUT T1 T2 

0.00 0.71 1494.37 721 81 1585. 697. 

0.00 0.71 1490.89 721.98 1585. 697. 

1615. 1666. 1491. 

1585. 1615. 1671. 1664. 1413 795. 721. 

1672. 1672 1672. 1672. 1413. 796. 646. 

QIN= 399 382 QOUT= 260.451 MRKEXP= 395.311 MRKCMP =-20 9 . 6 1 5 MRKTOT= 160 
QRlGEN= -4.931 PWRHP= 19.465 FREQ= 66.780 MRKLT= 25.379 UTOTAL = 

E H 3 r H - 6310.710 AVGU3P= 2268.878 

Qf*>TN = -4.664 QEXPP: 0.137 QHEATN=-388 .539 QHEATP= 2.543 QCOOLN = 

QREGN=-1317 369 QREGP=1812 . 433 QCNDRI = 5.178 QCNPRO= 5.178 


PE PC PRIN PROUT FRIN FROUT 
2568.83 2568.83 2568 83 2568.83-0.05616-0.23831 
2563.12 2563.12 2568.12 2568.12-0.05622-0.23828 
722. 705. 707. 697. 


704. 700. 

656. 656. 

317 EFFTOT= 

0.000 


707. 697. 

646. 594. 

0.401 REFF1= 1. 


001 REFF2= 


1 .000 


-5.600 QCOOLP= 200.579 QCOMPN= 0.000 QCOMPP= 


:ndcl= 


1.338 QCNDD= 


0.353 OCNDCN= 0.000 TGEXPA = 1534 . 17 9 T GCMPA = 


QINB = 390.523 MRKBAS = 

185. 

,696 MALT 

= 185.696 

QOUT B 

\'LCMR= -0 000 MLEXE = 

0 000 MLCME= 

-0.000 MR 

KLR = 

* CYCLE NO 7 X 

TIME ANGLE X(l) 

X( 2 ) 

TRIN 

TROUT 

T 1 

0 0898 360.0 0.00 

0.71 

1490.89 

721.98 

1585. 

0 1048 360.0 0.00 

0.71 

1488.24 

722.23 

1585. 

TG= 1584. 1614 


1665. 1488 


TGA= 1585 1584. 

1615. 

1671. 

1663. 

1411 . 

TM= 1643. 1672. 

1672. 

1672 

1672. 

1411. 


668.912 QSHTL= 1.939 
0.000 MRKLCM= -0.000 MLEXR= 


0.000 


6.430 URKLH= 6.940 MRKLC= 1.423 WRKLE= 10.587 MLALT- 


0.000 


T 2 PE PC 

697. 2568.12 2568.12 

697. 2564 18 2576.51 

722. 705. 

795. 721 704 

796. 646. 656 


Q I N = 382.149 QOUT = 267.250 MRKEXP= 374.903 MRKCMP=-215 . 062 MRKTOT= 159.840 


PRIN 
2568 . 12 
2565.07 

701 

656 

EFFTOT= 


PROUT FRIN 
2568.12-0.05622-0 
2571.16-0.05626-0 
707. 697. 

707. 697. 

646. 594. 

0.418 REFF1 = 1 


FROUT 

.23828 

.23824 


001 REFF2= 0.999 


QRCGEN= -4 851 PWRHP= 19 408 FREQ = 66.780 MRKLT= 25.350 UTOTAL- 

ctr>rn= 6310.710 avgmsp= 2268.445 

Q EXPN = -4 680 QEXPP= 0.136 QHEATN=-391 . 21 0 QHEATP= 2.510 QCOOLN = 

9FCGN=-1812.870 QREGP=1S08. 019 QCNDRI = 5.178 QCNDRO= 5.178 

QCNDCL= 1.388 QCNDD= 0.361 QCNDCN= 0 000 TGEXPA=1 5 33 . 7 92 TGCMPA= 
QINB = 393.243 MRKBAS = 135.190 MALT = 182.044 QOUTB= 199 . 907 MP.KLEX = 
ML CT1R= 2 649 ML EXE = 9.261 MLCME= 1.317 MRKLR= 


0.000 


-5.469 QCOOLP= 201.870 QCOMPN= 0.000 QCOMPP= 


668.810 QSHTL= 
19.961 MRKLCM= 


1.939 
5.389 MLEXR= 


3.768 


6 . 4 i 7 MRKLH= 6.931 MRKLC= 1.424 MP.KLE= 10.578 MLALT= 22.204 


•< CYCLE 
TIME 
0 1048 
0.1193 
TG = 

T G A - 
TM = 

QIN= 40 
QRLGEN= 
ei; 3 r n= 

QEXPN= 
QPEGN=- 
QCNDCL= 
QINB = 
MLCI1R = 


NO 8 
ANGLE 
360.0 
360.0 
1582 
1583. 
1643. 


XC1) X(2) TRIN TROUT T1 

0.00 0.71 1488.24 722.23 1535 

0.00 0.71 1491.98 708.23 1583 

1613. 1665. 1492. 

1582. 1613. 1670. 1663. 1413 

1672. 1672. 1672. 1672. 1413 


T2 PE PC 

697. 2564.18 2576.51 

689. 2543.65 2543.65 

708. 692. 

783. 707. 691, 

785. 646. 642, 


PRIN PROUT FRIN 
2565.07 2571.16-0.05626- 
2543.65 2543.65-0.05597- 
697. 689. 

689. 697. 689. 

642. 646. 592. 


FROUT 

.23824 

.23882 


0.120 QOUT = 288.393 URKEXP= 397.424 MRKCMP = -2 07 . 7 5 9 MRKTOT= 164.316 EFFTOT= 0.411 REFF1= 0.996 REFF2= 0.999 


45.184 PWRHP= 19.951 FREQ= 66.780 WRKLT= 25.350 UTOTAL= 
6310.710 AVGMSP= 2253 392 

-4.692 QEXPP= 0.124 QHEATN=-388 . 856 QHEATP= 2.170 QCOOLN= 
1809.218 QREGP=1854 402 QCNDRI= 5.178 QCNDRO= 5.178 

1.388 QCNDD= 0.361 QCNDCN= 0.000 TGEXP A =1532 . 46 4 TGCMPA= 
391 254 MRKBAS = 189.665 MALT = 189.665 QOUTB= 228.034 WRKL EX = 
0.003 ML EXE= 0.000 MLCME= 0.003 MRKLR= 6.417 MRKLH = 


0.000 

-3.186 QCOOLP= 227.849 QCOMPN= 0.000 QCOMPP= 


664.988 QSHTL= 
0.000 MRKLCM= 


1 . 939 
0.008 MLEXR= 


0.000 


6.931 MRKLC= 1.424 MRKLE= 10.578 MLALT= 0.000 


* CYCLE NO. 9 * 


TIME 

ANGLE 

XC1) 

X(2) 

TRIN TROUT 

T 1 

0 1198 

360.0 

0.00 

0.71 

1491.98 708.28 

1583. 

0 1348 

360 0 

0 .00 

0.71 

1490.10 709.14 

1584. 

TG = 

1583 

1614 


1665. 1490 

. 

0 A - 

1584 

1583. 

1614. 

1671. 1663. 

1*11. 


T2 PE PC PRIN PROUT FRIN FROUT 

689. 2543.65 2543.65 2543.65 2543.65-0.05597-0.23882 
686. 2543.24 2543.24 2543.24 2543.24-0.05608-0.23857 
709. 692. 696. 686. 

787. 708. 691. 688. 696. 686. 


.509 


.506 


.371 



14 PAGES OF LISTING OMITTED 



683. 


0 1967 360.0 0.00 0.71 

TG = 1586. 1617. 

TOA = 1587. 1586. 1617. 

TM= 1663. 1672. 1672. 


1686.88 708.68 lrc 

1668. 1687 

1676. 1666. lido. 

1672. 1672. ltf'S. 


2666.02 2658.08 2666.86 2652.88-0.05327 
709. 691. 693. 683. 

786. 708. 690. 687. 693. 683 

785. 666. 662. 662. 666. 592 

0.617 REFF1= 


-0.23035 


QIH= 386.193 Q0UT= 258.590 WRKEXP= 363.166 WRKCMP --2 02 . 837 WRKT0T= 160.307 EFFT0T= 

QREGEN= 1.387 PWRHP= 19.666 FREQ= 66.780 WRKL1 * ^6.836 UT0TAL= 0.000 

EN3PM= 6310.710 AVGWSP= 2170.861 

QEXPN= -6.357 QEXPP= 0 132 QHEATN=-373 .808 QHt. H = 2.712 QC00LN= -6.873 QC00LP= 200.910 QCOMPN- 

QRr.GN=-1776.058 QREGP = 1775 . 665 QCNDRI= 5.178 QCr.uNJ- 5.178 , 

QCNDCL= 1.388 QCNDD= 0.368 QCNDCN= 0.000 TO t~. 1 5 38 . 628 TGCNPA= 655.607 QSHTL= 1-939 

QINB = 375.320 WRKBAS = 185.161 WALT = 182.160 QC'jIl- 198.893 WRKLEX= 19.586 WRKLCM= 5.250 WLEXR = 


0.997 REFF2= 0.996 


0.000 QC0MPP= 


WLCMR= 2.606 MLEXE= 8.988 WLCME= 1.272 WRKLR= 6.610 WRKLH= 6.791 WRKLC : 


3.806 

375 WRKIE= 10.259 WLALT= 21.833 


* CYCLE NO. 16 X 


TIME 
0.1967 
0 2096 
TG = 

TGA= 1587 
TM= 1663 


ANGLE 
360.0 
360 . 0 
1586 . 


X(l) X(2) 
0.00 0.71 

0.00 0.71 

1617., 
1587. 1617. 

1672. 1672. 


T2 

683. 

633. 

786 
736 


PE PC PRIN PROUT FRIN FROUT 

2666.02 2658.08 2666.86 2652.88-0.05327-0.23035 
2669.86 2669.86 2669.86 2669.86-0.05327-0.23035 
709. 691. 693. 683. 

707. 690. 687. 693. 683. 

666. 662. 662. 666. 592. 

0.617 REFF1= 0 


TRIN TROUT T1 

1686.88 708.68 1587 

1686.99 708.67 1587 

1668. 1687. 

1676. 1666. 1608 

1672. 1672. 1608 

Q 1 N = 386.062 QOUT= 259.159 WRKEXP= 382.789 WRKCMP = -197 . 619 WRKTOT= 160.335 EFFTOT = 

QREGEN= 0.993 PWRHP= 19.668 FREQ= 66.780 WRKLT= 26.836 UTOTAL= 0.000 
EN3PM= 6310.710 AVGWSP= 2171.016 „ 

Q EXPN = -6.355 QEXPP= 0.133 QHEATN=-373 .677 QHEATP= 2.730 QCOOLN= -6.826 QCOOLP= 201.623 QCOMPN- 

Q D EGH = -1776 . 633 QREGP=1775 . 625 QCNDRI= 5.178 QCNDRO= 5.178 

QCNDCL = 1.388 QCNDD= 0.368 QCNDCN= 0.000 TGEXPA = 1538 . 6 91 TGCMPA= 655.675 QSHTL= 1.939 

QJND = 375.169 WRKBAS = 185.170 WALT = 185.170 Q0UTB= 199 . 65S_WRKLEX= 0.000 WRKLCM= 0.008 WLEXR= 

WLCMR= 0.003 WLEXE= 0.000 WLCME= 0.003 WRKLR= 


997 REFF2= 0.996 


0.000 QCOMPP: 


6.610 WRKLH= 6.791 WRKLC= 


0.000 

1.375 WRKLE= 10.259 WLALT= 


0.000 


*» CYCLE NO. 15 x 


TIME 

ANGLE 

xm 

X( 2 ) 

TRIN 

TROUT 

T 1 

0 2096 

360 . 0 

0.00 

0.71 

1686 . 99 

708.67 

1587 . 

0.2266 

360.0 

0.00 

0.71 

1687 . 08 

708.70 

1587. 

TG = 

1586 

1617 


1668. 1687 

. 

TGA = 

1587 . 

1587. 

1617. 

1676. 

1666 . 

1608. 

7 n= 

1663. 

1672. 

1672. 

1672. 

1672. 

1608. 


T2 PE PC PRIN 

683. 2669.86 2669.86 2669.86 

683. 2650.00 2650.00 2650.00 

709. 691. 

786. 708. 690. 687 

786. 666. 662. 662 


PROUT FRIN FROUT 
2669.86-0.05327-0.23035 
2650.00-0.05327-0.23036 
693. 683. 

693. 683. 

666. 592. 

0.618 REFF1= 0.997 REFF2= 


66.780 WRKLT = 26.836 UTOTAL = 


0.000 


QIN= 383.931 QOUT= 259.362 WRKEXP= 

QRFGEN= 0.728 PWRHP= 19.670 FREQ 

QEXPN= -6 1 356 1 QEXPP= SP 0 2 133 QHEATN=-373 . 575 QHEATP= 2.739 QCOOLN= -6.811 QCOOLP= 201.586 QCOMPN= 
QREGN =-1776. 656 QREGP=1775 . 386 QCNDRI= 5.178 QCNDRO= 5.178 


0.996 


0.000 QCOMPP= 


QCNDCL= 1 388 QCNDD= 0.368 QCNDCN= 0.000 TGEXPA=1 538 . 527 TGCMPA= 655.506 QSHTL= 
QINB = 375.058 WRKBAS = 185.193 WALT = 185.193 QOUTB= 199.635 WRKLEX= 0.000 WRKLCM= 
WLCMR= -0.000 WLEXE= 0.000 WLCME= -0.000 WRKLR= 6.610 WRKLH= 6.791 WRKLC= 


1.939 

0.000 WLEXR= 0.000 
1.375 WRKLE= 10.259 WLALT= 


0.000 


AVEPAGE VALUES OVER LAST 5 CYCLES xx 


QIN= 386.380 QOUT= 258.507 WRKEXP= 378.788 WRKCMP =-l 98 . 6 26 WRKTOT= 160.277 EFFTOT= 0.617 REFF1= 0.997 REFF2= 0.996 


0.000 


QREGEN= 0.561 PWRHP= 19.661 FREQ= 66.780 WRKLT= 26.851 UTOTAL- 
EN3PM= 6310.710 AVGWSP= 2170.770 

QEXPN= -6.359 QEXPP= 0.131 QHEATN=-373 . 933 QHEATP= 2.656 QCOOLN= 

QRCGN=-1776 776 QREGP=1775 . 335 QCNDRI= 5.178 QCNDRO= 5.178 

QCNDCL = 1.388 QCNDD= 0.368 QCNDCN= 0 . 000 TGEXPA = 1538 . 285 TGCMPA= 655.329 QSHTL- 


6.875 QCOOLP= 200.866 QCOMPN= 


1.939 


0.000 QCOMPP= 


.855 


.861 


.862 


.852 



QINB = 375.507 WRKBAS = 185.129 WALT = 184.529 QOUTB= 198.821 WRKLEX= 3.917 WRKLCM= 1.053 WLEXR= 0.761 

WLCMR= 0.522 WLEXE= 1.798 WLCME= 0.256 WRKLR= 6.410 WRKLH= 6.796 WRKLC= 1.376 WRKLE= 10.270 WLALT= 4.367 


AVG TEMPS. --LAST 5 CYCLES, AVG. H. T. COEFS. S HEAT FLUXES— LAST CYCLE X* 


TGCYC= 1530. 1528. 1524. 1482. 702. 690. 669. 659. 


TGACYC= 

1539 

1527. 1545. 

1542. 1505. 

1400. 

780. 

696. 679. 

668. 

665. 

656 . 

TMCYC= 

1643. 

1672. 1672. 

1672. 1672. 

1400 . 

780 . 

646. 642. 

642. 

646 . 

592. 

HACYC= 

0.057 

0.000 0.808 

0. 

810 0.000 

1.922 

1.879 

0.000 0.622 

0.671 

0.000 

0.075 

HMX = 

0 444 

0.000 1 412 

1.435 0.000 

2.949 

2.779 

0.000 1.027 

1.099 

0.000 

0.245 

Q 0 A A V G = 

5.66 

0.00 69.19 

72 

.96 0.00 

4.45 

4.38 

0.00 22.62 

19.19 

0.00 

5.01 

QOAM\= 

22 99 

0.00 132.72 

128 

.79 0.00 

8.45 

8.52 

0.00 46.24 

59.02 

0.00 

25.80 

* CYCLE 

NO. 16 

X 









TIME 

ANGLE 

X(l) X( 2 ) TRIN 


TROUT 

T1 T2 

PE PC PRIN PROUT FRIN FROUT 

0 2246 

360.0 

0.00 0.71 1487 

.08 

708.70 

1587. 683. 2450 

.00 2450.00 2450 

.00 2450 

.00-0.05327-0.23034 

0 2396 

360.0 

0.00 0.71 1487 

. 16 

708.72 

1587. 683. 2446 

.43 2458.49 2447 

.26 2453 

.28-0.05327-0.23034 

TG = 

1586 

1617. 

1668. 1487 


709. 

691. 

693. 

683. 


T C A = 

1587 . 

1587. 1617. 1674. 

1666 . 

1408. 

784. 

708. 690. 

687. 

693. 

683. 

TM = 

1643. 

1672. 1672. 1672. 

1672. 

1408. 

786. 

646. 642. 

642. 

646 . 

592. 

QIN= 383.845 Q0UT= 259.495 WRKEXP 

= 363.255 WRKCMP = -2 02 .883 WRKTOT= 160.368 EFFTOT= 0.418 REFF1 

= 0.997 REFF2 

QREGEN= 

0.598 

PWRHP = 19.472 FREQ = 

66.780 WRKLT = 24. 

837 UTOTAL= 0.000 





EN3PM= 6310.710 AVGWSP= 2171.243 

QEXPN= -4 354 QEXPP= 0.133 QHEATN=-37 3.496 QHEATP= 2.744 QCOOLN= -4.799 QCOOLP= 201.722 QCOMPN= 0.000 QCOMPP= 2.863 

QRFGN =-1774. 759 QREGP=1775 357 QCNDRI= 5.178 QCNDRO= 5.178 

QCNDCL= 1.388 QCNDD= 0.368 QCNDCN= 0 000 TGEXPA=153S . 553 TGCMPA= 655.527 QSHTL= 1.939 

►-.QINB = 374.972 WRKBAS = 185.205 WALT = 182.203 QOUTB= 199.786 WRKLEX= 19.587 WRKLCM= 5.251 WLEXR= 3.807 

'— * WLCilR= 2.604 WLEXE= 8.988 WLCME= 1.272 WRKLR= 6.412 WRKLH= 6.791 WRKLC= 1.375 WRKLE= 10.260 WLALT= 21.836 

* CYCLE NO. 17 x 

TIME ANGLE X(l) X(2) TRIN TROUT T1 T2 PE PC PRIN PROUT FRIN FROUT 

0 2396 360.0 0.00 0.71 1487.16 708.72 15S7. 683. 2446.43 2458.49 2447.26 2453.28-0.05327-0.23034 

0 2546 360.0 0.00 0.71 1487.22 708.81 1537. 683. 2450.29 2450.29 2450.29 2450.29-0.05327-0.23034 

TG= 1586. 1617. 1668. 1487. 709. 691. 693. 683. 

TGA= 1587. 1587. 1617. 1674. 1666. 14U9. 784. 708. 690. 687. 693. 683. 

TN= 1643. 1672. 1672. 1672. 1672. 1408. 786. 646. 642. 642. 646. 592. 

QIN = 383 768 QOUT= 259.249 WRKEXP= 382.843 WRKCMP=-197 . 659 WRKTOT= 160.347 EFFTOT = 0.418 REFF1= 0.997 REFF2= 0.996 

QREGEN= 0.554 PWRHP= 19.469 FREQ= 66.780 l!RKLT= 24.837 UTOTAL= 0.000 
EN3FM= 6310.710 AVG'1SP= 2171.364 

Q EXPN = -4.353 QEXPP= 0.134 QHEATN = -373 . 424 QHEATP= 2.749 QCOOLN= -4.822 QCOOLP= 201.491 QCOMPN= 0.000 QCOMPP= 2.863 

QRLGN=-1774 . 784 Q R EGP = 1 7 7 5 . 3 38 QCNDRI= 5.178 QC‘JURO= 5.178 

QCNDCL = 1.388 QCNDD= 0 368 QCNDCN= 0.000 TGEYPA=1538 . 579 TGCMPA= 655.563 QSHTL= 1.939 

QINB = 374 895 WRKBAS = 185.184 WALT = 185.184 QOU.b= 199.532 WRKLEX= 0.000 WRKLCN= 0.008 WLEXR= 0.000 

WICMR= 0.003 WLEXE= 0.000 WLCNE= 0.003 WRKLR= 6.412 WRKLH= 6.791 WRKLC= 1.375 WRKLE= 10.260 WLALT= 0.000 

_ < CYCLE NO. 18 « 

" TIME ANGLE X(l) X(2) TRIN TROUT T1 

0 2546 360.0 0 00 0.71 1487.22 708.81 1587. 

0.2695 360.0 0.00 0 71 1487.27 708.82 1587. 

TG= 1586. 1617. 1668. 1487. 

T GA= 1587. 1587 . 1617. 1674. 1666 . 1409. 

TM= 1643. 1672. 1672. 1672. 1672. 1409. 


T2 PE PC PRIN PROUT FRIN FROUT 

683. 2450.29 2450.29 2450.29 2450.29-0.05327-0.23034 
683. 2450.38 2450.38 2450.38 2450.38-0.05327-0.23034 
709. 691. 693. 683. 

784. 708. 690. 687. 693. 683. 

786. 646. 642. 642. 646. 592. 



1 


QIN= 383.717 QOUT = 259.456 WRKEXP= 382.868 WRKCMP=-1 97 . 6 6 7 WRKTOT = 160.363 EFFTOT = 

QRFCEN- 0.554 PWRHP= 19.471 FREQ= 66.780 WRKLT= 24.837 UTOTAl= 0.000 

EN3PM= 6310.710 AVGWSP= 2171.462 , „ , „„„„„„ 

QFXf N = -4.353 QEXPP= 0.134 QHEATN=-373 . 378 QHEATP= 2.753 QCDOLN= -4.808 QCOOLP= 201.681 QCOMPN- 

0CTCN=-1774 769 QREGP = 1775 . 323 QCHDRI = 5.178 QCNDRO= 5 178 

QCNDCL = 1.388 QCNDD= 0.368 QCNDCN= 0 000 TGF.XP A = 1 5 38 . 5 9 5 TGCMPA= 655.604 QSHTL= 1.939 

q j MB = 374.844 WRKBAS = 185.200 UALT = 185.200 Q OUT B = 199.739 WRKL EX= 0.000 WRKLCM= -0.000 WLEXR- 
UICI1R = -0 000 WLEXE= 0.000 WLCME= -0.000 WRKLR= 


0.418 REFF1= 0.997 REFF2= 0.996 


0.000 QCOMPP= 


6.412 WRKLH= 6.791 WRKLC= 1.375 WRKLE= 


0.000 

10.260 WLALT= 


0.000 


" CYCLE 
T I ME 
0 26 95 
0 284 5 
TG = 

TGA = 


NO 19 
ANGLE 
360 0 
360.0 
1586 
1587. 


XC 1 ) X(2) TRIN TROUT T1 

0.00 0.71 1487.27 708.82 1587. 

0.00 0.71 1487.31 703.84 1587. 

1617. 1668. 1487. 

1587. 1617. 1674. 1666. 1409. 

1672. 1672. 1672. 1672. 1409. 


T2 

683. 

683. 


PE 

2450.38 
2446.78 
709. 

785. 708 


786 . 


646. 


PC 

2450.38 

2458.84 

691. 

690 . 
642. 


.590 WRKEXP= 363.296 WRKCMP =-20 2 . 92 9 WRKTOT= 160.367 EFFTOT= 


PP.IN PROUT FRIN 
2450.38 2450.38-0.05327-0 
2447.61 2453.64-0.05327-0 
693. 683. 

687. 693. 683. 

642. 646. 592. 

0.418 REFF1= 0 


FROUT 

,23034 

,23034 


997 REFF2= 0.996 


lfl= 1643 

QIN = 383 677 QOUT = 259. 

0 REGEN- 0.443 PWRHP= 19.472 FREQ= 

QEXPN = -^SS^QEXPP^ 515 0 f 1 34 ' QHEATN=-373 . 34 1 QHEATP= 2.756 QCOOLN= -4.797 QCOOLP= 201.802 QCOMPN= 0.000 QCOMPP= 2 
QPEGN=-1774 .848 QREGP=1775 291 QCHDRI = 5.178 QCNDRO= 5.178 

QCNDCL = 1.388 QCNDD= 0.368 QCNDCH= 0.000 TGEXPA = 1538 . 6 06 TGCNPA= 

Q J MB = 374.804 WRKBAS = 185.206 WALT = 182.205 QOUTB= 199.872 WRKLEX = 

WLCMR= 2.605 WLEXE= 8.989 WLCME= 1.272 WRKLR= 


66.780 WRKLT = 24.839 UTOTAL = 


6.413 WRKLH : 


0.000 


655.631 QSHTL= 
19.588 WRKLCM= 
.791 WRKLC= 1 


1.939 

5.251 WLEXR= 3. 
375 WRKLE= 10.261 


808 

WLALT= 


21.837 


“ CYCLE NO. 20 * 
TIME ANGLE 
0 2845 360.0 

0 2995 360.0 

TG- 1586 . 

r\j TGA= 1587. 

TM= 1643. 


XC 1 ) 

0.00 

0.00 0.71 

1617. 

1587. 1617 

1672 1672 


X( 2 ) TRIN 
0.71 1487.31 

1487.34 


TROUT T1 

708.84 1587. 

708.90 1587. 

1668. 1487. 

1674. 1666. 1409. 

1672. 1672. 1409. 


T2 

683. 

683. 


PE 


PC 


PRIN 


2446.78 2458.84 2447.61 


PROUT FRIN 
2453.64-0.05327 


FROUT 

-0.23034 


2446.90 2458.96 2447.73 2453.75-0.05327-0.23033 


709 . 

785. 708 

786. 646 


691. 


690. 
642 . 


687 

642 


QIN= 383 640 QOUT = 259.469 WRKEXP= 363.297 WRKCMP = -202 . 948 WRKTOT= 160.350 EFFTOT = 

0.000 


694. 683. 

693. 683. 

646. 592. 

0.418 REFF1= 0. 


997 REFF2= 


0.996 


QRFGEN= 0.404 PWRHP= 19.469 FREQ= 66.780 WP.KLT= 24.848 UTOTAL = 

EN3PM= 6310.710 AVGWSP= 2171.620 

QEXPN- -4.352 QEXPP= 0 134 QHEATN=-373 . 307 QHEATP= 2.758 QCOOLN= 

QREGN =-1774. 855 QREGP=1775 260 QCHDRI = 5.178 QCHDRO= 5.178 

QCNDCL= 1.383 QCNDD= 0.368 QCNDCN= 0.000 TGEXP A =1 5 38 . 6 2 0 TGCMPA= 655^661 QSHTL- 
QINB = 374 767 WRKBAS = 185.197 WALT = 182.188 QOUTB= 199.746 WRKL EX= 

WLCMR= 2 608 WLEXE= 8.989 WLCME= 1.275 WRKLR= 6.416 WRKLH= 6 


-4.809 QCDOLP= 201.687 QCOMPN= 0.000 QCOMPP= 


1.939 

19.588 WRKLCM= 5.259 WLEXR= 
791 WRKLC= 1.377 WRKLE= 10. 


3.808 

264 WLALT= 21.838 


-rx AVERAGE VALUES OVER LAST 5 CYCLES ** 

Q1N= 383.729 QOUT= 259.452 WRKEXP= 371.112 WRKCMP=-200 . 818 WRKT0T= 160.359 EFFTOT= 0.418 REFF1= 0.997 REFF2= 0.996 


QPrGEN= 0.511 PWRHP = 19.470 FREQ= 66.780 WRKLT= 24.840 UTOTAL = 

EN3PM= 6310 710 AVGUSP= 2171.445 

QEXPN- -4.353 QEXPP= 0.134 QHEATN=-373 . 389 QHEATP= 2.752 QCOOLN= 
QREGN=-1774 803 QREGP=1775 . 314 QCNDRI= 5.178 QCNDRO= 5.178 
QCNDCL - 1.383 QCNDD= 0.368 QCNDCN= 0.000 TGEXPA=1538 . 591 TGCMPA= 

QINB = 374 856 WRKBAS = 185.199 WALT = 183.396 Q0UTB= 199.735 WRKL EX = 
WLCMR= 1.564 WLEXE= 5.393 WLCME= 0.764 WRKLR= 6.413 WRKLH= 6. 


0.000 

-4.807 QCOOLP= 201.677 QCOMPN= 0.000 QCOMPP= 


655.597 QSHTL= 
11.753 WRKLCM= 


1.939 
3.154 WLEXR= 


2.285 


791 WRKLC= 1.375 WRKLE= 10.261 WLALT= 13.102 


.866 


.867 


.868 


.865 



xx AVG. TEMPS. — LAST 5 CYCLES, AVG. H. T. COEFS. & HEAT FLUXES— LAST CYCLE XX 


TCCYC= 1530. 1528. 1524. 1482. 703. 691. 669. 659. 


TGACYC= 

1539. 

1527. 

1545. 

1543. 

1505. 

1400 . 

780. 

697. 

679. 

668. 

665. 

656. 

TMCYC= 

1643. 

1672. 

1672. 

1672. 

1672. 

1400 . 

780. 

646 . 

642. 

642. 

646 . 

592. 

HACYC= 

0.057 

0.000 

0.808 

0.810 

0.000 

1.922 

1.879 

0.000 

0.622 

0.672 

0.000 

0.075 

h'MX = 

0.444 

0.000 

1.412 

1.435 

0.000 

2.949 

2.779 

0.000 

1.027 

1.099 

0.000 

0.245 

QO A AVG= 

5.65 

0.00 

69.15 

72.90 

0.00 

4.45 

4.38 

0.00 

22.64 

19.20 

0.00 

5.02 

00 AMX= 

22.92 

0.00 

132.68 

128.64 

0.00 

8.45 

8.52 

0.00 

46.23 

59.02 

0.00 

25.83 


x CYCLE HO 21 x 


TIME 

AHGLE 

X(l) 

X( 2 ) 

TRIH 

TROUT 

T1 

T2 

PE 

PC 

PRIH 

PROUT FRIH FROUT 

0 2995 

360.0 

0.00 

0.71 

1487.34 

708.90 

1587. 

683. 

2446.90 

2458.96 

2447.73 

2453.75-0.05327-0.23033 

0 3001 

14.4 

0.03 

0.52 

1485.75 

700.03 

1614. 

695. 

2587.92 

2607.44 

2593.97 

2601.99-0.11060-0.26132 

0 . 3007 

28.8 

0.12 

0.35 

1483.80 

694.10 

1645. 

703. 

2691.60 

2720 .80 

2706.57 

2715.93-0.16240-0.26890 

0 3013 

43.2 

0.25 

0.21 

1481.42 

689.33 

1667. 

707 . 

2748.13 

2785.85 

2772.35 

2782.21-0.19812-0.24991 

0.3019 

57.6 

0.42 

0.10 

1478.76 

684.29 

1673. 

708. 

2756.66 

2798.16 

2786.50 

2795.99-0.21164-0.20731 

0 3025 

72.0 

0.62 

0.03 

1476.06 

678.27 

1663. 

705 . 

2722.45 

2761.48 

2752.22 

2760.56-0.20240-0.14794 

0.3031 

86.4 

0.82 

0.00 

1473.57 

671.10 

1643. 

699 . 

2653.89 

2685.04 

2678.34 

2684.87-0.17352-0.07952 

0 3037 

100.8 

1.01 

0.01 

1471.45 

662.95 

1618 . 

691 . 

2559.59 

2580.05 

2575.85 

2580.05-0.13060-0.00890 

0 .3043 

115.2 

1 19 

0.06 

1469.78 

692.98 

1592. 

683. 

2446.70 

2456.33 

2455.00 

2456.65-0.08194 0.06085 

0.3049 

129.6 

1.34 

0.15 

1468.14 

695.38 

1566. 

672. 

2324.08 

2324.60 

2326.67 

2325.61-0.03289 0.11694 

0 3055 

144.0 

1.45 

0.27 

1410.40 

696.81 

1541. 

660 . 

2196 . 98 

2191.65 

2197.13 

2193.47 0.01133 0.16000 

0.3061 

158.4 

1.53 

0.43 

1391.33 

698.39 

1517. 

647 . 

2072.12 

2062.96 

2071.59 

2065.55 0.05149 0.19115 

0.3067 

172 8 

1.57 

0.61 

1381.49 

700.19 

1495. 

635. 

1959.38 

1945.20 

1956.55 

1948.40 0.08421 0.21058 

0.3073 

187.2 

1.57 

0.81 

1384.14 

702.11 

1474. 

624. 

1861.41 

1842.01 

1855.23 

1845.49 0.10744 0.21681 

0 3079 

201.6 

1 53 

1.01 

1399.00 

704.04 

1457 . 

616 . 

1780.48 

1756.53 

1770.65 

1759.94 0.12198 0.21073 

0 3085 

216 . 0 

1.45 

1.20 

1422.41 

705.88 

1444. 

610. 

1718.92 

1691.59 

1705.65 

1694.58 0.12957 0.19441 

0 3091 

230.4 

1.34 

1.36 

1449.66 

707.58 

1435. 

607 . 

1679.10 

1649.65 

1662.86 

1652.00 0.13231 0.17031 

0 3097 

244.8 

1 . 19 

1.49 

1476.87 

709.08 

1431 . 

606 . 

1663.25 

1632.88 

1644.67 

1634.49 0.13185 0.14047 

0.3103 

259.2 

1.01 

1.56 

1502.24 

710.38 

1434. 

609 . 

1673.33 

1643.26 

1653.24 

1644.16 0.12896 0.10597 

0.3109 

273.6 

0.82 

1.57 

1525.99 

711.47 

1443. 

613. 

1711.01 

1682.57 

1690.50 

1682.90 0.12331 0.06687 

0.3115 

288.0 

0.62 

1.53 

1549.57 

712.50 

1458. 

621. 

1777.68 

1752.57 

1758.30 

1752.55 0.11346 0.02235 

0.3121 

302.4 

0.42 

1.43 

1574.24 

730.44 

1478. 

631 . 

1873.13 

1853.32 

1856.56 

1852.83 0.09899-0.02473 

0 3127 

316.8 

0.25 

1.28 

1600.17 

735.52 

1503. 

643. 

1992.59 

1980.71 

1980.69 

1979.20 0.07674-0.07476 

0.3133 

331.2 

0.12 

1.11 

1627.72 

731.75 

1530 . 

656 . 

2132.71 

2130.48 

2126.56 

2127.59 0.04431-0.12968 

0 3139 

345.6 

0.03 

0 . 91 

1656.73 

720.86 

1559. 

670 . 

2287.77 

2294.28 

2286.44 

2290.04 0.00136-0.18351 

0.3145 

360.0 

0.00 

0.71 

1487.36 

708.91 

1587. 

683. 

2446.97 

2459.03 

2447.80 

2453.82-0.05327-0.23033 

TG = 

1586 

1617. 

1668. 1487. 


709. 

691 . 


694. 683. 

T GA = 

1587. 

1587. 

1617 

1674. 

1666 . 

1409. 

785 

708 

690 

687 

693. 683. 

TM = 

1643. 

1672. 

1672 

1672. 

1672. 

1409. 

786 

646 

642 

642 

646. 592. 

QIH = 383.617 QOUT= 259. 

622 1 

WRKEXP= 363.312 WRKCMP=-202 . 959 

WRKTOT= ; 

160.353 

EFFTOT= 

0.418 REFF1= 0.997 REFF2 

QREGEH= 

0.398 

PWRHP= 

19. 

470 FREQ = 

66.780 

URKLT= 

24.848 

UTOTAL = 

0.000 



EH3PM= 6310.710 AVGUSP= 2171.688 

QEXPM= -4.352 QEXPP= 0.134 QHEATH=-373 . 286 QHEATP= 2.760 QCOOLN= -4.798 QCOOLP= 201.827 QCOMPN= 0.000 QCOMPP= 2.870 
QREGH=-1774.834 QREGP =1775 . 232 QCHDRI= 5.178 QCNDRO= 5.178 

QCMDCL= 1.388 QCNDD= 0.368 QCNDCN= 0.000 TGEXPA=1538 . 627 TGCNPA= 655.691 QSHTL= 1.939 
” QIHB = 374.744 WRKBAS = 185.201 WALT = 182.191 QOUTB= 199.899 WRKL EX = 19.589 WRKLCM= 5.259 WLEXR= 3.808 

ULCMR= 2.608 WLEXE= 8.989 WLCME= 1.275 WRKLR= 6.416 WRKLH= 6.791 WRKLC= 1.377 WRKLE= 10.264 WLALT= 21.838 

x CYCLE HO. 22 X 

TIME AHGLE X(l) XC2) TRIH TROUT T1 T2 PE PC PRIN PROUT FRIH FROUT 

0.3145 360.0 0.00 0.71 1487.36 708.91 1587. 683. 2446.97 2459.03 2447.80 2453.82-0.05327-0.23033 



0 3151 

14.4 

0.03 

0.52 

1485.75 

700.04 

1614. 

695 . 

2587.96 

2607.48 

2594.00 

0 3157 

28.8 

0 . 12 

0 35 

1483 81 

694.11 

1645. 

703 . 

2691.64 

2720.84 

2706.60 

0 3163 

43.2 

0.25 

0.21 

1481 42 

689.34 

1667 . 

707 . 

2748.17 

2785.89 

2772.38 

0.3169 

57.6 

0.42 

0.10 

1478.76 

684 . 30 

1673. 

70S. 

2756.69 

2798 .20 

2786.53 

0.3175 

72.0 

0.62 

0.03 

1476.07 

678.28 

1663. 

705 . 

2722.48 

2761.51 

2752.25 

0 3181 

86.4 

0.82 

0 .00 

1473.58 

671.11 

16 4 3 . 

699. 

2653 . 91 

2685.07 

2678.57 

0.3187 

100.8 

1.01 

0.01 

1471.46 

662.95 

1618. 

691. 

2559.61 

2580.08 

2575.87 

0.3193 

115.2 

1.19 

0.06 

1469 79 

692.99 

1592. 

683. 

2446.72 

2456.35 

2455.02 

0 3199 

129.6 

1.34 

0.15 

1468.14 

695.39 

1566. 

672. 

2324.10 

2324.62 

2326.69 

0 3205 

144.0 

1.45 

0.27 

1410.40 

696.82 

1541. 

660 . 

2197.00 

2191.67 

2197.15 

0 3211 

158.4 

1.53 

0.43 

1391.34 

698.40 

1517. 

647 . 

2072.14 

2062.98 

2071.60 

0 3217 

172.8 

1.57 

0.61 

1381.49 

700.20 

1495. 

635. 

1959.40 

1945.22 

1956.57 

0.3223 

187 .2 

1.57 

0.81 

1384.14 

702.12 

1474. 

624. 

1861.43 

1842.03 

1855.25 

0 3229 

201.6 

1.53 

1.01 

1399.00 

704.04 

1457 . 

616 . 

1780.49 

1756.55 

1770.67 

0 3235 

216.0 

1.45 

1 20 

1422.41 

705.89 

1444. 

610. 

1718.94 

1691.61 

1705.67 

0.3240 

230.4 

1.34 

1.36 

1449.66 

707.58 

1435. 

607 . 

1679.12 

1649.67 

1662.88 

0 3246 

244 8 

1.19 

1.49 

1476 .87 

709.09 

1431. 

606 . 

1663.26 

1632.90 

1644.69 

0 3252 

259.2 

1.01 

1.56 

1502.24 

710.38 

1434. 

609. 

1673.35 

1643.27 

1653.25 

0 3258 

273.6 

0.82 

1.57 

1525.99 

711.48 

1443 . 

613. 

1711.03 

1682.58 

1690.52 

0.3264 

288.0 

0.62 

1.53 

1549.57 

712.51 

1458 . 

621. 

1777.69 

1752.58 

1758.32 

0.3270 

302.4 

0.42 

1.43 

1574 24 

730.45 

1478 . 

631 . 

1873.15 

1853.34 

1856.58 

0 . 3276 

316 .8 

0.25 

1.28 

1600 . 17 

735.53 

1503. 

643. 

1992.61 

1980.73 

1980.71 

0 3282 

331.2 

0.12 

1.11 

1627.72 

731.76 

1530 . 

656 . 

2132.73 

2130.50 

2126.58 

0 . 3288 

345.6 

0.03 

0.91 

1656.73 

720.86 

1559. 

670. 

2287.79 

2294.30 

2286.46 

0.3294 

360.0 

0.00 

0.71 

1487.36 

708.92 

1587. 

683. 

2446.99 

2459.05 

2447.82 

TG = 

1586 

1617 


1668. 1487. 


709. 

691. 


TGA = 

1587. 

1587. 

1617. 

1674. 

1666 . 

1409. 

785 

708. 

690. 

687 

TM = 

1643. 

1672. 

1672. 

1672. 

1672. 

1409. 

786 

646. 

642. 

642 


2602 02-0.11060-0 
2715.97-0.16240-0 
2782.25-0.19812-0 
2796 02-0.21164-0 
2760.59-0.20240-0 
2684.90-0.17352-0 
2580.07-0.13060-0 
2456.67-0.08194 0 
2325.63-0.03289 
2193.49 0.01133 
2065 57 0.05149 
1948.42 0.08421 

1845.51 0.10744 
1759.95 0.12198 
1694.60 0.12957 
1652.02 0.13231 

1634.51 0.13186 
1644.17 0.12896 


1682.91 0.12331 


QIN= 383.616 QOUT= 259.672 WRKEXP= 363.317 WRKCMP=-202 . 964 WRKT0T= 160.353 EFFTOT= 


1752.57 0.11346 
1852.85 0.09899-0 
1979.22 0.07674-0 
2127.61 0.04431-0 
2290.06 0.00136-0 
2453.85-0.05327-0 
694. 683. 

693. 683. 

646. 592. 

0.418 P.EFF1= 0 


.26132 
.26890 
.24990 
.20730 
.14794 
.07952 
. 00890 
.06085 
.11694 
. 16000 
.19115 
21058 
.21681 
.21073 
.19441 
. 17031 
.14047 
.10597 
.06687 
.02235 
.02473 
.07476 
.12968 
. 18351 
.23033 


997 REFF2= 0.996 


0.000 


QRCGEH= 0.517 PWRHP= 19.470 FREQ= 66.780 WRKLT= 24 848 UTOTAL= 

QEXPN= -4 1 352 1 QEXPr= SP Q 2 134 QHEATN=-373 . 285 QHEATP= 2.760 QCOOLN= -4.793 QCOOLP= 201.875 QCOMPN= 0.000 QCOMPP= 2 
QREGM=-1 7 7 4.713 QREGP = 177 5 . 23 0 QCHDRI= 5.178 QCf(DR0= 5.178 
QCMDCL = 1.388 QCNDD= 0 368 QCNDCN= 0.000 TGEXPA = 1538 . 626 

OIHB = 374.743 WRKBAS = 185.202 WALT = 182.192 Q0UTB= 

WI.CHP.= 2.608 WLEXE= 8.989 WLCME= 1.275 WRKLR = 


_ TGCMPA= 
199.952 WRKL EX= 
f. <11 6 URKI H= 6 


655.708 QSHTL 
19.589 WRKLCM 
7 01 i.ipvi r = 


1.939 

5.260 WLEXR= 


7 777 I.IPIM Tr 


3.808 

264 WLALT= 21.838 


CYCLE 

TIME 

NO. 23 # 
ANGLE 

X(l) 

X(2) 

TRIN 

TROUT 

T1 

T2 

PE 

PC 

PRIN 

0 . 3294 

360.0 

0.00 

0.71 

1487.36 

708.92 

1587 . 

683. 

2446 . 99 

2459.05 

2447.82 

0 3300 

14.4 

0.03 

0.52 

1485.75 

700.04 

1614. 

695. 

2587 . 98 

2607.50 

2594 . 02 

0 3306 

28.8 

0.12 

0.35 

1483.81 

694.11 

1645 . 

703 . 

2691.66 

2720 86 

2706.62 

0 3312 

43.2 

0.25 

0.21 

1481.42 

689.35 

1667 . 

708. 

2748.19 

2785.91 

2772.41 

0 3318 

57.6 

0.42 

0.10 

1478.76 

684.30 

1673. 

708. 

2756.71 

2798.22 

2786.55 

0 .3324 

72.0 

0.62 

0.03 

1476.07 

678.28 

1663. 

705. 

2722.50 

2761.53 

2752.27 

0 3330 

86.4 

0.82 

0.00 

1473.58 

671.12 

1643. 

699. 

2653.93 

2685.09 

2678.39 

0 3336 

100.8 

1.01 

0.01 

1471.46 

662.96 

1618. 

691 . 

2559.63 

2580.09 

2575.89 

0.3342 

115 2 

1.19 

0.06 

1469.79 

692.99 

1592. 

683. 

2446.74 

2456.37 

2455.04 

0.3348 

129.6 

1.34 

0.15 

1468.15 

695.40 

1566 . 

672. 

2324.12 

2324.64 

2326.70 

0 3354 

144.0 

1.45 

0 27 

1410.40 

696.83 

1541 

560 . 

2197.01 

2191.69 

2197.16 

0.3360 

158.4 

1.53 

0.43 

1391.34 

698.41 

1517. 

647. 

2072.15 

2062.99 

2071.62 

0 3366 

172.8 

1.57 

0.61 

1381.49 

700.21 

1495. 

635. 

1959.41 

1945.23 

1956.58 

0 3372 

187.2 

1.57 

0.81 

1384.15 

702.12 

1474. 

624. 

1861.44 

1842.04 

1855.26 

0 3378 

201.6 

1.53 

1.01 

1399.00 

704.05 

1457 . 

616. 

1780.50 

1756.56 

1770.68 


PROUT FRIN FROUT 
2453.85-0 . 05327-0.23033 

2602.04- 0.11060-0.26132 
2715.99-0.16240-0.26890 
2782.27-0.19812-0.24990 

2796.04- 0.21164-0.20730 
2760.61-0.20240-0.14794 
2684.91-0.17352-0.07952 
2580.09-0 .13060-0.00890 
2456.69-0.08194 0.06085 
2325.65-0.03289 0.11694 

2193.51 0.01133 0.16000 
2065.58 0.05149 0.19115 
1948.43 0.08421 0.21058 

1845.52 0.10744 0.21681 
1759.96 0.12198 0.21073 


.870 



1 


0 3386 

216.0 

1.65 

1.20 

1622.62 

705.90 

1666. 

610 . 

1718.95 

1691.62 

1705.68 

1696.61 

0.12957 0.19661 

0 3390 

230 6 

1.36 

1 . 36 

1669.66 

707.59 

1635. 

607 . 

1679.13 

1669.68 

1662.89 

1652.03 

0.13231 0.17031 

0 3396 

266.8 

1.19 

1.69 

1676.87 

709.10 

1631. 

606 . 

1663.27 

1632.91 

1666.70 

1636.52 

0.13186 0.16067 

0 3602 

259.2 

1.01 

1.56 

1502.26 

710.39 

1636. 

609. 

1673.36 

1663.28 

1653.26 

1666.18 

0.12896 0.10597 

0 . 36 08 

273.6 

0.82 

1.57 

1525.99 

711.68 

1663. 

613. 

1711.06 

1682.59 

1690.53 

1682.92 

0.12331 0.06687 

0.3616 

288.0 

0.62 

1.53 

1569.57 

712.51 

1658. 

621. 

1777.70 

1752.59 

1758.33 

1752.58 

0.11366 0.02235 

0 3620 

302.6 

0.62 

1.63 

1576.26 

730.66 

1678. 

631. 

1873.16 

1853.35 

1856.59 

1852.86 

0.09899-0.02672 

0 . 3626 

316.8 

0.25 

1.28 

1600.17 

735.56 

1503. 

663. 

1992.62 

1980.76 

1980.72 

1979.26 

0.07676-0.07676 

0.3632 

331.2 

0 . 12 

1 .11 

1627.72 

731.76 

1530 . 

656 . 

2132.76 

2130.51 

2126.59 

2127.62 

0.06631-0.12968 

0.3638 

365.6 

0.03 

0.91 

1656.73 

720.87 

1559. 

670 . 

2287.80 

2296.32 

2286.67 

2290.08 

0.00136-0.18351 

0 3666 

360.0 

0.00 

0.71 

1687.37 

708.92 

1587. 

683. 

2667.00 

2659.07 

2667.86 

2653.86 

-0.05327-0.23033 

TG = 

1586 

1617. 

1668. 1687. 


709. 

691. 


696. 

683. 

CA = 

1587. 

1587. 

1617 

1676. 

1666. 

1609. 

785 

708 

690 

687. 

693 

683. 

1M = 

1663. 

1672. 

1672 

1672. 

1672. 

1609. 

786 

666 

662 

662. 

666 

592. 

!IN= 383.615 QOUT= 259. 

706 WRKEXP= 363 

.320 WRKCMP=-202.966 

WRKTOT = 

160.356 

EFFTOT= 

0.618 1 

REFF1= 0.997 REFF2 

!REGEN= 

0.508 

PWRHP= 

19.' 

670 FREQ= 

66.780 

WRKLT= 

26.868 

UTOTAL = 

0.000 




cN3PM= 6310.710 AVGWSP= 2171.726 

QEXPN= -6.352 QEXPP = 0.136 QHEATN=-373 .286 QHEATP= 2.760 QCOOLN= -6.791 QCOOLP= 201.906 QCOMPH= 0.000 QCOMPP= 2.871 

QREGN=-1776 .715 QREGP=1775 . 223 QCNDRI= 5.178 QCNPRO= 5.178 

9CNDCL = 1.388 QCNDD= 0.368 QCNDCN= 0.000 TGEXP A = 1538 . 626 TGCMPA= 655.718 QSHTL= 1.939 

9 1 MB = 376.762 WRKBAS = 185.202 WAIT = 182.192 QOUTB= 199.986 WRKLEX= 19.589 WRKLCM= 5.260 WLEXR= 3.808 

Wl CMR = 2.608 WLEXE= 8.989 WICME= 1.275 WRKLR= 6.616 WRKLH= 6.792 WRKLC= 1.377 WRKLE= 10.266 WLALT= 21.838 


x CYCLE NO. 26 * 


TIME 

ANGLE 

REEXP 

RECD1 

REIH 

REOH 

RECD2 

REIR 

REOR 

RECP3 

REIC 

REOC 

RECD6 

RECMP 

0 0000 

360.0 

6367. 

5293. 

8689. 

18023. 

12117. 

63. 

195. 

113300. 

36677. 

61650 . 

265363. 

289173. 

0 0006 

16.6 

38789. 

39066. 

61602. 

68668. 

27650 . 

78. 

228. 

128622. 

60636. 

66361. 

276889. 

296081 . 

0 0012 

28.8 

69966. 

69510 . 

71616. 

79639. 

63093. 

110. 

261. 

131861. 

60759. 

63210. 

260688. 

273661. 

0 0013 

63.2 

93816. 

93206 . 

96035. 

106260 . 

55156. 

130. 

230 . 

122256. 

37179. 

38161. 

223116. 

229201. 

_ 0.0026 

57.6 

108831. 

109025. 

112581. 

122653. 

61106. 

136. 

197. 

101161. 

30260. 

29872. 

168017. 

167510. 

tn 0 0030 

72.0 

116026. 

115395. 

118557. 

126500. 

60076. 

128. 

167 . 

71883. 

20976. 

19561 . 

102625. 

96200. 

0 0036 

86.6 

108288. 

110273 . 

111998. 

112560. 

52731. 

108. 

87. 

38072. 

10693. 

8299. 

33226. 

22633. 

0 0062 

100.8 

92612. 

96163. 

96107 . 

89802. 

60723. 

79. 

23. 

2926. 

696. 

2931. 

36102. 

68160. 

0 .0068 

115.2 

69722. 

70678. 

69133. 

61661. 

26567 . 

66. 

39. 

32631. 

11260. 

13899. 

96826 . 

111939. 

0 0056 

129.6 

63559. 

63655. 

61006. 

31266. 

11866. 

16. 

92. 

60872. 

20331. 

23639. 

152090 . 

167070. 

0 0060 

166.0 

17063. 

16696. 

13229. 

2826. 

1775. 

15. 

133. 

81679. 

26950. 

31670 . 

197611 . 

211830. 

0.0066 

158.6 

5576 . 

6581 . 

10117. 

21960. 

16092. 

60. 

166. 

95771. 

31233. 

36711. 

229916. 

266222. 

0 .0072 

172.8 

26325. 

25595. 

29581. 

62208. 

26261. 

61 . 

186. 

106363. 

33508. 

39100. 

267052. 

262059. 

0 0078 

187.2 

60225. 

61230 . 

65062. 

56681. 

31616. 

75. 

191. 

106669. 

33822. 

38816. 

267676. 

263366. 

0 0086 

201.6 

53292. 

53968. 

56976. 

66278. 

35786 . 

83. 

188. 

103087. 

32375. 

36667. 

232925. 

268666. 

0.0090 

216.0 

63800. 

66217. 

66163. 

72303. 

37938. 

87. 

175. 

96597. 

29657. 

32519. 

206563 . 

220671. 

0 0096 

230.6 

72167 . 

72376 . 

73297. 

75987. 

38765. 

87. 

156. 

82393. 

25616. 

27660 . 

171578. 

182058. 

0.0102 

266.8 

78680 . 

78720 . 

78735. 

78116. 

38616 . 

86 . 

131. 

67666. 

20561. 

21512. 

130629. 

136366. 

0 0108 

259.2 

83658. 

83296. 

82687. 

78853. 

37825. 

82. 

103. 

50352. 

15006. 

16916. 

85336. 

85365. 

0 0116 

273.6 

86116. 

85760. 

86161. 

77806 . 

36272. 

77. 

70. 

31123. 

8850. 

7730. 

36583. 

30397. 

0 0120 

288.0 

85692. 

85105. 

82801. 

76001 . 

33521. 

69. 

32. 

9509. 

2069. 

530 . 

15385. 

28671. 

0.0126 

302.6 

81258. 

80689. 

77605. 

66936. 

29621. 

58. 

16. 

13110. 

5172. 

8695. 

73532. 

92661. 

0.0132 

316.8 

70862. 

69935. 

66617. 

56775. 

23123. 

62. 

51. 

37151. 

13035. 

18169. 

133129. 

156779. 

“ 0 0138 

331.2 

52996. 

51996. 

68372. 

36166. 

16027 . 

20 . 

100. 

63755. 

21596. 

27595. 

189179. 

215261. 

0.0166 

365.6 

27635. 

26398. 

22623. 

11007 . 

2133. 

9. 

150. 

90199. 

29781. 

35767. 

235039. 

261207. 

1G = 

1586 

1617 

. 

1668 

1687 

, 

709 

691 


696 

683 



TCA = 

1587. 

1587 . 

1617. 

1676. 

1666. 

1609. 

785. 

708. 

690. 

687. 

693. 

683. 


111 = 

1663. 

1672. 

1672. 

1672. 

1672. 

1609. 

786 . 

666 . 

662. 

662. 

666 . 

592. 



QIM = 383.616 QOUT= 259.727 WRKEXP= 363.322 WRKCNP=-202 . 968 WRKTOT= 160.356 EFFTOT= 0.618 REFF1= 0.997 REFF2= 0.996 



0.000 QCOMPP= 


QREGEN= 0.495 PWRHP= 19.470 FREQ= 66.780 WRKLT= 24.848 UTOTAL= 0.000 

El,3FM = 6310 7 1 0 AVGWSP= 2171.734 

QEXPN= -4.352 QEXPP= 0.134 QHEATN=-373 . 283 QHCATP= 2.760 QCOOLN= -4.789 QCOOLP= 201.925 QCOMPN- 
QREGN=-1774 . 719 QREGP=1775 . 215 QCNDRI = 5.178 QCNL’RG = 5.178 

QCNDCL = 1.388 QCNDD= 0.368 QCNDCN= 0.000 T GEXP A= 1 5 38 . 6 26 TGCMPA= 655.723 QSHTL= 1.939 

QINB = 374 741 WRKBAS = 185.202 WALT = 182.193 QUUTB= 200.008 WPKLEX= 19.589 WRKLCM= 5.260 WLEXR= 3.808 

Ml CMR = 2.608 WLEXE= 8.989 WLCME= 1.275 IJRKLP= 6.416 URKLH= 6.792 WRKLC= 1.377 WRKLE= 10.264 WLALT- 21.839 


x CYCLE NO. 25 x 


TIME 

0 . 

.0150 

0 

0156 

0 

0162 

0 . 

.0163 

0 

0174 

0 

0180 

0 

0186 

0 

0192 

0 

0198 

0 

0204 

0 

0210 

0 

0216 

0 

0222 

0 

0228 

0 

0234 

0 

0 24 0 

0 

0246 

0 

0252 

0 

0 258 

0 

0264 

0 . 

.0270 

0 

0276 

0 . 

. 0282 

0 

0288 

0 

0294 


ANGLE FOEXP FIHTR FOHTR FIREG FOREG FICLR FOCLR 
36 0.0-0.00746-0 .0106 2-0.0 3578-0 .05327-0 3033-0 . 24378-0 . 29628 

14.4-0.06709-0.06869-0.09021-0 . 1 1 06 0-0 . 2o 1 32-0 . 27237-0 . 3 1 49 9 
28.8-0 . 12272-0 .12343-0 .14399-0.16240-0 26890-0.27658-0.30599 
43.2-0 . 16625-0 .16688-0 . 18588-0.19812-0 .24990-0.25399-0.26872 
57.6-0.19347-0 .19425-0.20692-0.21164-0.20730-0 .20803-0.20854 
72.0-0. 20185-0.20221-0.20437-0 .20240-0.14794-0.14580-0.13407 
86 4-0.18996-0.18935-0.18051-0 
100.8-0 16030-0.15863-0.14092-0 
11951-0.11706-0.09387-0 
07377-0.07090-0.04510-0 
00022 0 


TG- 
1GA = 


115.2-0 
129.6-0 
144 .0-0 

158.4 0 

172.8 

187.2 
201.6 
216 . 0 
230 4 

244.8 

259.2 

273.6 
288 . 0 

302.4 

316.8 

331.2 

345.6 
1586 

1587 . 


02856-0 .02558-0 
00922 0.01215 0 


03980 
06517 
08560 
10178 
11457 
12471 
13243 
13723 
13756 
13176 
11625 
08811 
04624 
1617 
1587. 


04292 
. 06774 
. 08746 
.10305 
.11527 
.12482 
.13190 
.13602 
. 13566 
. 12926 
. 11326 
. 08478 
. 04274 


1617. 


03931 
07132 
09421 
10948 
11905 
12472 
12773 
12844 
12632 
11981 
10812 
08800 
05689 
01424 
1668 
1674. 


. 17352-0 
, 13060-0 
.08194 0 
.03239 0 
.01133 0 
.05149 0 
. 03421 
.10744 
.12198 
. 12957 
.13231 
.13186 
. 12396 
.12331 
. 11346 
.09899 


07674-0 .07475-0 
04431-0 . 12968-0 
00136-0 . 18351-0 
1487 . 

1666. 1409. 


07517-0.05401 
00291 0.02461 
07035 0.09881 
12896 0.16362 
17189 0.21519 
20171 0.25001 
21962 0.26645 
22440 0.26502 
0.21682 0.24916 
0.19887 0.22232 
0.17302 0.18733 
0.14132 0.14631 
0.10490 0.10051 
0.06392 0.05046 
0 01781-0.00342 
02472-0 .03107-0 . 06634 
08396-0.13317 
14189-0 .19867 


07952-0 
00390-0 
06 C 3 5 0 
0.11694 
0 . 16000 
0.19115 
0.21058 
0.21681 
0.21073 
0 19441 
17031 
14047 
10597 
06687 
02236 


FICMP 
-0.34836 
-0.35756 
-0.33471 
-0.28162 
-0.20596 
-0.11802 
-0 . 02739 
0 . 05841 
0.13494 
0 .19968 
0.25066 
0.28590 
0.30362 
0.30239 
0.28348 
0.25009 
0.20596 
0.15417 
0.09669 
0 .03457 
-0 .03257 
-0 .10687 
-0.18268 
-0.25361 


.19747-0.25515-0.31130 
709. 691. 


785. 


708. 

TM= 1643. 1672. 1672. 1672. 1672. 1409. 786. 646. 

Q I N = 383.612 QOUT - 259.746 WRKEXP= 363.323 WRKCMP=-202 . 96 9 WRKT0T= 160. 
OREGEN= 0.482 PWRHP= 19.470 FREQ = 66 780 WRKLT= 24.848 UTOT AL = 

EN 3PM- 6310.710 AVGMSP= 2171.742 

QEXPN= -4.352 QEXPP= 0.134 QHEATN=-373 . 28 1 QHEATP= 2.760 QCOOLN= 
QREGN=-1774 .724 QREGP=1775 . 207 QCNDRI= 5.178 QCNDRO= 5.178 
QCNDCL = 1.388 QCNDD= 0.368 QCNDCN= 0.000 TGEXP A = 1538 . 627 TGCMPA = 

Q1NB = 374.739 WRKBAS = 185.203 WALT = 182.193 QOUTB= 200.026 WRKLEX= 

8 989 WLCME= 1.275 WRKLR= 6.416 WRKLH= 6. 


690. 

642. 

354 EFFTOT= 

0.000 


694. 683. 

687. 693. 683. 

642. 646. 592. 

0.418 REFF1= 0 


997 REFF2= 0.996 


WLCMR : 


2.608 WLEXE= 


-4.787 QCOOLP= 201.942 QCOMPN= 0.000 QCOMPP= 

655.727 QSHTL= 1.939 

19.589 WRKLCM= 5.260 WLEXR= 3.808 

792 WRKLC= 1.377 WRKLE= 10.264 WLALT= 21.839 


xx AVERAGE VALUES OVER LAST 5 CYCLES xx 

QIN= 383.615 QOUT= 259.694 WRKEXP= 363.319 WRKCMP=-202 . 965 WRKTOT= 160.354 EFFTOT= 
QREGEN= 0.480 PWRHP= 19.470 FREQ= 66.780 WRKLT= 24.848 UTOTAL= 0.000 

EN3PM- 6310.710 AVGWSP= 2171.719 D 

QEXPN= -4.352 QEXPP= 0.134 QHEATN=-373 . 284 QHEATP= 2.760 QCOOLN= -4.792 QCOOLP 


0.418 REFF1= 


0.997 REFF2= 


0.996 


201.895 QCOMPN= 0.000 QCOMPP= 


.871 


.871 


.871 



t 


QREGN — 1774.741 QREGP =1 7 7 5 . 22 1 QCNDRI = 
<5 C N D C L = 1.388 QCNDD= 0.368 QCHDCN = 
QIND = 374 742 WRKBAS = 185.202 WALT = 
WLCMR= 2.608 WLF.XE= 8.989 WLCME= 


5.178 QCNDRO= 5.178 

0.000 TGEXPA=1538.626 TGCMPA= 655.714 QSHTL= 1.939 
182.192 Q0UTB= 199.974 WRKLEX= 19.589 WRKLCM= 5.260 WLEXR= 3.808 
1.275 WRKLR= 6.416 WRKLH= 6.792 WRKIC= 1.377 WRKLE= 10.264 WLALT= 


21.838 


** AVG. TEMPS.— LAST 5 CYCLES, AVG. H. T. COEFS. 8 HEAT FLUXES— LAST CYCLE *X 


TGCYC = 1530. 1528. 1524. 1482. 703. 691. 669. 659. 


TC-ACYC = 

1539 . 

1527. 

1545. 

1543. 

1505. 

1400 . 

780. 

697. 

679. 

668. 

665. 

656 . 

1MCYC = 

1643. 

1672. 

1672. 

1672. 

1672. 

1400 . 

780 . 

646. 

642. 

642. 

646. 

592. 

HACYC = 

0.057 

0.000 

0.808 

0.310 

0.000 

1.923 

1.879 

0.000 

0.622 

0.672 

0.000 

0.075 

!IMX = 

0.443 

0.000 

1.413 

1.435 

0.000 

2.950 

2.779 

0.000 

1.027 

1.099 

0.000 

0.245 

0 0 A A V G = 

5.65 

0.00 

69.15 

72.89 

0.00 

4.45 

4.38 

0.00 

22.67 

19.22 

0.00 

5.03 

OOAMX= 

22.90 

0.00 

132.68 

128.62 

0.00 

8.45 

8.53 

0.00 

46.29 

59.10 

0.00 

25.85 


** PPESSURE CALCULATIONS— EXPANSION AND COMPRESSION SPACE ** 

AVGPE-2172.765 AVGPC=21 7 1 . 588 PEMAX=2758 . 990 PEMIN=1663. 097 
PCMAX=2S00 .018 PCMIN= 16 32 . 7 0 1 

APFMAX= 53.280 APEMIN= 246.240 APCMAX= 54.000 APCNIN= 246.960 


-* ENGINE POWER AND EFFICIENCY CALCULATIONS XX 

BASICP= 22.487 FRLOSS= 4.294 BRKP= 15.176 BRKEFF= 0.326 
^ AUXLOS= 2.692 AUXPWR= 12.485 AUXEFF= 0.268 


COOLER AND APPENDIX GAP PUMPING CALCULATIONS XX 
8CHKPOM 

TH20AV= 585.98948208376 
RH20= 3.0311088482520 
RTUBE= 0.27440452269027 
QBTUPS= 22.310718811145 
U= 0.45779056053055D-02 
QHGPS= 15.373266780887 
0CGPS= 0.11240588565159 
? END 




TABLE VIII. - SHORT FORM OF OUTPUT DATA 
(PRODUCED WHEN IOUT = 1, JIP = 1) 


RUN IDENTIFICATION x# 
READ #123BR, 


“X CALCULATED CONTROL VOLUME AND ENGINE PARAMETERS, AND INPUT DATA xx 
ICPARAM 

UDEDV = 2.1788355358357 
RDEDV= 7.0 A 20962338519 
CDEDV- I . 4725243G37283 
DGAPDV= 0.38090003924925 

DHO= 2X2.1650, 3x0.11810, 2 2640, 5X0 . 27 1 90 952380 952D-02 , 0.22440D01, 3X0 . 39370D-01 
2-'0.169260D01, 13X0.0 

ACSO= 3.6813379069113, 1.8406689534556, 3X0 . 98589845512928D-01 , 0 . 39549007266381D01 
5-0.22938424214501D01, 0 . 39549007266381D01 , 3X0 . 23373402598861D0 , 0 . 17530335352219D01 
0.35060670704439D01, 13X0.0 

AHIO= 2X0.0, 3X11.033826006996, 0.0, 5X1035.9469775373, 0.0, 3X20.945228440128 
15X0.0 

XLO- 0.0, 0.73098424215499D-01, 3X0 . 36833333333333D0 1 , 0 . 128574655S8833D0 

5X0.3070D0, 0.39343591850981D-01, 3X0.1050D01, 0 . 52566022354121DO , 14X0.0 

VO= 0.0, 0.26910, 3X0.72627852861190, 1.0170, 5X1.4084192467704, 0.31120 

3X0.49084145457609, 1.8430, 14X0.0 

VCL E= 0.52610003924925 

VCLC= 0.17599000392493 

CLRLOD= 80.010160020320 

IITRLOD= 93.564775613887 

RGAREA= 3.9549007266381 

AP= 3.6813379069113 

AR= 0.17527083646743 

APMAR= 3.5060670704439 

VHAPGP= 0.38090003924925 

VCAPGP= 0.38090003924925D-01 

7HGP= 1373.6666666667 

TCGP= 646.40894661755 

QUGPS= 0.0 

QCGPS= 0.0 

I PUMP = 1 

ICOND= 0 

t END 


xx INPUT DATA— ENGINE PARAMETERS xx 
J ENGINE 

EID= -0 . 47377992506120D28 
ETYPE= 1.0 
DISPD= 2.1650 
DISPRD= 0.47240 
DSPGAP= 0.15750D-01 
DSPHGT = 3.530 
RODL= 3.9370 



il 


RCPAh'K = 0.78740 
E = 0.0 


PIIASE = 

90 . 0 


HTCOD= 

0 .17720 


HTBID- 

0.11810 


HTBPCN 

= 18.0 


HTDL = 

11.050 


EM 1 BL = 

9.9130 


REG°CN 

= 2.0 


REOID= 

2.2440 


REGL = 

1 .5350 


RUIRED 

= 0 19690D 

-02 

FROSTY 

= 0.580 


Rf1DEN = 

0 282 0 


crn= 0 

.110 


CTBOD= 

0 .59060D- 

01 

CTBID= 

0 . 39370D- 

01 

CTBPCN 

= 384.0 


CTBL = 

3.150 


ECT BL = 

2.6460 


CNDS S= 

0 .27780D- 

02 

CTH20= 

1.0 


RII0H20 

= 62.40 


AEf H2U 

= 5.0 


EXPSCL 

= 0.14520 


EXPMDV 

= 0.26910 


HRDV = 

1.0170 


RCDV = 

0 31120 


CCf'PDV 

= 1.8430 



ro CJ1I'SCL= 0.13790 
° CY10RT= 1.360 
CYL0RM= 1.310 
CYL0RB= 1.250 
C Y l D FM = 0.7870 
CYLD!1B = 0 7870 
D5PW f H = 0 750D-01 
REG0RT= 1.420 
REG0RM= 1.370 
REG0RB= 1.30 
REGDTM= 0.5940 
PECDMB= 0.5940 
STRUKE= 1.5750 
CAN0R= 0.0 
CANIR= 0.0 
CUHDTB = 0 0 
DFTEQ = 66.670 
DFL0SS= 17.1650 
DAL0SS= 10.060 
” i END 


v* INPUT DATA— OPTION SWITCHES AND MULTIPLYING FACTORS xx 


2 STRLNG 



F.EALGS- 1.0 
TACT1= 0.90 
F ACT2 = 10 0 
NCCYC= 25 
NS I RT = 1 
HO END= 20 
f.WGAS = 2 
RIICFAC = 1.0 
HHCFAC= 1.0 
CHCF AC= 1.0 
IPCV= 0 
H1UL T= 1.0 
FMUL TR= 1.0 

imx= 0 

VH2 = 0.900 
IPUMP= 1 
ICOND= 0 
1 0 U T = 1 
JIP= 1 
IFRINT= 500 
ITMPS= 0 
MAPLOT= 1 
& END 


'5* INPUT DATA — ENGINE OPERATING CONDITIONS XX 
2 INDATA 

I DRUN = -691351228, 1081869690, -205339165 

P = 2635.5939539073 

Of1EGA= 66.780 

TMEXP: 1693.0 

TMHFR= 1672.0 

TKUBR= 1672.0 

ICYL = 1693 0, 1991.0, 1239.0 

I CAN - 1191.0, 999.0 

TRO= 1972.0 

TR1= 1219.0 

TR2= 956.0 

Gt MH20= 13.630 

TU20IN= 580.10 

X END 

xx METAL TEMPERATURES FOR CONDUCTION CALCULATIONS *x 
SIN TEMP 

TCYL = 1693.0, 1991.0, 1239.0 

TCAM= 1191.0, 999.0 

TRO= 1972.0 

TR1= 1219.0 

TP,2= 956.0 

SEND 


QIN= 378.131 QOUT= 260.212 URKEXP= 358.513 WRKCMP=-202 . 013 WRKTOT= 156.500 EFFTOT= 


0.919 REFF1= 


0.997 REFF2= 


0.995 



QP.EGEN = -0.653 PWRHP= 19.002 FREQ= 66.780 WRKLT = 24.639 UTOT AL = 

tN3m= 6310 710 AVGWSP= 2157.025 

QEXPN= -4.589 QEXPP= 0.125 QHEATN=-386 .310 QHEATP= 2.176 QCOOLN= 
QREGN=-1761 219 QREGP = 1760 . 561 QCNDRI = 5.178 QCHDRO= 5.178 
QCNDCL = 1.388 QCNDD= 0.365 QCNDCN= 0 . 000 TCEXPA = 1532 . 013 TGCMPA = 

QINB = 388 598 WRKBAS = 181.134 WALT = 178.041 QOUTB= 195.212 WRKLEX = 
ULCKR= 2.635 WLEXE= 8.881 ULCnE= 1.278 WRKLR= 6.388 WRKLH= 6 


0.000 

-4.616 QCOOLP= 197.033 QCOMPN= 0.000 QCOMPP= 

654.260 QSHTL= 1.939 

19.338 WRKLCM= 5.296 WLEXR= 3.753 

704 WRKLC= 1.383 WRKLE= 10.159 WLALT= 21.541 


2.795 


xx LAST CYCLE xx 


0.414 REFF1= 0.997 REFF2= 0.995 


QIN= 378.131 QOUT= 260.212 WRKEXP= 358.513 WRKCMP=-2 02 . 0 1 3 WRKTOT= 156.500 EFFTOT= 

QREGEN= -0.653 PWRHP= 19.002 FREQ= 66.780 WRKLT= 24.634 UTOTAL= 0.000 

QEXPN= -4l589 1 QEXPP= SP 0 2 125 * QHEATN=-386 . 31 0 QHEATP= 2.176 QCOOLN= -4.616 QCOOLP= 197.033 QCOMPN= 0.000 QCOMPP= 
4REGN=-1761 .219 QREGP=176 0 . 56 1 QCNDRI= 5.178 QCNDRO= 5.178 

QCNDCL= 1.388 QCNDD= 0.365 QCNDCN= 0.000 TGEXP A= 1 532 . 0 1 3 TGCMPA= 654.260 QSHTL- 

QINB = 388.598 WRKBAS = 181.134 WALT = 178.041 Q0U1B= 195.212 WRKLEX= 19.338 WRKLCM= - - - , . ... 

WLCMR= 2.635 WLEXE= 8.881 WLCME= 1.278 WRKLR= 6.388 WRKLH= 6.704 WRKLC= 1.383 WRKLE= 10.159 WLALT- 21.541 


2.795 


1.939 

5.296 WLEXR= 


3.753 


AVERAGE VALUES OVER LAST 5 CYCLES XX 


QIM= 378 176 QOUT= 260.410 WRKEXP= 358.522 WRKCMP=-202 . 027 WRKTOT= 156.495 EFFTOT= 

QREGEH= -0.648 PWRHP= 19.001 FREQ= 66.780 WRKLT= 24.635 UTOTAL= 0.000 

EH 3PM- 6310.710 AVGWSP= 2157.094 , 

QEXPH= -4.590 QEXPP= 0.125 QHEATN=-386 . 352 QHEATP= 2.172 QCOOLN= -4.599 QCOOLP= 197.216 QCOMPN- 
0REGN=-1761 .119 QREGP=1760 .471 QCNDRI= 5.178 QCNDRO= 5.178 

QCHDCL= 1.388 QCNDD= 0.365 QCNDCN= 0.000 TGEXPA=1531 . 996 TGCMPA= 654.312 QSHTL= 1.939 
M QIHB = 388 644 WRKBAS = 181.130 WALT = 178.036 QOUTB= 195.414 WRKLEX= 19.338 WRKLCM= 5.296 WLEXR= 
FS WLCMR= 2.635 WLEXE= 8.881 WLCME= 1.278 WRKLR= 6.388 WRKLH= 6.704 WRKLC= 1.383 WRKLE= 10 


0.414 REFF1 = 0.997 REFF2= 0.995 


0.000 QCOMPP= 2.797 


3.753 

159 WLALT= 21.541 


XX AVO 

TEMPS.-- 

LAST 5 

CYCLES, 

AVG. H. T. COEFS. 

i HEAT 

FLUXES- 

-LAST CYCLE xx 



TGCYC- 

1523. 1522. 

1519. 1477. 


701. 689. 

668. 

658. 

TGACYC= 

1532. 

1521. 

1540 . 

1538. 1500. 

1395. 

778. 

695. 678. 

667. 663. 

654 . 

TI'CYC = 

1643. 

1672. 

1672. 

1672. 1672. 

1395. 

778. 

646. 641. 

641 . 646 . 

592. 

HACYC= 

0 . 056 

0.000 

0.805 

0.807 0.000 

1 . 915 

1.878 

0.000 0.622 

0.671 0.000 

0 . 075 

HMX- 

0 447 

0.000 

1.415 

1.439 0.000 

2.957 

2.778 

0.000 1.025 

1.096 0.000 

0.242 

QOAAVG= 

5.93 

0.00 

71.82 

74.96 0.00 

4.41 

4.36 

0.00 22.19 

18.66 0.00 

4.90 

QOAt1X = 

25.24 

0.00 

137.42 

131.73 0.00 

8.43 

8.53 

0.00 45.16 

57.51 0.00 

25.05 


PRESSURE CAL CLU. AT IONS — EXPANSION AND COMPRESSION SPACE XX 

AVGPE=2158.025 AVGPC=2156 .872 PEMAX=2735.649 PEMIN=1653 . 789 
“ PCMAX=2777 . 046 PCMIN=1623 . 724 

APEMAX= 52.560 APEMIN= 246.240 APCMAX= 53.280 APCMIN= 246.240 


xx ENGINE POWER AND EFFICIENCY CALCULATIONS XX 



3 A S I CP = 21.993 FRLOSS= 4.272 BRKP= 14.730 BRKEFF= 

AUXLOS= 2.686 AUXPWR= 12.045 AUXEFF= 0.262 


0.321 


xv COOLER AND APPENDIX GAP PUMPING CALCULATIONS xx 
iCHKROM 

TU20AV= 586.01304542543 
RH20 = 3.0307937576335 
RTUBE= 0.27440452269027 
9BTUPS= 22.399982193349 
U' = 0.45617150644983D-02 
CHGPS= 15.544561636230 
QCGPS= 0.10697760560297 
SEND 


QIN= 383.612 QOUT = 259.746 WRKEXP= 363.323 WRKCMP=-202. 969 WRKTOT= 160.354 EFFTOT= 0.418 REFF1= 0.997 REFF2= 0.996 
QRCGCN= 0.482 PWRHP= 19.470 FREQ= 66.780 WRKLT= 24.848 UTOTAL= 0.000 
EN3PM= 6310.710 AVGWSP= 2171.742 

QEXPN= -4.352 QEXPP= 0.134 QHEATN=-37 3.281 QHEATP= 2.760 qCOOLN= -4.787 QCOOLP= 201.942 QCOMPN= 0.000 QCOMPP= 2.871 
QREGN =-1774. 724 QREGP=1775 . 207 QCNDRI= 5.178 QCNDRO= 5.178 

QCNDCL= 1.388 QCNDD= 0.368 QCNDCN= 0.000 TGEXPA=1538 . 627 TGCMPA= 655.727 QSHTL= 1.939 

QINB = 374.739 WRKBAS = 185.203 WALT = 182.193 QOUTB= 200.026 WRKLEX= 19.589 WRKLCM= 5.260 WLEXR= 3.808 

WLCMR= 2.608 WLEXE= 8.989 WLCME= 1.275 WRKLR= 6.416 WRKLH= 6.792 WRKLC= 1.377 WRKLE= 10.264 WLALT= 21.839 


xv LAST CYCLE xx 

co QIN= 383.612 QOUT= 259.746 WRKEXP= 363.323 WRKCNP=-202 . 96 9 WRKTOT= 160.354 EFFTOT= 0.418 REFF1= 0.997 REFF2= 0.996 
QREGEN= 0.482 PWRHP= 19.470 FREQ= 66.780 WRKLT= 24.848 UTOTAL= 0.000 
EN3PM= 6310.710 AVGWSP= 2171.742 

QEXPN= -4.352 QEXPP= 0.134 QHEATN=-373 . 281 QHEATP= 2.760 QCOOLN= -4.787 QCOOLP= 201.942 QCOMPN= 0.000 QCOMPP= 2.871 
QREGN =-177 4. 724 QREGP=1775 . 207 QCNDRI= 5.178 QCNDRO= 5.178 

CCNDCL = 1.388 QCNDD= 0.368 QCNDCN= 0 . 000 TGEXPA = 1538 . 627 TGCMPA= 655.727 QSHTL= 1.939 

QINB = 374.739 WRKBAS = 185.203 WALT = 182.193 QOUTB= 200.026 WRKLEX= 19.589 WRKLCM= 5.260 WLEXR= 3.808 

WLCMR= 2.608 WLEXE= 8.989 WLCME= 1.275 WRKLR= 6.416 WRKLH= 6.792 WRKLC= 1.377 WRKLE= 10.264 WLALT= 21.839 


xx AVERAGE VALUES OVER LAST 5 CYCLES #x 

QIN= 383.615 QOUT= 259.694 WRKEXP= 363.319 URKCMP=-202 . 965 WRKTOT= 160.354 EFFTOT= 0.418 REFF1= 0.997 REFF2= 0.996 
QREGEN= 0.480 PWRUP= 19.470 FREQ= 66.780 WRKLT- 24.848 UTOTAL= 0.000 

EN3PM= 6310.710 AVGWSP= 2171.719 

QEXPN= -4.352 QEXPP= 0.134 QHEATN=-37 3 . 284 QHEATP= 2.760 QCOOLN= -4.792 QCOOLP= 201.895 QCOMPN= 0.000 QCOMPP= 2.871 
QREGN=-1774 .741 QREGP=1775 . 221 QCNDRI= 5.178 QCNDRO= 5.178 

QCNDCL = 1.388 QCNDD= 0.368 QCNDCN= 0.000 TGEXPA = 1538 . 626 TGCMPA= 655.714 QSHTL= 1.939 

“ QINB = 374.742 WRKBAS = 185.202 WALT = 182.192 QOUTB= 199.974 WRKL EX= 19.589 WRKLCM= 5.260 WLEXR= 3.808 

ULCMR= 2.608 WLEXE= 8.989 WLCME= 1.275 WRKLR= 6.416 WRKLH= 6.792 WRKLC= 1.377 WRKLE= 10.264 WLALT= 21.838 


v* AVG. TEMPS.— LAST 5 CYCLES, AVG. H. T. COEFS. i HEAT FLUXES— LAST CYCLE XX 



E 


TGCYC = 1530. 1528. 1524. 1482. 703. 691. 669. 659. 


TGACYC= 

1539. 

1527 . 

1545. 

1543. 

1505. 

1400 . 

780 . 

697 . 

679. 

668. 

665. 

656. 

1MCYC = 

1643. 

1672. 

1672. 

1672. 

1672. 

1400 . 

780 . 

646 . 

642. 

642. 

646 . 

592. 

HACYC= 

0.057 

0.000 

0.808 

0.810 

0.000 

1 . 923 

1.879 

0.000 

0.622 

0.672 

0.000 

0.075 

UMX = 

0.443 

0.000 

1.413 

1.435 

0.000 

2.950 

2.779 

0.000 

1.027 

1.099 

0.000 

0.245 

00 A A VG = 

5.65 

0 00 

69.15 

72.89 

0.00 

4.45 

4.38 

0.00 

22.67 

19 22 

0.00 

5.03 

QOAMX= 

22.90 

0.00 

132.68 

128.62 

0 .00 

8.45 

8.53 

0.00 

46.29 

59.10 

0.00 

25.85 


PRESSURE CALCUL AT IONS- - EXPANSION AND COriPRESSION SPACE X* 

AVGPE=21 72 . 765 AVGPC=21 7 1 . 588 PENAX=2758 . 990 PEMIN=1663 . 097 
rC!'.AX = 280 0 . 018 PCNIN = 16 32 . 7 0 1 

APEH/'X= 53.280 APEMIN= 246.240 APCMAX= 54.000 APCMIN= 246.960 


“X ENGINE POWER AND EFFICIENCY CALCULATIONS *x 

BASICP = 22.487 FRLOSS= 4.294 BRKP = 15.176 BRKEFF= 0.326 

AUXLOS = 2.692 AUXPWR= 12.485 AUXEFF= 0.268 


** COOLER AND APPENDIX GAP PUMPING CALCULATIONS *x 
s rnKJ?nM 

T H 2 0 A V = 585 98948208376 
RH20= 3.0311088482520 
£ RTUBE= 0.27440452269027 
QBTUPS= 22.310718811145 
W= 0 .45779056053055D-02 
QHGFS= 15.373266780887 
QCGPS= 0.11240588565159 
ifcND 


n 



TABLE IX. - FINAL SUMMARY PRINTOUT 


ENGINE OPERATING CONDITIONS 


ENGINE SPEED 

4006.8 

RPM 

66 

.78 

HZ 


MEAN PRESSURE 

2171.7 

PSI 

14 

. 98 

MPA 


HEATER TUBE OUTSIDE TEMPERATURE 

1672.0 

R 

1212.3 F 


928.9 K 

655.7 C 

COOLANT INLET TEMPERATURE 

580 . 1 

R 

120.4 F 


322.3 K 

49.1 C 

RESULTANT INSIDE COOLER TUBE TEMP. 

641.9 

R 

182.2 F 


356.6 K 

83.5 C 

ENGINE PERFORMANCE SUMMARY 







BPAKE POWER-1 CYLINDER 

12.485 

HP 

9.310 

KW 

102.824 

FT-LBF/CYCLE 

-4 CYLINDERS 

49.939 

HP 

37.239 

KW 

411.296 

FT-LBF/CYCLE 

BRAKE EFFICIENCY 

0.268 






INDICATED POWER-1 CYLINDER 

19.470 

HP 

14.519 

KW 

160.355 

FT-LBF/CYCLE 

-4 CYLINDERS 

77.880 

HP 

58.075 

KW 

641.418 

FT-LBF/CYCLE 

INDICATED EFFICIENCY 

0.418 






HEAT RATE TO ENGINE-1 CYLINDER 

46.577 

HP 

34.733 

KW 

383.612 

FT-LBF/CYCLE 

-4 CYLINDERS 

186 . 310 

HP 

138.931 

KW 

1534.447 

FT-LBF/CYCLE 

HEAT RATE FROM ENGINE 







HEAT RATE TO COOLANT-1 CYLINDER 

31.538 

HP 

23.518 

KW 

259.746 

FT-LBF/CYCLE 

-4 CYLINDERS 

126.151 

HP 

94.071 

KW 

1038.982 

FT-LBF/CYCLE 

AUXILIARY LOSS (4 CYLINDERS) 

10.767 

HP 

8.029 

KW 

88.677 

FT-LBF/CYCLE 

TOTAL RATE (4 CYLINDERS) 

136.918 

HP 

102.100 

KW 

1127.659 

FT-LBF/CYCLE 

'< ERROR IN ENERGY BALANCE 

0.294 

% 





£ MECHANICAL LOSS (4 CYLINDERS) 

17.174 

HP 

12.807 

KW 

141.445 

FT-LBF/CYCLE 

FLOW FRICTION LOSS (4 CYLINDERS) 
INDICATED WORK PER CYCLE SUMMARY (1 

12.068 

HP 

8.999 

KW 

99.393 

FT-LBF/CYCLE 

OIL INUtK ) •• 






EXPANSION SPACE (WRKEXP) 

363.323 

FT- 

LBF/CYCLE 




COMPRESSION SPACE (WRKCMP) 

-202.969 

FT- 

LBF/CYCLE 




NET (WRKTOT ) 



160.: 

355 

FT-LBF/CYCLE 



HEAT FLOW SUMMARY (1 CYLINDER) 

HEAT RATE TO ENGINE 
EXPANSION SPACE HEAT RATE 
METAL TO GAS (QEXPN) -4.352 

GAS TO METAL (QEXPP) 0.134 

NET (QEXP ) 

HEATER HEAT RATE 

METAL TO GAS CQHEATN) -373.281 

GAS TO METAL ( QHEATP ) 2.760 

NET ( QHEATR ) 

CONDUCTION LOSSES 

THROUGH REGENERATOR HOUS ING( QCNDRI ) 5.178 

THROUGH CYLINDER HOUSING! QCNDCL ) 1.388 

DIRECTLY THROUGH PISTON(QCNDD) 0.368 

SHUTTLE LOSS VIA PISTON ( QSHTL ) 1.939 

NET (QCNDTI) 

NET HEAT RATE TO ENGINE (QEIN) 

(- SIGN MEANS FLOW INTO ENGINE) 


FT-LBF/CYCLE 
FT-LBF/CYCL E 

-4.218 FT-LBF/CYCLE 

FT-LBF/CYCLE 

FT-LBG/CYCLE 

-370.521 FT-LBF/CYCLE 

FT-LBF/CYCLE 

FT-LBF/CYCLE 

FT-LBF/CYCLE 

FT-LBF/CYCLE 

-8.873 FT-LBF/CYCLE 

-383.612 FT-LBF/CYCLE 


HEAT RATE FROM ENGINE 



II 


COOLER HEAT RATE 

GAS TO METAL (QCOOLP) 201. 942 

METAL TO GAS (QCOOLN) -A. 787 

NET ( QCOOLR ) 

COMPRESSION SPACE HEAT RATE 
GAS TO METAL (QCOMPP) 2.871 

METAL TO GAS (QCOMPN) 0.000 

NET ( QCOMP ) 

APPENDIX GAP PUMPING LOSSES 
HOT GAP ( QHGPS ) 15.373 

COLD GAP ( QCGPS ) 0.112 

NET ( QAPGAP ) 

CONDUCTION LOSSES (QCNDTO) 


TOTAL HEAT FLOW TO COOL ANT , EXCLUDING MECHANICAL 
MECHANICAL LOSSES (1 CYLINDER) 


FT-LBF/CYCLE 

FT-LBF/CYCLE 

197.155 


FT-LBF/CYCLE 

FT-LBF/CYCLE 

2.871 


FT-LBF/CYCLE 

FT-LBF/CYCLE 

15.486 

8.873 

LOSSES 

35.361 


FT-LBF/CYCLE 


FT-LBF/CYCLE 


FT-LBF/CYCLE 

FT-LBF/CYCLE 

224.384 FT-LBF/CYCLE 

FT-LBF/CYCLE 


NET HEAT RATE TO COOLANT(QCLOUT) 
AUXILIARY LOSSES (1 CYLINDER) 

NET HEAT RATE FROM ENGINE ( QEOUT ) 


22.169 


259.746 FT-LBF/CYCLE 


FT-LBF/CYCLE 

281.915 FT-LBF/CYCLE 


REGENERATOR HEAT FLOW 

METAL TO GAS ( QREGN ) -1774.724 FT-LBF/CYCLE 

GAS TO METAL (QREGP ) 1775.206 FT-LBF/CYCLE 

NET ( QREGEN ) 

X ERROR REG. ENERGY BALANCE(PREGER) 

( QREGEN/ (MINIMUM OF ABS. VALUE OF QREGN 8 QREGP)) 


0.482 FT-LBF/CYCLE 
0.027 X 


REGENERATOR EFFECTIVENESS CALCULATION (BASED ON ENTHALPY FLOW PER CYLINDER) 


NET ENTHALPY FLOW REG. TO HTR . ( ENFRTH) 
NET ENTHALPY FLOW HTR. TO REG . ( ENFHTR) 
REG. EFFECT. (REFF1=ENFRTH/ENFHTR) 

NET ENTHALPY FLOW CLR . TO REG. (ENFCTR) 
NET ENTHALPY FLOW REG. TO CLR . ( ENFRTC) 
REG. EFFECT. (REFF2=ENFCTR/ENFRTC) 


PRESSURE DROP LOSS SUMMARY (PER CYLINDER) 
HEATER (WRKLH) 

REGENERATOR 
HOT SIDE (WLEXR) 

COLD SIDE (WLCMR) 

NET (WRKLR) 

COOLER (WRKLC) 

CONNECTING DUCTS (END EFFECTS) 

HOT SIDE (WLEXE) 

COLD SIDE (WLCME) 

‘ NET (WRKLE) 

NET HOT SIDE (WRKLEX) 

NET COLD SIDE (WRKLCM) 

NET ENGINE PRESSURE DROP LOSS (WRKLT) 


3284.703 FT-LBF/CYCLE 
3294.560 FT-LBF/CYCLE 
0.9970 

2108.668 FT-LBF/CYCLE 
2117.549 FT-LBF/CYCLE 


0.9958 



6.792 

FT-LBF/CYCLE 

3.808 

2.608 

FT-LBF/CYCLE 

FT-LBF/CYCLE 

6.416 

1.377 

FT-LBF/CYCLE 

FT-LBF/CYCLE 

8.989 

1.275 

FT-LBF/CYCLE 

FT-LBF/CYCLE 

10.264 

FT-LBF/CYCLE 


19.589 FT-LBF/CYCLE 
5.260 FT-LBF/CYCLE 

24.848 FT-LBF/CYCLE 


PRESSURES 

EXPANSION SPACE 
MAXIMUM (PEMAX) 

~Z MINIMUM (PEMIN) 

MEAN ( AVGPE) 

RATIO (PEMAX/PEMIN) 
COMPRESSION SPACE 


2759.0 PSI 

1663.1 PSI 
2172.8 PSI 

1.659 


19.028 MPA 
11.470 MPA 
14.985 MPA 



SYMBOL 

IDRUN 

P 

OMEGA 

TMEXP 

TMHFR 

TMHBR 

TCYL(l) 

T CYL ( 2 ) 

TC YL ( 3 ) 

TCAN(l) 

TCAN(2) 

TRO 

TRl 

TR2 

GPMH20 

TH20IN 


TABEL VI. - SYMBOL DEFINITIONS FOR ENGINE 
OPERATING CONDITIONS 
(NAMELIST /INDATA/) 

DEFINITION 

Alphanumeric run identifier 
Mean pressure, Ibf/in , (MPa) 

Engine frequency, hz 

Expansion space wall temperature, °R (K) 

Outside temperature of front row (flame side) portion of heater 
tubes, °R (K) 

Outside temperature of back row portion of heater tubes, *R (K) 
Cylinder housing temperature, top, °R (K) 

Cylinder housing temperature, middle, °R (K) 

Cylinder housing temperature, bottom, *R (K) 

Insulation container temperature, top, °R (K) 

Insulation container temperature, bottom, °R (K) 

Regenerator housing temperature, top, °R (K) 

Regenerator housing temperature, middle, °R (K) 

Regenerator housing temperature, bottom, °R (K) 

Cooling water flow rate per cylinder, gal. /min (liter/sec) 
Cooling water inlet temperature, °R (K) 
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TABLE V. - SYMBOL DEFINITIONS (AND TEST CASE SETTINGS) FOR MODEL 
OPTION SWITCHES AND MULTIPLYING FACTORS 
(NAMELIST /STRLNG/) 


SYMBOL 

REALGS 

FACT1 

FACT2 

NOCYC 

NSTRT 

NOEND 

MWGAS 

RHCFAC 

HHCFAC 

CHCFAC 

IPCV 


FMULT 
FMULTR 
I M I X 

VH2 

I PUMP 

ICOND 


I OUT 
JIP 
I PRINT 
ITMPS 


MAPLOT 


SETTING 

1. 

0. 

0. 4, 

10 . or 

25 

1 

20 

2 

4 

1. 

1. 

1. 

0 

1 

1.0 

1.0 

1 

0 

0.99 

1 

0 

1 


0 


1 

0 

0 

1 

500 

1 

0 

1 

0 


DEFINITION 

Use real gas equation of state 
Use ideal gas equation of state 

Empirical factors used in regenerator matrix temperature 
convergence procedure 

Number of engine cycles to be calculated (per pass) 

Cycle number at which regenerator matrix temperature 
convergence procedure begins 

Cycle number at which regenerator matrix temperature 
convergence procedure ends 
Use hydrogen working gas 
Use helium working gas 

Regenerator heat transfer coefficient multiplying factor 

Heater heat transfer coefficient multiplying factor 

Cooler heat transfer coefficient multiplying factor 

Make second pass through calculations to improve 

prediction of effect of pressure drop 

Eliminate second pass 

Overall pressure drop multiplying factor 

Regenerator pressure drop multiplying factor 

Use mixture of hydrogen and carbon dioxide working gas 

Pure hydrogen or helium working gas 

Volume fraction of hydrogen in hydrogen-carbon dioxide 

mixture (used only if IMIX=1) 

Calculate pumping loss due to piston-cylinder gap 
Omit pumping loss calculation 

Calculate cylinder and regenerator housing temperatures 
from TM(1), TM(4) and TH20IN (Input hot and cold end 
temperatures) 

Use the specified input values of the cylinder and 
regenerator housing temperatures for conduction 
calculations 

Write out Table VII or VIII data 

Don't write out Table VII or VIII data 

Write out Table VII data if I0UT=1 

Write out Table VIII data if I0UT=1 

Number of time steps between variable printouts in 

Table VII data 

Write out instantaneous gas temperatures at each time 
step in Table VII (for debugging) 

Don't write out instantaneous gas temperatures at each 
time step 

Store variables for plotting 
Don't store variables for plotting 
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I 


s 


MAXIMUM ( PCMAX) 
MINIMUM CPCMIN) 

MEAN ( AVGPC ) 

RATIO (PCMAX/PCMIN) 


2800.0 PSI 
1632.7 PSI 
2171.6 PSI 
1.715 


19.310 MPA 
11.260 MPA 
14.976 MPA 


> 


MEAN PRESSURE, CTR. OF REG. CAVGWSP) 2171.7 PSI 14.978 MPA 


AVERAGE PRESSURE RATIO ( ( PEMAX+PCMAX)/ ( PEMIN+PCMIN ) ) 1.687 


ro 




1-H CM CO ^ in 


TABLE X. - READ/WRITE UNIT NUMBERS 
USED FOR INPUT/OUTPUT 



INPUT/OUTPUT 

STATEMENT 

SUBROUTINE 

UNIT 

# 

INPUT/OUTPUT DATA 
(TEST CASE) 

1. 

READ 

ROMBC 

4 

TABLE II 

2. 

READ 

ROMBC 

5 

TABLE IV 

3. 

WRITE 

ROMBC, CYCL 

6 

TABLE VII (if IOUT=l and JIP=0) 

4. 

WRITE 

CYCL 

16 

TABLE IX 

5. 

WRITE 

ROMBC 

13 

BINARY OUTPUT (if MAPLOT=l) 
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TABLE XI. - ARRAYS SET VIA EQUIVALENCE 
STATEMENT IN SUBROUTINE ROMBC 


COMMON/RESET/ 
EQUIVALENCE ARRAY 
FOR COMMON/RESET/ 

COMMON/TIMEO/ 


COMMON/TIME1/ 


COMMON/TIME2/ 


COMMON/PSET/ 


COMMON /ASET/ 


WRKLT, WRKLR, WRKLE, WRKLTO, WRKLRO, WRKLHO, WRKLEO, 

SET ( 1 ) , SET(2) , SET (3) , SET(4), SET{5), SET(6), SET(7) 

VARIABLES 

WEPO, WCPO, WEPEO, WCPCO, WEDURO, WCDURO, WEDUHO, 
WO(l), WO ( 2 ) , WO(3), W0(4) , WO(5), WO(6), WO(7), 
WCDUCO, WEDUEO, WCDUCO, WEALTO, WTPCO 

W0(8), WO(9), WO(IO), WO(ll), WO ( 12) 

WEP1, WCP1, WEPE1, WCPC1, WEDUR1, WCDUR1, WEDUH1, 

Wl(l), Wl(2), Wl(3), Wl(4), Wl(5) , Wl(6) , Wl(7). 

WCDUC1, WEDUE1, WCDUC1, WEALT1, WTPC1 
WI(8), WI( 9) , WI(IO), WI(ll), WI( 12) 

WEP2, WCP2, WEPE2, WCPC2, WEDUR2, WCDUR2, WEDUH2, 

W2( 1) , W2(2) , W2(3) , W2(4) , W2(5), W2(6), W2(7), 

WCDUC2, WEDUE2, WCDUC2, WEALT2, WTPC2 
W2(8), W2(9) , W2( 10) , W2( 11) , W2(12) 

PRIN, PROUT, PEXP, PCMP, PEDUR, PCDUR, PEDUH, PCDUC, 

PS( 1) , PS ( 2 ) , PS(3) , PS ( 4 ) , PS( 5) , PS(6), PS(7), PS(8), 

PEDUE, PCDUE, PEALT 

P(9) P(10) P( 11) 

AEP, ACP, AE, AC, AEDUR, ACDUR, AEDUH, ACDUC, 

AS( 1) , AS ( 2 ) , AS (3) , AS (4) , AS(5), AS(6), AS(7), AS(8), 

AEDUE, ACDUE, AEALT, ATPC 

AS ( 9 ) , AS(10) , AS (11) , AS (12) 
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-Preheater-combustor assembly 


Insulation 


Fuel injector- 
Turbulator 


Drive shaft 


Crankshaft 


Lube oil pump 
Sump 


Ignitor 

Combustion 
chamber 


Heater head 


Piston assembly 
Regenerator 

Exhaust 



Simplified schematic of square- 
four cylinder arrangement 


Figure L - P40 Stirling engine cross section. 
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Figure 10 - Generalized control volume 


No 


26 

Reset 

time-0 

and 

make 

2nd 

pass 

using 

Zip 

Information 
from 1s t pass 


No 


Inputs, preliminary calculations, initializations 

*-L Update time and crank angle, ROMBC 
2 Update expansion and compression space volumes, ROMBC 

3. Update thermal conductivity and viscosity for gas 
control volumes, HEATX 

4. Updat pressure level, HEATX 

5 Update gas specific heats. HEATX 

6 Update gas temperatures for effect of change in specific 
heats, HEATX 

7. Update gas temperatures for effect of change in 
pressure, HEATX 

8 Update mass distribution, HEATX 

9 Update flow rates, HEATX 

10 Update gas temperatures for effect of 
flow between control volumes, HEATX 

1L Update heat transfer coefficients, HEATX 
12 Update gas temperatures for effect of 
metal-gas heat transfer, HEATX 

13. Update regenerator matrix temperatures, HEATX 

14. Update friction factors and pressure drops for 
each control volume, XDEL 

15 Update conduction and shuttle losses 
once each cycle, CNDCT 

16 Sum up heat transfers for each component, HEATX 

17 Calculate work and sum vp for cycle, ROMBC 
-E Is cycle complete 7 

J Yes 

19 Revise regenerator matrix temperatures?, ROMBC 

I Yes, I No 

| Make revision J 

20 Calculate indicated power and efficiency, CYCL 
21. Calculate mechanical friction losses, CYCL 

22 Calculate auxiliary losses and brake power and efficiency, CYCL 
(if specified number of cycles has been completed) 

23 Revise cooler tube temperatures?, ROMBC 

j Yes, I No 

| Make revision I 

24. Have specified no of cycles been completed?, ROMBC 
JYes 

•25 Is this the second pass through nocyc cycles?, ROMBC 
| Yes 

27. Write summary of predictions and terminate run 


Figure 1L - Outline of calculation procedures 




Regenerator node 

(a) Sample regenerator control volume 

(b) Control-volume temperature profile. 

Figure 12 - Sample regenerator control volume and temperature profile 



Engine speed, rpm 

Figure 13 - Mechanical power loss as a function of engine speed and 
mean pressure 
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Engine speed, rpm 

Figure 14 - Auxiliary power requirement as a function of engine speed 
and mean pressure 


Nominal operating conditions 
15 MPa mean pressure 
720° C heater set temperature 



Figure 15 - P-40 brake power and efficiency as functionsof engine 
speed. 


Nominal operating conditions- 


15 MPa mean pressure 



Engine speed, rpm 


Figure 16 - P-40 brake power and efficiency as functions of engine 
speed 
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